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PREFACE 


THis book is intended primarily for the student of applied 
mathematics, physics, chemistry or engineering who wishes 
to use the “ special ’ functions associated with the names of 
Legendre, Bessel, Hermite and Laguerre. It aims at provid- 
ing in a compact form most of the properties of these 
functions which arise most frequently in applications, and 
at establishing these properties in the simplest possible way. 
For that reason the methods it employs should be intelli- 
gible to anyone who has completed a first course in calculus 
and has a slight acquaintance with the theory of differential 
equations. Use is made of the theory of functions of a 
complex variable only very sparingly, and most of the book 
should be accessible to a reader who has no knowledge of 
this theory. Throughout the text an attempt is made to 
show how these functions may be used in the discussion 
of problems in classical physics and in quantum theory. 
A brief account is given in an appendix of the main 
properties of the Dirac delta “ὁ function ”’. 

I should like to record my debt of gratitude to the late 
Sir John Lennard-Jones, and to my colleagues Mr. B. Noble 
and Dr. J. G. Clunie for their generous help in reading the 
first draft of the manuscript and making valuable sugges- 
tions for its improvement. I am indebted to Miss Janet 
Burchnall for her assistance in the preparation of the final 
manuscript, to Mr. J. S. Lowndes for help in correcting 
proof sheets and to Miss Elizabeth Gildart for preparing 
the index. I should also like to thank Dr. D. E. Rutherford, 
general editor of the series, for his advice and criticism 
throughout the preparation of the book. 
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My debt of gratitude to Professor T. M. MacRobert is 
a much more general one. It was at his lectures that I first 
acquired a taste for the subject, and it will be obvious to 
anyone who knows his published writings how much I have 
been influenced by them. 


KEELE, STAFFORDSHIRE 
20th August 1955 
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CHAPTER I! 


INTRODUCTION 


81, The origin of special functions. The special functions 
of mathematical physics arise in the solution of partial 


differential equations governing the behaviour of certain 


physical quantities. Probably the most frequently occur- 
ring equation of this type in all physics is Laplace’s equation 

V7 =0 (1.1) 
satisfied by a certain function y describing the physical 
situation under discussion. The mathematical problem 
consists of finding those functions which satisfy equation 
(1.1) and also satisfy certain prescribed conditions on the 
surfaces bounding the region being considered. For 
example, if y denotes the electrostatic potential of a system, 
w will be constant over any conducting surface. The shape 
of these boundaries often makes it desirable to work in 
curvilinear coordinates g,, 4.5, 83 instead of in rectangular 
Cartesian coordinates x,y,z. In this case we have relations 


Χ = X(G1> ἢ» 43). VY =V(41> 42: 3)» 2 = 2(G1> Ta, 3) (1.2) 
expressing the Cartesian coordinates in terms of the curvi- 
linear coordinates. Jf equations (1.2) are such that 

Ox Ox - ὃν ὃν Jj. Oz 0z = 

0q;0q; 6q;0q; 94; δῆ; 


when i #7 we say that the coordinates q;, 712, q3 are 
orthogonal curvilinear coordinates.—| The element of 
+t D. E. Rutherford, Fector Methods (Oliver and Boyd, 1939), 


pp. 59-63. 
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length d/ is then given by 
dl? = hi dq +h? dq? +h? dq (1.3) 


. δεν γῶν" [ὅτι 
ἐπι )} (Ὁ) ὦ» 
and it can sk! be shown that 
ὦ. Ga ov) Ἢ 3 (ΠΝ zt) 
~ hy : hs (@q,\ hy Gq;/ 0q2\ hz 86, 


ὃ (hyh, sa 
4 ES 
(Fa h; 043 δὼ 


One method of solving Laplace’s equation consists of 
finding solutions of the type 


w= 0,(91)02(92)03(43) 
by substituting from (1.5) into (1.1). We then find that 


1 @ 2 (ints 22s) 4 2 #9 J (te 202) 
οι eq; hy 04; Q,0q2\ hy, 642 
yedrd 2 (he ih, 20») 
| J 9: 0q3\ hs 04, 
If, further, it so happens that 
Ἢ" =Si (4 Fi, 43) 


etc., then this last equation reduces to the form 


δία... 43) > | nan aut 4 


where 


QO, dq, 


+ F(q35 41) ok ik -{ fila) a 
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Now, in certain circumstances, it is possible to find three 
functions g;(q1), 92(92), 93(¢3) with the property that 


Fi(425 93)91(91) + F243, 91)92(92) + F3(41, 92)93(43) = 0. 


When this is so, it follows immediately that the solution 
of Laplace’s equation (1.1) reduces to the solution of three 
self-adjoint ordinary linear differential equations 


ra =0, (i=1, 2,3 1.6 
11} dq; —g,Q; (i ). ( * ) 


It is the study of differential equations of this kind which 
leads to the special functions of mathematical physics. 
The adjective “ special” is used in this connection because 
here we are not, as in analysis, concerned with the general 
properties of functions, but only with the properties of 
functions which arise in the solution of special problems. 

To take a particular case, consider the cylindrical polar 
coordinates (0, @, z) defined by the equations 

x=gcos ᾧ, y=esind, z=z 


for which A, = 1, hz = 0, h3 = 1. From equation (1.5) 
we see that, for these coordinates, Laplace’s equation is 
of the form 


Oe 1 6p Hee 1 dy εἰς δὺς 
do ode a? ad? 
If we now make the substitution 
Ψ = R(e)®(4)Z(z), (1.8) 
we find that equation (1.7) may be written in the form 
‘(4 it) Ls Lee 
do» odo) odd? Z dz? 
This shows that if ®, Z, R satisfy the equations 
εἰ ῷ 
dq? 


= 0. (1.7) 


+n?o = 0, (1.9a) 
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Ze 
ΓΞ —m?Z =0, (1.9b) 
d*R, 1dR 
teat (ma S)R=0 (190) 


respectively, then the function (1.8) is a solution of Laplace’s 
equation (1.7). The study of these ordinary differential 
equations will lead us to the special functions appropriate 
to this coordinate system. For instance, equation (1.9a) 
may be taken as the equation defining the circular functions. 
In this context sin(m@) is defined as that solution of 
(1.9a) which has value 0 when ¢ = 0 and cos (n@) as that 
which has value 1 when ¢ = Ὁ and the properties of the 
functions derived therefrom, cf. ex. 4 below. Similarly, 
equation (1.9b) defines the exponential functions. In actual 
practice we do not proceed in this way merely because 
we have already encountered these functions in another 
context and from their familiar properties studied their 
relation to equations (1.9a) and (1.9b). The situation with 
respect to equation (1.9c) is different; we cannot express 
its solution in terms of the elementary functions of analysis, 
as we were able to do with the other two equations. In 
this case we define new functions in terms of the solutions 
of this equation and by investigating the series solutions 
of the equations derive the properties of the functions so 
defined. Equation (1.9c) is called Bessel’s equation and 
solutions of it are called Bessel functions. Bessel functions 
are of great importance in theoretical physics; they are 
discussed in Chapter IV below. 


§ 2. Ordinary points of a linear differential equation. 
We shall have occasion to discuss ordinary linear differen- 
tial equations of the second order with variable coefficients 
whose solutions cannot be obtained in terms of the elemen- 
tary functions of mathematical analysis. In such cases 
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one of the standard procedures is to derive a pair of linearly 
independent solutions in the form of infinite series and from 
these series to compute tables of standard solutions. With 
the aid of such tables the solution appropriate to any 
given initial conditions may then be readily found. The 
object of this note is to outline briefly the procedure to 
be followed in these instances; for proofs of the theorems 
quoted the reader is referred to the standard textbooks.+ 

A function is called analytic at a point if it is possible 
to expand it in a Taylor series valid in some neighbourhood 
of the point. This is equivalent to saying that the function 
is single-valued and possesses derivatives of all orders at 
the point in question. In the equations we shall consider 
the coefficients will be analytic functions of the independent 
variable except possibly at certain isolated points. 

An ordinary point x = a of the second order differential 


equation 
y"+a(x)y'+Bx)y = 0 (2.1) 
is one at which the coefficients «, # are analytic functions. 
It can be shown that at any ordinary point every solution of 
the equation is analytic. Furthermore if the Taylor expan- 
sions of «(x) and B(x) are valid in the range |x—a [<R 
the Taylor expansion of the solution is valid for the same 
range. As a consequence, if a(x) and f(x) are polynomials 
in x the series solution of (2.1) is valid for ail values of x. 
When, as is usually the case, a(x) and f(x) are poly- 
nomials of low degree, the solution is most easily found 
by assuming a power series of the form 


i] 


y= 2, c(x—a)’ (2.2) 


r= 
for the solution and determining the coefficients cy, c,, Co, ..., 


+ See, for example, E. L. Ince, Ordinary Differential Equations, 
(Longmans, 1927), Chap. VIL; E. Goursat, A Course in Mathematical 
Analysis, Vol. Π, Part ἢ (Ginn, 1904), Chap. TT; J. C. Burkill, 
ΠΡ of Ordinary Differential Equations Chap. IV (Oliver and Boyd: 
1956). 
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by direct substitution of (2.2) into (2.1) and equating 
coefficients of successive powers of x to zero. 
The simplest equation of this type is 


y’+y =0. (2.3) 


Substituting a solution of the type (2.2) with a = 0 into this 
equation we find that, if the equation is to be satisfied, 


2, r(r—1)c,x"~? + Σ. cx" τα Q, 


The series on the left is equivalent to 


Σ, (r+1)(r+2)c,, 5x" 


so that, equating coefficients of x", we see that the equation 
is satisfied by a solution of type (2.2) provided that the 
coefficients are connected by the relation 


(r+1)(r+2)c,..+0¢, = 0. (2.4) 
The coefficients cp, ὧι are determined by the prescribed 
values of y, y’ at x = 0, and the others are determined 


by equation (2.4). From this relation it follows that the 
solution is 


Fe | ee 
Feet ss 4.2 re ee i 
y οἱ Tha )+ea(s 31° 5] ) (2.5) 

An equation of the kind (2.4) which determines the 


subsequent coefficients in terms of the first two is called a 
recurrence relation. 


§ 3. Regular singular points. If either of the functions 
a(x), B(x) is not analytic at the point x = a, we say 
that this point is a singular point of the differential equation. 
When the functions a(x), f(x) are of such a nature that 
the differential equation may be written in the form 


(x—a)*y"+(x—a)p(x)y'+q(x)y = 0 (3.1) 
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where p(x) and g(x) are analytic at the point x = @ we 
say that this point is a regular singular point of the differ- 
ential equation. 

If x = ais a regular singular point of the equation (3.1) 
it can be shown that there exists at least one solution of 
the form 


vi Σ ehxe=ay** (3.2) 


which is valid in some neighbourhood of x = a. More 
specifically, if the Taylor expansions for p(x), g(x) are 
valid for |x—a|<R, the solution (3.2) is valid in the 
same range. | 

Putting 


δ © 
P= 2 pdx—a)', a= Σ, adx—aY (3.3) 
r= = 
and substituting the expansions (3.2) and (3.3) into (3.1) 
we see that for the equation (3.1) to be satisfied we must 
have 


Σ οὐρτηζεν-ηρ αν" 


+ Ρίχ-- αὐ" Σ, ele x—ay 


Equating to zero the coefficient of (x—a)e we have the 
relation 


CoQ(Q—1)+ Polo P+ Golo = 9 
so that if co # 0 we have the quadratic equation 


0° +(Po—De+do = 0 (3.5) 


for the determination of ρ. This is known as the indicial 
equation. Similarly if we equate to zero the coefficient of 
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(x—a)**" we obtain the relation | }ι(χ) +(x) where 
Ρ oo 
cf{ot+tretr—1)+ 2. {p{o+r—s)+q,}¢,-, = 0 | yi) =(x—-a)" Σ᾽ o(x—ay’, 
is | tis (3.8) 
which may be written in the form yo(x) = y4(x).log(x—a)-+(x—a)* 2: (2 (x—a)! 
. ᾿ e001 


c,{(g+r\(et+r—1)+ pole t+r)+4o} 
᾿ Case (iii) 02 = @,—n where n is a positive integer. 
+ >} {p{ot+r—s)+q,}c,--, = 0. (3.6) In this case all the coefficients in one of the solutions 
s=1 | from some point onwards are either infinite or indeter- 
Equation (3.5) gives the two possible values @;, @2 of @. minate. It can be shown that the appropriate solutions are 
If we take one of these values, ὁ, say, and substitute it . ” 
in the recurrence relation (3.6) we obtain the corresponding y,(x) = (x—a)" 2 c{x—a)’, 
value of the coefficients c, and hence the solution | (3.9) 


© ya(x) = Jny1(x) log(x— a) Ἐ(χ -- αὐ Σ ",(χ--αὐ, 
yi(x) a > ο(χ-- αὐγαὶ, =0 
‘=e where 4, is the coefficient of x” in the cece of 


In a similar way the root g, of the indicial equation leads ΧΡ ἔπε ΚῈ - 4 y(u)d 
to the solution / fy,(x+a)}? P + 
y2(x) = 2, onx— 475. It may happen that g, = 0 in ae, case y2(x) does not 


contain a logarithmic term. 

Three distinct cases arise according to the nature of the 
roots of the indicial equation. § 4. The point at infinity. In many problems we wish 
to find solutions of differential equations of the type (2.1) 
Case (i) 0; —@2 neither zero nor an integer which are valid for large values of x. We seek solutions in 
In these circumstances the solutions y,(x) and y.(x) 
are linearly independent and the general solution of 


equation (3.1) is of the form 


the form of infinite series with variable τ If we make 
the transformation 


ae 
ξ 
the “ point at infinity” is taken into the origin on the 
é-axis. With this change of variable pagers (2.1) becomes 


So G-do( gees ὦ 
eu ge \e i+ pee 0 


y= 1, CRA) + ¥. el(x—ayt@, (3.7 


r 


is 


Case (ii) Q; = 02. 

If @; = @2 the solutions y,(x) and y,(x) are identical 
(except, possibly, for a multiplicative constant), The 
general solution of the equation can be shown to be 


a πὶ ...ὄ δου, ῖ».... τῷ5 ὁ τῶν 
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If 2€-*—&-7a(€-*), E~*B(E~*) are both 0(1) as €>0 and 
analytic in € then € = 0 is an ordinary point of equation 
(4.1) and we say that x = o is an ordinary point of 
equation (2.1). Returning to the original independent 
variable we see that the conditions ἴοι the point at infinity 
to be an ordinary point of equation (2.1) are that 


a(x) = : +0O(x~*), B(x) = O(x~*)as χοῦ. (4.2) 
The corresponding solutions are of the form 
7 po ‘a 
Similarly if, as x co 
a(x) = 9. +0(x"%), px) = FS ὁ ἐστ 3) 


where %, βὺ are constants we say Re the point at infinity 
is a regular singular point of the equation (2.1). The 
corresponding indicial equation is 


0° +(1—a)o+p =0. 


If the roots of this equation are @,, 02 the solutions of (2.1) 
valid for large values of x are of the form 


y= Yee, ylx)= Sere er, (4.4) 


§5. The gamma function and related functions. In 
developing series solutions of differential equations and 
in other formal calculations it is often convenient to make 
use of properties of gamma and beta functions. The 
integral 


r(n) = | "tS Ye (5.1) 
0 
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converges if n>0O and defines the gamma function. Similarly 
if m>0, n>0 the beta function is defined by the equation 


B(m, n) = | * "11 x)" Nd. (5.2) 


0 
It is then easily shown that ἢ 
(i) Td) = 1, 
Gi) [(@@+1) = aT), 
(iii) P'(n+1) = n! ifn is a positive integer, 
ἀπ 
(iv) B(m, n) =2 | sin?™~!@ cos*""' @ dd, 
0 


| oy Pen) ™) 
(v) B(m, n) = Tim-+n)’ 
(vi) TG) = Jz, 
(vii) [(p)T(i—p) = ncosec(px), O<p<1, 
(viii) Γ Γ(2η) = 2?7"-'T(n)T(n+4)—the duplication 
formula, 


nin? 
(ix) [(z+1)= Low πα Ὁ {πὲ (z>0). 


When ἡ is a negative fraction I'(m) is defined by means 
of equation (ii); for example 
Γ Γ 4Τ' 
πον τους Ὁ τῷ 
-}  (-3.-Ὁ }) 3 
By means of the result (ix) we can derive an interesting 


expression for Euler’s constant, γ, which is defined by the 
equation 


γ = lim (1+ lt Soli n) = 0.5772... (5.3) 
πεν 2 n 


+ For proofs of these results the reader is referred to R. P. Gillespie, 
Integration (Oliver and Boyd), 1951, pp. 90-95. 
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From (ix) we have 


d . ᾿ 1 1 1 
— tlog I(z+1)} = lim | logn— -Ξ..-.ὄ.--.-ὄ -.,, — 
az G+} tim ( » z+1 z+2 = 
so that letting z+0 we obtain the result 
d 
= —| — logI(z+1) ‘ 
ΙΖ og ΓίΖ Ὁ | (5.4) 


and from (5.1) we find 
ΡΝ Ἂ e~' log tdt (5.5) 
/ ἢ 
Integrating by parts we see that 
} e~' log tdt = +e77 log z+ Ε ἐς at 
z : 


so that 


z-~0 


—y =lim == — dt-+log 2) (5.6) 


Closely related to the gamma function are the 
exponential-integral ei(x) defined by the equation 


ei(x) = [ -- du (x>0), (5.7) 
and the logarithmic-integral li(x) defined by 
; * du 
li(x) = | 
(x) Ι τ (5.8) 


which are themselves connected by the relation 
ei(x) = —li(e™*), (5.9) 
Other integrals of importance are the sine and cosine 
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integrals Ci(x), Si(x); which are defined by the equations 


Ci(x) = -- = “du, Si(x)= (= δα, (640) 


and whose variation with x is shown in Fig. 1. 
In heat conduction problems solutions can often be 
expressed in terms of the error-function 


erf(x) = |. e~ du (5.11) 


Fig. 1. Variation of Ci(x) and Si(x) with x. 


whose variation with x is exhibited graphically in Fig, 2.7 
Similarly in problems of wave motion the Fresnel 
integrals 


C{x) = | cos (4nu?)du, S(x)= | sin (4nu7)du (5.12) 
0 0 
occur. The variation of these functions with x is shown 
in Fig. 3. 


+ A.C. Aitken, Statistical Mathematics (Oliver and Boyd, Seventh 
Edition, 1952), p. 62, gives a short table of values of erf (x). 


1-0 


> 
ΕΝ 


— erf(x) > 


φῳ 
ro 


0 μὰ 2.0 3-0 
—_— C—e 


Fig. 2, Variation of erf (x) with x. 


0 1 2 3 4 δ 


—= —s 
Fig. 3. Variation of the Fresnel integrals, C(x) and S(x), with x. 
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The importance of these functions lies in the fact that 
it is often possible to express solutions of physical problems 
in terms of them. The corresponding numerical values can 
then be obtained from works such as E. Jahnke and 
F. Emde, “* Funktionentafeln*’ (Teubner, Leipzig, 1933) in 
which they are tabulated. 

Examples | 


(1) Show that, in spherical polar coordinates r, 0, ᾧ 
defined by 


x = rsin@cos ᾧ, y = rsin θ sin ¢, z = rcos 8, 


Laplace’s equation becomes 


8 (2 dw 1 Μὴ Le 

-- --[ι᾿ὲ — sin θ) + 0, 

ar (, *) sin 0 G5 (* 80) sin? ap? 
and prove that it possesses solutions of the form r"e'"*© 
(cos 0), where O(,1) satisfies the ordinary differential equation 


2 

ἀμ): -μϑ εἰ be = 0. 
dy? -μ 

(2) Show that if 


dy 
x =acosh €cosy, y=asinhésiny, z=2z 


Laplace’s equation assumes the form 
2 ᾿ 2 
uty + nt +a*(cosh? €—cos? n) = = 0. 


Deduce that it has solutions of the form /(ié)f(n)e"* 
where /() satisfies the equation 


a +(G+16q cos 2η)[Ξ Ὁ 


in which G is a constant of separation and g = —a’y?/32. 


| 
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(3) Parabolic coordinates €, ἡ, @ are defined by 
x = /(&y) cos φ, y = ./(En) sin φ, z = (E—n). 


Show that in these coordinates Laplace’s equation becomes 


4 ὃ Μὴ ἐν 2( “4 1 ay 
aca aé)* En on\" Τ᾿ ag? 


Prove that if F,(x) is a solution of the equation 


_@P dF 
cP —_— -- — 
Bt ides (n 4x 7) P= : 


then F,(¢)F_,(y)e*'"® is a solution of Laplace’s equation. 
(4) Defining cos x, sin x to be the solutions of 
d*y 
= 0 
dx? i 
which respectively are 1,0 when x = 0, prove 


(i) cos(—x) = cos x, sin(—x) = —sin x; 
(1) cos (x+.x') = cos x cos x’—sin x sin x’; 
(iii) sin (x+x') = sin x cos x’+cos x sin x’: 
(iv) cos* x+sin? x = 1; 
ae ee 
(v) —(sin x) = cos x, < (cos x) = —sin x. 
A 


dx 
(5) The only singularities of the differential equation 


y"+ p(x)y’ +q(x)y = 0 


are regular singularities at x = 1 of exponents a, α' and 
at x = —1 of exponents β, β' the ΡΝ, αἱ infinity being 
an ordinary point. Prove that B = —a, β' = --α' and 
that the differential equation is 


(x? — 1)? yp” +2(x—1)(x—a—«a’')y’ +4an'y =0. 
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Show that the solution is 


y=ec =) +c ay 
"\ x41 *\x+1 


where c, and c, are constants. 
(6) Apply the method of solution in series to the 
equation 
“2 dy 
ἜΣ Awl 


showing that, near x = 0, y = Au+Bv where wu is a 
Maclaurin series and νυ = x'~“e*. (a is not an integer.) 


τα) > —y=0 


(7) Find two solutions of the equation 


2 | 
(x?4+2x) 2% 4 Y —ek-+ ty =0 
dx dx 


in the form 


i] 
= 2, ans, 
Show that, if k is a positive integer, one of these solutions 
is the polynomial 
ΚΕῚ (k+n—1)! (2x)" 
1+k(k+1 a 
( ) ὦ, (k—n+1)! (2n)! 


(8) Prove that if s is an integer and a is fractional 


fa-)=(-f τ 9. ey 


(9) Show that 
(ὃ («@—1)(@),-1 = (α--1}.; 


? ee eee 
(ii) (@),-, = (—1) ἄπ 
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. 


/ χὰ ἡ i 
(iii) ey 


=(-1)(-m),: 


where 
(x), = a(a+1)(a+2)...(x+r—1). 


(10) Prove that CHAPTER II 


2 a 
ei(x) = —y—log x+x— ὟΝ + ΧΩ ae HYPERGEOMETRIC FUNCTIONS 
and deduce that § 6. The hypergeometric series. The series 
x? ok ) (Pam 3 
Ci(x) aaa 7+log i 5.2] + 441 waar ys 1+ αν β κα α(α Ὁ ἡβί(β: 1) x? (6.1) 
2! 4! 1.7 1. 2y(y+1) 
Si(x) = χ -- x ἣν δ". ως is of great importance in mathematics. Since it is an obvious 
ls Del 


generalisation of the geometric series 


1+x4-x*+..., 


it is called the hypergeometric series. It is readily shown 
᾿ | that, provided y is not zero or a negative integer, the 
series is absolutely convergent if |x|>1, divergent if 
|x |>1, while if | x | = 1 the series converges absolutely 
if y>a+f.t It is convergent when x = —1, provided 
that y>oa+f—1. 
If we introduce the notation 


a ee ~1) = τὴ 
(a), = a(a+1)...(a+r—1) r@ (6.2) 
| we may write the series (6.1) in the form 
2Fila, B35 x)= Ye x", (6.3) 


the suffixes 2 and 1 denoting that there are two parameters 
of the type « and one of the type y. We shall generalise 


+ See J. M. Hyslop, Infinite Series, Fifth Edition (Oliver and Boyd, 
1954), p. 50. 
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this concept at a later stage (§ 12 below) but it is advisable 
at this stage to denote the “ ordinary” hypergeometric 
function by the symbol 2} instead of simply F, if we are 
to avoid confusion later. From the definition (6.3) it is 
obvious that 


2F,(B, αν y; x) = 2F (α, B; y; x). (6.4) 


A significant property of the hypergeometric series 
follows pearorsey from the definition (6.3). We have 


2 δία, B37; x)= Σ OMB yes 
Σ 1 (r—1)!(), 
— τ @r+iB)+1 «χ᾽ 
r=0 γ](γ)),.: 
Now (α),... = «(a+1), so the right-hand side of the last 
equation becomes 
ap x (a+1),B+1), Ὁ 
ἡ ταῦ r\(y +1), 
showing that 


d . 
= 2F y(t B73 3) = 2Fi(a+1, B+1; +1; x). (6.5) 


It should also be observed that 
2F (α, B; y; 0) Ξ 1 (6.6) 
so that 
d ap 
— oF (a, B; ¥; Ξε --- : 
4. (a, Bs 7 εἰ - ' (6.1) 


Several well-known elementary functions can be 
expressed as hypergeometric series; examples of them are 
given in ex. 1 below. 

It should be noted that, if we adopt a certain convention, 
a hypergeometric series can stop and start again after a 
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number of zero terms. For example, consider the hyper- 
geometric series ,F,(—n; δ᾽ —n—m; x) where both m 
and m are positive integers and ὃ is neither zero nor a 
negative integer. Because of the occurrence of (—v), in 
the numerator in the expansion in powers of x it is obvious 
that the (n+1)th term of the expansion will be zero, and 
we are tempted to think that every subsequent term is 
also zero. If we note that, as a result of ex. 9 (ili) of 
Chapter I, 


(—n), a τὰς 
ΩΣ en (n+m—r)(n+m—r—1)...(n—r+1) 
(6.8) 


when the form on the left is not of type 0/0, and if, further, 
we assume that it still has the value on the left when it 
is indeterminate, we see that we may write 


2F,(—n, δὲ —n—m; x) 


CS ee ae ey Pe (b),x’ 
DAC ΕΝ πεπρί ΘῈ r! ala 


so that although the series stops at the mth term it starts 
up again at the (2+m+1)th term. For instance, 


oF; (—2,.1; —5; x) 


2 1 τι ee | 
=1 χὰ —x?— —x®°—=x'—x* +... 
τς 5 * 10 10 5 
According to a different convention, however, the 
hypergeometric function does not restart after a set of 
zero terms. 


§ 7. An integral formula for the hypergeometric series. 
In order to derive some further properties of the hyper- 
geometric series we shall first of all establish an expression 
for the series in the form of an integral. It is readily shown 
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that 

(6), _ BB+r, y—B) _ 1 

2), 8Β(β,γ-βὺὶ Βιβ,γ--β) 
from which it follows that 

2F (a, B; γ; x) 


| (1 — tyr F- 1btr-1 dt 
0 


geses © ΒΥ Ὡς ay-t~1ypte- 
BB veo ΤΡ |, OO 


Interchanging the order in which the operations of sum- 
mation and integration are performed we see that 


2F (a, B; γ; x) 


1 1 | 2s Ὁ < (a) | 
= -ς--.-.--.---ς- 1—1t)?~8-148-! ——(xt)"> dt. 
B(B, y—B) x ΙΣ r! crf μὰ 
Using the fact that 


Σ ω, (xt)’ = (1—x1)™, 
r=o r! 


we have the integral formula 
2F\(a, B; y; x) 


1 ᾿ | 
gee eee | _ y-B-1,8-1 a —@ | ) 
B(B, »— B) | 0 oils Nee rd 


valid if | x|<1, y>f>0. The results hold if x is complex 
provided that we choose the branch of (1—x?)~* in such a 
way that (1—xf)~*-1 as 0 and &(y)>A(p)>0. 

The first application of (7.1) is the derivation of the value 
of the hypergeometric series with unit argument. Putting 
x = 1 in (7.1) we have 

1 


2 F(a, β: " 1) ΞΞ BB. γ-β) |. (L— gyre Fi et dt 
it Π 0 


ms B(B, γ--ἅ-- ἢ 
Β(β, γ-- βῸ 
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if y-«a—B>0, B>0. Ifwe express the beta function in terms 
of gamma functions we have Gauss’s Theorem 


ΓΩΓίγ) --α-- β) : 
Pit $33 De Στ -τι (7.2) 
3 T(y—a)T(y—B) 
Now if « = —n, a negative integer, we have 


so that equation (7.2) reduces to 


2F(—n, B37; D= O-Bn 
(7), 
which is known, in elementary mathematics, as Vander- 
monde’s theorem. 
Again, if we put x = —1 and « = 1+/f-—y we have, 
from equation (7.1) 


Βα, B; B—a+1; -1) = ἘΠ ΈΡΞ ΟΣ | aay tae 
ΞΞ 0 


If we write € = 12 in this integral we see that its value is 
4B(4B, 1—«). Using this result and the relation 41(4f)/T(f) 
= T(1+4)/l'(1+) we have Kummer’s theorem 


Γ(Ἐβ-- ΓΕ Ἐ2β) (7.3) 
ΓΑΙ ἘΡΊΓΑ - 3 --αὶ : 


Further, we can deduce from the formula (7.1) relations 
between hypergeometric series of argument x and those 
of argument x/(x—1). Putting τ = 1—7 in equation (7.1), 
and noting that 


(1—x(1~9)"*= 0-3" f1- —— x 
x-—1 


2F (a, B, B—a+1; —1)= 
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we see that 
2F (a, B; y; x) 


ap), |, d-orteret fi oh "ας 
δ | 


~ BOB, y—B) χ--ὶ 
.« τῶ κι. τ δι ae, 
Fy Β-θ, Dak (αν 7-85 =), 


whence we have the relation 


2F,(a, β; 73 x) = (1—x)7*2F, (« γ--β; γ; 2) (7.4) 


and, by symmetry, the relation 
2 (α, B; y; x) =(1—x)~* 2 F, (-α β:γ: 4) (7.5) 


Using the symmetry relation (6.4) and equation (7.4) with 
x replaced by x/(x—1) we see that 


» ἡ νει ον ἐκ a 2 
2 ( y-B; 7 4) 2 ι βια;γ; =) 
=(1-x)’"* ,F,(y—B, γ-- αὶ γ; x), 
so that 
2F (a, B; γ; x) = (1—x)’"*"* ,F(y—a, y—B; 7; x). (7.6) 
If we put x = 4 in equation (7.4) we obtain the relation 
ΔΕ ία, B; 7; 3) = 2°. F(a, y—B; y; — 1). 


The series on the right-hand side of this equation can be 
derived from equation (7.3) provided either that 


y= y—B—a+1, ie. B = 1—«, 
or that 
y= a—(y—f)+1, i.e. y = 4(a+8+1). 
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We then obtain the formulae 
| ray)Gr+4) 
F(a, 1-—«; γ; 4) = ———  --- 7. 
Mie I-Gnd= nate βδνκῶνι 


| aay P@UG +4448) 


§ 8. The hypergeometric equation. In certain problems 
it is possible to reduce the solution to that of solving the 
second order linear differential equation 


2 
x(1—x) oo +{y—(1+a+)x} dy —aBy=0 (8.1) 
dx dx 
in which ἃ, 6 and y are constants. For instance, the 
Schrédinger equation for a symmetrical-top molecule, 
which is of importance in the theory of molecular spectra, 
can, by simple transformations, be reduced to this type. An 
equation of this type also arises in the study of the flow of 
compressible fluids. In addition certain other differential 
equations (such as that occurring in ex. 1 of Chapter ἢ 
which arise in the solution of boundary value problems in 
mathematical physics can, by a simple change of variable, 
be transformed to an equation of type (8.1). Indeed it can 
be shown that any ordinary linear differential equation 
of the second order whose only singular points are regular 
singular points, one of which may be the point at infinity, 
can be transformed to the form (8.1). For that reason it 
is desirable to investigate the nature of the solutions of 
this equation, which is called the hypergeometric equation. 
We may write the hypergeometric equation in the form 
x?y"+x(1+x+x7+...){p—(a+B+1)x}y’ 
—aPx(1+x+x7+...)y =0, 
+ See, for example, L. Pauling and E. B. Wilson, Introduction to 
uantum Mechanics, with Applications to Chemistry (McGraw-Hill, 


OQ 
New York, 1935), pp. 275-280, and ex. 10 below. 
σ 
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so that, in the notation of § 3, we see that near x = 0 


Po=?, Go=9, 
and the indicial equation is 
oe’ +(y—1)e =0 
with roots g@ = 0 and g = 1 --. 
Similarly, the equation can be put in the form 
(x ~1)?y”—(x—1) {py -—a—B—1—y(x-1)+...}y’ 
+aP(x—1){1—(x—1)+...}y = 0, 
with indicial equation 
eo’ +(a+f—y)e = 0, 
of which the roots are 9 = 0, 9 = y—a—f. 
Finally in the notation of § 4 we have for large values 
of x 


a(x)~ StP+D gob 
x x 
and so the indicial equation appropriate to the point at 
infinity is 
0’ —(a+ B)o+ap = 0, 

with roots a, β. 

Thus the regular singular points of the hypergeometric 
equation are:— 

(i) x = 0 with exponents 0, 1—y; 

(ii) x = oo with exponents ἃ, f; 

(iii) x = 1 with exponents 0, y—a—f. 
These facts are exhibited symbolically by denoting the 
most general solution of the hypergeometric equation by 
a scheme of the form 


0 oO 1 
y=P) 0 α 0 "ee (8.2) 
at hele i es 
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The symbol on the right is called the Riemann-P-function 
of the equation. 

We shall now consider the form of the solutions in the 
neighbourhood of the regular singular points. 

(a) x = 0: Corresponding to the root g@ = 0 we have 
a solution of the form 


yo > Ox. 
r=0 
Substituting this series into equation (8.1) we obtain the 
relation 
(i—x) ¥ ¢,r(r—1)x"™? 
r=0 

+{y—(a+B+1)x} 2, c,rx’~1—aB 2 cx’ =0 

which is readily seen to be equivalent to 


Σ {erasfr(r+1)+(rt1y]—o(r-+a)(r-+B)}x" = 0, 
r=0 


so that 
Cr+ i= (r-+a)(r+B) δι, (8.3) 
(r+1)(r+7) 
from which it follows that 
ῦ, = (2),(B), Co: (8.4) 
(γ),.} 

It follows that the solution which reduces to unity when 

x = Oils 

θ᾽... Met OPE ΤῸ 025. 

y1! (7+ 1)2! 

i.€. y = 2F,(@, B; γ; x). (8.5) 
Similarly, if --- is not zero nor a positive nor negative 


y=1+ 


πὶ τ τας 
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integer, the solution corresponding to the root 9g = 1—y is 


y= % Re ei Ba 
r=0 
where 


a) 
(1—x) Σ, σ(Ρ Ἐ1.--- }7ὁὦ.7ὕν -- 7) χη 7. 
οὐ ΠῚ 
+{9-@+B+)x} Σ o(r+i—yx"?—ap Σ οχ'ττλ = 0, 
r= r=0 
which is equivalent to 


᾿ o{(r-+1—y)\(r—y) +9(r+1-y)}x°77- 


r= 


- = ¢{(r+1—7)(r—y) + (e+ B+ ΧΡ Έ1.-- 7} αβ)χη Ὑ}} =0, 


implying that 
(r+a—y+1(r+h—y+1) Ξ 
(r+1)(r+2—y) J 
Comparing this relation with (8.3) and taking ῳ = 1 
we see that this solution is 
x*~?, F,(a—y+1, B—y+1; 2-7; x). (8.6) 
Combining equations (8.5) and (8.6) we see that the 
general solution valid in the neighbourhood of the origin is 
y=A2F,(a, B; y; x) 
+Bx'~?,F,(a—y+1, B—y+1; 2-7; x), (8.7) 
provided that 1—y is not zero or a positive integer. 


If y = 1, the solutions (8.5) and (8.6) are identical. If 
we write 


Cr+i ᾿Ξ. 


}ι(Χ) = 2 (a, B; γ; x) 
and put 


Y2(x) = γι) log x+ > Gx" 
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we find on substituting in (8.1), with y = 1, that 
2 (@)(B) (aB—a—B—r) _ 
σε θο, rat B+ Net SET 0 


from which the coefficients c, may be determined. 
A similar procedure holds when 1 —} is a positive integer. 
(b) x = 1: If we let € = 1— x, equation (8.1) reduces to 


ξ(! -- -9 52 +{a+p- —7+1—(a+ p+ Iie} 2 —apy = =0 


which is identical with equation (8.1) with » replaced by 
a+B—y+1, and x by €=1—x. Hence it follows from 
equation (8.7) that the required solution is 
y=A2F,(a, B; a+B—y+1; 1—x) 
+ B(i—x)'"*~* ,F,(y—a, y—B; y—a—B+1; 1—x). (8.8) 
(c) x = οὐ: Corresponding to the root o = a, we put 
y= 2. oo ede 


which gives 
(l—x) Σ᾽ ¢(r+a)(r+a4+))x-"*"! 
r=0 


—{y-(@+P+1)x} Yo (rtadex ttf YY δ,» "πε == 0, 


1.€., 

oo oo 

Σ ο(ν τα) α--γῈ 0) χ 5 t= Σ, corr α-- β)χ ᾿ 
whence it follows that 
(r+a)r+o-y+1) 
(r+1)r+a—B+1) ” 


Cr+1 st 


; 
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which in turn is equivalent to 
δ. = (a) (a—y+ 1), δ 
———— ἢ γος 2 ῦ" 
r W(x -β Ἔ 1), 
Taking ᾧ = 1 we obtain the solution 


x *4F, (« a—y+1;«¢—B+1; +) 


From the symmetry we see that the other solution is 


xf (, B-y+1; p-a+;+), 
᾿ * 


so that the required solution is 


y= Ax *4F; ( a—y+1;a—B+1; 1) 
. (8.9) 


+Bx~’,F, ( ᾿ B-y+1; p—a+1;+) 
x 


§ 9. Linear relations between the solutions of the hyper- 
geometric equation. The series in the solution (8.7) are 
convergent if ah lbhe ie. in the interval (—1, 1) whereas 
those in the solution (8.8) are convergent in (0, 2). 
There is therefore an interval, namely (0, 1), in which all 
four series converge, and since only two solutions of the 
differential equation are linearly independent it follows 
that there must be a linear relation, valid if 0<x<1, 
between solutions of type (8.7) and those of type (8.8). 

Let 


2F (4, B; y; x) = A, F,(@, B; 2+ B—y+1; 1—x) 
+ B(1—x)’~*"* 5 Fy(y—a, y—B; y—-a—B +1; 1—»), 
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then putting x = 0 we have 
1=A>,F,(a, B; a+ B—y+1; 1) 
+B,F,(y—«, y—B; yp—a—B+1; 1), 
and putting x = 1 we have 
2F (a, B; y; 1) = A, 
if we assume that 
I>y>at+f. (9.1) 
Substituting for the series with unit argument from equa- 
tion (7.2) we see that 
4 — EO G-a—B) 
I(y—a)P(y — B) 


and that 
1=A Va+B-7+ DPU—y) +B Γὼγ --οα-- β ὁ ἩΓ( --Ὁ 
Γ(β--γῈΓα - γε)  6Ὸ ΓΑ ρΓᾺ --αὶ 
so that 


- ΓΟ Γ(αἘβ--) 
Γ(Γ(β) ᾿ 
whence we find that 
2F (a, B; }> x) 
ΤΩ ΓΟ--α-- βὺ pp. 
r-ar@—p) λον 
Γ Γ(α-ἐ β-- Ὁ) 4. .νν--ατῇ 
δ. - 
2F(y—a, y—B; yp—a—B+1; 1—x), (9.2) 
provided that the condition (9.1) is satisfied and 0<x<1. 


a+B—y+1; 1—x) 
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If we replace x by : in equation (9.2) we have 


2F; (« β:γ; 
Χ 


ἘΝῚ A 2 eat oF; (« B; a+ B—y+1; 1- ἢ 


Γ(γ--ὐΓ( --β) 
ὡ ΓΟ Γ(α Ἐβ--νῦ} [ - Ἵππον 
Γ(αΓ(β) Χ 


1 
2 Τὶ (- γ.-}.γ7--ἀ-ΒῈ1.1- +} 


and from equation (7.4) 
1 ᾿ 
2F (« B37; 1- 2) = x" 2 F(a, y—B; y; 1—x), 


so that 


F(a pi 7:4) 


(yl Q—«—-B) « : 

= F(a, 2— 7 +1; 0+ B—p +1; 1—x) 

Py—a)T—fp) * 
Γ(Γί(α -Ἐ β-- } x(x ae 1)“: 
Γ(ΟΓ(β) 
2Fi(y—a, 1—a; yp-a—-B+1;1-x), (9.3) 
where 1<x<2 and l>y>a+. 

These relations are typical of a larger number which 
exist between the solutions of the hypergeometric equation 
(8.1). If we change the independent variable in this 
equation to any one of 


§ 10 HYPERGEOMETRIC FUNCTIONS 33 


the equation transforms to one of the same type (but, of 
course, with different parameters). The equation (8.1) 
therefore has twelve solutions of the types (8.5) and (8.6) 
—two for each independent variable—each convergent 
within the unit circle. Any one of these can be expressed 
in terms of two fundamental solutions. In addition twelve 
more solutions of the kinds 

(1—x)""9"? 2 F,(y—a, γ-- βὶ 75 >); 

χ᾽ "--χ) τ 2Ρ.(1--α, 1-- β; 2-7; x) 
can be derived. The relations between these twenty-four 
solutions of the hypergeometric equation are of the types 
(9.2) and (9.3); for a full discussion of them the reader is 
referred to T. M. MacRobert, Functions of a Complex 
Variable (Macmillan, 2nd edition), pp. 298-301. 


ἢ 10. Relations of contiguity. Certain simple relations 
exist between hypergeometric functions whose parameters 
differ by +1. For example if the parameters « and f 
remain fixed and y is varied we can prove that 

γίγ--1--(2)--1--α-- β)χΧ} 2" Πα, B; γ; x) 

+(y—a)(y— β)χ 2Fi (a, B; y+1; x) 
-- ( -- 1)( --χ) 2 Ε Π(α, B; y-1; x) =0. (10.1) 
The proof follows from the definition (6.3). For the co- 
efficient of x” in the expansion of the function on the 
left of (10.1) is 
(ct),( 2), (α), -- (β), -ἴἰ 
1-2 ee Oa PD a 
(γ)),! ())»..-- τοι -- 1}} 
. (α), = 1(B), a-1 (x), (B),, 
+(y—a)(y—p) Sams Dams fy — 4) SSI 
(y+1),-1(2—-D)! (y—1)qn! 
(),-1(B),—1 
εγο τ ἢ enn — 
(y—1),-1(n—1)! 
and it is not difficult to show that this is zero. 
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In another kind 8 and y are kept constant and « is 
varied. One such is 


{y—a—B+(B—a)(1—x)} 2 (α, B; γ; x) 
+a(1—x)2F,(a+1, B; y; x)-(y—@) 2F(a—1, B; γ; x) =0 
(10.2) 
the proof of which is similarly direct. 
In the third type of relation y is kept constant and ἃ 
and f vary. One of the simplest among these relations is 


(α-- β)2 Γῷ(α, B; γ; x) =a2F, (+1, B; 7; x) 


—B2F (a, B+1;7;x) (10.3) 


The proof of these relations is left to the reader; further 
examples are given below (exs. 3, 4). 


§ 11. The confluent hypergeometric function. If we 


replace x by x/B in equation (8.1) we see that the hyper- 
geometric function 
2F,(a, B; 7; x/B) 


is a solution of the differential equation 


d*y | | a+1\ | dy 
Ἐπ, Χ es — —_ <= — ἘΞΕΣ 
. ( *) dx? a (1 Ἐ β ) xf dx oii 


so that letting Boo we see that the function 
lim 2 F(a, B; γ; x/B) (11.1) 
fa 


is a solution of the differential equation 


From the definition of (6),, we see that 


hs (ῥ), . 
al οσύς 
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so that the function (11.1) is the series 


= (a), x 
.:- [11.3 
r=o(y), Τ' sates 
and this series we denote by the symbol ,F,(a; y; x). This 
function is called a confluent hypergeometric function, and 
the equation (11.2) is the confluent hypergeometric equation. 

Equations of the type (11.2) occur in mathematical 
physics in the discussion of boundary value problems in 
potential theory, and in the theory of atomic collisions 
(see examples 13, 14 below). 

It is readily verified that the point x = 0 is a regular 
point of the differential equation (11.2) and that, in the 
notation of ὃ 3, pp = y and 40 = 0. The indicial equation 
is therefore 

o(e+y—1)=0 


with roots oe = 0 and 9 =  --. 
Corresponding to the root g = 0 there is a solution of 


the form 
oD 


yi ™= 2 Cx; 


,s 
substituting this solution in equation (11.2) and equating 
to zero the coefficient of x" we find that 

_ (e+r)je, 

ΕΣ 

Putting co = 1 we see that 


Cr+4 


and if y is neither zero nor a negative integer the solution is 


γι) = ,F,(@; γ; x). (11.4) 
Similarly, the root 0 = 1—7, leads, if 1—+y is neither zero 
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nor a positive integer, to a solution of the type 
‘oo 
vax? Σ ον 
r=0 
If we write 


y2(x) = x'~7u(x), 

and substitute in (11.2) we find that u(x) satisfies the 
equation 

d*u du 

πεν Ὅν, oh tie —(a—y+1)u = 0, 

which is the same as equation (11.2) with y replaced by 
2—y and ἃ replaced by «—7+1. We know from equation 
(11.4) that the solution of this equation which has value 
unity when x = Oisu = ,F,(a—y+1; 2—y; x) so that 

y2(x) = x'~?  F,(a—y+1; 2—7; x). (11.5) 
Thus if » is neither 0 nor an integer the general solution 
of equation (11.2) is 
y(x) = As Fy(@; 93 x) + Bx'~? Fy(a—y +1; 2—7; x), (11.6) 


where A and 8 are arbitrary constants. 
In the exceptional case y = 1 we have 


yi(x) = 1 F(a; 1; x), (11.7) 
obtained simply by putting y = 1 in equation (11.4). For 
the second solution we write 

y2(x) = y,(x) log x+ 2 Cx. (11.8) 
Substituting this expression in equation (11.2) we find that 
the unknown coefficients c, must be such that 
a —yyt+ Σ, { Ἐ1)26, ει το, χ' +c, =0. 


Inserting the value of y,(x) from equation (11.7) we see 
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that these coefficients are determined by the recurrence 
relation 


ὃ, = 1--αἰ (r+1)*c,4,—re, = (1—a) oe (11.9) 


The complete solution is therefore given by y = Ay,(x)+ 
By,(x) where A and B are arbitrary constants and the 
functions y,(x), y.(x) are defined by equations (11.7), 
(11.8) and (11.9). The complete solution when y is an 
integer may be found by a similar method. 

If in equation (11.2) we put 


y(x) = x7 37 e** W(x) (11.10) 


we find that the function W(x) satisfies the differential 
equation 


2 2 
ἘΝ \-#4 a τς Ως =0, (11.11) 
x > x 
where we have written k for 4y—a and m for (4—4)). 
The solutions of this equation are known as Whittaker’s 
confluent hypergeometric functions. 

If 2m is neither 1 nor an integer the solutions of the 
confluent hypergeometric equation corresponding to equa- 
tion (11.11) are given by equation (11.6) with » = 1+2m 
and «a = 4—k-+m. Thus the solutions of equation (11.11) 
are the Whittaker functions 


Mg, w(x) = x2 πε , F,(4—k-+m; 142m; x), (11.128) 
My, - κ(χ) = χε τ. F,(4—k—m; 1-2m; x). (11.120) 
Several of the properties of 2} functions have analogues 


for the ,F,; functions. Corresponding to equation (7.1) 
there is the integral formula 


gee seh en 1 ” (1 1rr8-1e- prt 
ΚΕ{(α; γ; χ) = rarer i) *-1y""Lextgr (11.13) 
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from which Kummer’s relation 
1F(a; γ; x) = & ΡΠ) -τ-αἰ γ; —x) (11.14) 


may be obtained by a simple change of variable. The 
analogue of equation (6.5) is 


d 
= Ὁ ι; 75 9)} = ; (Fiat+i;y+i;x), (11.15) 


while corresponding to the contiguity relations of § 10 we 
have relations of the type 
αι Fy(a+1; y+1; x)+(Q—2) 1 F,(@; γῈ1; x) 
—7 F(a; 73x) =0, (11.16) 
(x+a) Fy(at+1; y+; x)+(—o) Fy; y +1; x) 
—7,Fi@+1;7;x)=0, (11.17) 
%Fy(a+1; γγ; x)+(y—2a—x) , Πα; γ; x) 
+(4—y) Fi(a—-15 γ; x)=0, = (11.18) 


(a—y)x F(a; y+1; x) +9(x+7—1) 1 Fi(@; 9; x) 
+y(y—1), F(a; y—-1;x)=0. (11.19) 


§ 12. Generalised hypergeometric series. There are 
two ways by which we may approach the problem of 
generalising the idea of a hypergeometric function. We 
may think of such a function as being the solution of a 
linear differential equation which is an immediate general- 
isation of the equation (8.1) or we can define the function 
by a series which is analogous to the series (6.1). 

At first sight it is difficult to see how the differential 
equation (8.1) can be generalised immediately, but if we 
introduce the operator 
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and notice that (8.1) is equivalent to 
{9(9+y—1)—x(9+a)(9+f)}y Ξε 0, (12.1) 


an obvious generalisation is 


{9(9+0;—1)...(9+0,—1)—x(9-+0)...9+ap44)}y = 0, 
(12.2) 


Where 1, ἀν...» ἄρξει» Qi, ---» @, are constants. Further- 
more it is readily shown that this equation is satisfied by 
the series 

τ᾿ (αιλι(α))....(αρ: wn : x (12.3) 

n=0 (ρ,),(02),...(0.), 45! 

which is, itself, a generalisation of the series (6.1). Such 
a series is called a generalised hypergeometric series and 
is denoted by the symbol 4, F (1, ..-, %pi13 Q1s ...-» Op} X). 
It is left as an exercise to the reader to show that, if no 
two of the numbers 1, 0, 02» ..., 0, differ by an integer 
(or zero), the other p linearly independent solutions of 
equation (12.2) are 


x πο (1 +0, — Qj... 14+0,41-0:3 2-0, 
1+21—02, rey 1+0,—0;; x), (i = 1, 2, 448} n). 

As it stands (12.3) is a generalisation of the series (6.1) 
but it is not sufficiently wide to cover a simple series of 
the type (11.3). To cover such cases we generalise, not the 
differential equation, but the series defining the function. 
The generalisation of (6.1) which includes (12.3) is the series 

bd | ἢ 
(α.),(α)),. ᾿ (ἀρ), ς x . (12. 4) 
n=0 (@ t)n(O2)n--(Og)n n! 
which we denote by the symbol 
gE Ais +<Ggs Oss ---Oes Ἅ}, 


or, if we wish particularly to throw into relief the difference 
between the numerator and the denominator parameters, 
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αι [αν -.. αν; *| 

pi is sey Og; 

The suffix p in front of the F denotes that there are p 

numerator parameters ἃ), ..., ἄρ; similarly the suffix q 

indicates the number of denominator parameters. 
Generalised hypergeometric series do not usually arise 

in mathematical physics because we have to solve equations 

of the type (12.2). Their use is more indirect. Such series 

occur normally only in the evaluation of integrals involving 

special functions. In certain cases these series reduce to 


series of the type | 
Gis G5 
ele aes 
Oty vers Og; 


which have unit argument. For this reason it is desirable 
to have information about sums of this type. An account 
of the theory of such sums is given in W. N. Bailey, 
Generalised Hypergeometric Series (Cambridge University 
Press, 1935). Here we shall consider only one such calcula- 
tion because it illustrates the use of the theorems of Gauss 
and Kummer proved above (equations (7.2) and (7.3) 
respectively). Other results of this kind are given in 
examples 18 and 20 below. 
By expanding the ,F, series involved we see that 


Γ(Γ(ΡΊΓ() | i ἢ oy ‘| 
Γ(--α-- β)Γ((--α--γ)}ἢ : l+a—B, 1l+a-—7; | 
- 5 —PetolGtarotn) 

n=o niT(1l+a—6+n)0+a—y+n) 
= T(at+nl(B+n)0(y+n) 
n=o AID(1+a+2n)(i+a—B—y) 
T(1+«+2n)I'(1+a—B-y) 
[Mii+e—B+n)+a—-y+n) 


by the symbol 


5 
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Now by Gauss’s theorem (7.2) the expression inside the 
curly brackets is equal to 2 ΕΠ(β Ἐπ, p+n; 1+a%4+2n; 1) 
which may be written 
= [(p+n+ml(yt+n+m0(i+e+2n), 
m=o T(B+n)l(y+n)(i+a+2n+m)m! ἢ 
whence we find that 


αἱ δ S Vet+tn(p+n+mlyt+n+m) 
nso m= 0 T(1+a+2n+m)(1+a—B—y)n!im! 


Interchanging the order of summation and putting p = n-+m 
we see that 


5. 5 τὠτργοτὴ ᾧ. Γ(α τη) ἱ 
p=0 T(i+a—fB-—y) n=0 n\(p—n) IT (1+a+n-+ p) 


Now by example 9 (iii) of Chapter I 


(-1yp CP» 
car 1) i 


so that the inner sum is equal to 


Γ(α) a, —Pp; ry 
pir(it+a+p)° ἢ Os 
which by Kummer’s theorem (7.3) is equal to 
P@PA+se) 
piT(1+ae)P(1+4a+ p) 
Therefore 
᾿ ToT (B+ P+ PTA Ἐ3.). 
p=o p\(ita—B—y)T+o)l(1+40+4 p) 


_ __ T@r@ro) F Be γ; ‘| 
x οὐ τα βΞ)" *| 1440; 1 
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This ,F; series with unit argument can be summed by (2) By transforming the equation y"+n’y =0 to 
Gauss’s formula (7.2) and the expression for S found. hypergeometric form by the substitution ἔ = sin* z, prove 
It then follows that that, if -4inszS4z, 

a fp, 9 | (i) cos (nz) = 2F,(4n, —4n; 4; sin? 2); 
ὅ2|1.α--  , 1+0—9; β (ii) sin (nz) = nsin z 2: -- 3, $+40; 3; sin? 2); 


and that, if O<z<7, 
(iii) cos (nz) = cos ($num) 2F,(4n, —4n; 4; cos” z) 
+n sin (4nn) cos (2) 2F\(4—4n, $+4n; 3; cos? 2); 
(iv) sin (nz) = sin (4nn) 2F,(4n, —4n; 4; cos? z) 
—n cos ($n) cos (z) 2F,(4—4n, 4+4n; 3; cos? 2); 
(3) Prove the relations: 
(i) (α-- β)( -- x) 2Ε,(α, B; 7; x) 
= (y—f)2F (a, B—1; 7; x)—(y—%)2F(a—1, B; 7; x); 
(ii) (γ -- β -- 1). (α, B; 7; x) 
=(y—a—B—1).F,(@, B+1; 7; N+oeU—x)FiG@+1, βῈῈ1;γ; x); 
=(%—B—1)(1—x).F,(@,B +137; )+Q—-®2Fi(¢—-1,8+1; 9; x); 


(iii) (7-a—f)2F (4, B; γ; x) 


_ T+40)M(1+4a—f-y) (1 +a—f) i +a—y) 
Γ( ἘΟΓ( Ἐα-- β--)Γ( -Ε3α-- βΓ( Ἐ3α--ν} 


a result which is known as Dixon’s theorem. 
Examples II 
(1) Show that 
(i) 2 (α, B; B; 2) = (1—z)~*; 
(il) 2F, Go, 4a+4; 1; 2) = 53(({--2. Ἐ{Γ{Ἐ2) ἢ); 


(iii) 2}, 6.2, 3. Ἐ; 3; 2%) = εἰ {(1—z)"*-(1+2)"*}; 


(iv) 2F,(1, 1; 2; 2) = — : log (1-- 2); 


ΚΗ τ 2. ὦ 
ἡ μεν τω ary = (y—a)2Fy(a—1, B; γ; x)-BU—)2Fi@, B+13 95»); 
(iv) o.F,(a+1; B3 93 O—-G—D2Fi@, Bs y—-13 x) 


1, 
(vi) AiG, 2; 3; 2°) = a = (2+1-y).F,(«, Bs γ; χ); 
: ) (1—x) 2Fi(@, Bs γ; χ)-- 2Fi(@—1, B-1; 9; 
(vii) »Fy(d, 1; 3; -- δῦ) = 82, (v) (1-x) F(a, Bs 93 x)— 2Fi@—1, B-13 93») 
= at B-7—-! ΟΡ, Bs y+1; x); 
? 


(viii) 2}, 3; 1; 2) = : κρὸ; 
(vi) ΞΞΕΕ 2F (a Bs γε1; χ) 


Pa 2 
(x) 2F\(—4, 4s 5B) = = BO. = ,F,(a-1, B—1; 93 x)— 2F,(a, B-13 9; χ); 
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(vii) 1—x)2Fi(@%, Bsy;x) 
= 2F,(a, B—1; 7; x)+ i ΝΒ 2F (a, B; y+1; x). 
4. Prove that ; 
(ὃ 2Fi(, B+1; y+1; 2)-- ὃ" ,(α, B; γ; 2) 
. α(γ--β) 
r(y+1) 
(ii) 2Fi(@, B; γ; 2) = 2F,@+1, B—1; 7; 2) 
a—B+1 
? 


Z2F,(a+1, B+1; y+2; 2); 


+ ——- 2 ΒΑ ΤΊ, B; y+1; 2). 


Deduce a simple expression for the hypergeometric series 
γ᾿, τία, β; β- 1; 2). 


(5) If n is a positive integer, prove that 
2F\(—n, a+n; 7; x) 
x-"(1— xT) do tent 
ἑω ae AY fa κ« amiey a—-ytn . 
I'(y+n) dx" x Rie 
and deduce that 
oF (-- χη; 3. Ὁ}; te) 
_ - 1)*~ **°T(4a+4) a" — (4? — Sak oi . 
2"T(40+4+n) ἀμ" 
(6) If π is a positive integer, and | x |>1, prove that 
JF, (ΞΞ m+2.5._ 5)- Gye Ff 1 
2. 3 x? nt dx" L/(x? +f 
(7) Establish the following formule: 
(ὃ 2Fi(@; Bs a+B+e; x) x 2F,(y; δ; y+5—Q; x) 
= 2F,(a+0,B+0;4+8+0;x) x 2Fi(y—@, 5-0; y+ 5-03); 
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(ii) Dx{x* F(a, B; y; kx)} = re πὶ 2F,(a+1; B; y; kx), 
where «0; 


(iit) BCA, y—A) 2 F(a, B; γ; x) 
= [ t*-*(1—1)’"*-! , F(a, B; A; xddt, 
0 


where | x |<1, A>0, y—A>0. 
(8) Prove that if f>0, 
2F (4, B; 2B; 2) 
; Gs ap-1| {1+0cos p}-*] | 
276-1B(B, B) Jo “@ wks Pre cos $} SA 49, 


where ζ = z(2—2). 
Deduce that 


2F (a, B; 2B; 2) = (1 -- 22), 52 Ga, 4a +4; B+4; C’). 
(9) Prove that 
| " cos m@ cos" θ ἐθ = eer. ἢ eee 
0 2"**T(4n+4m 4+ 10 (4n—4m +1) 
and evaluate 


ἀπ 
| cos 2mé@ sin" 6 dé 
ἢ 


where m is a positive integer. 


(10) Schrédinger’s equation for the rotation of a sym- 
metrical-top molecule is 


ee ot) 4 1 ay 
sin θ 80 80) ᾿ sin? 6 802 


δ'ψν 2.ο50θ d*~ | 8n7AW 
+ | cot? 04+ —)—— — ———— τ = = 0, 
( 4) dy? _— sin? 8 dyad hea: ha 4 
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where A, C, W, A are constants. Show that it possesses 
solutions of the form 

ya eee — λα g8O-™) (0, B; 75 2), 
where n=m, x = 4(1—cos 0), y = n—m+1, and a, f are 
the roots of the equation 


A 827AW 
2_(2n+1 —n*+n— = 0, 
z*—(2n+1)z+ ra n 3 


(11) Prove that: 
(i) ,F,(@; a; x) =e; 


(ii) ,F\(a+1; α; x) = (1+ =) δ! 

Gi) ,F,0; 4; -x)= VE ert (x) 

(iv) ,Fy\(a+1; 7; x)—, F(a; y; x) = F(a 1374152); 
(v) {δι (- 53; —x?) = ε  -- πα erf (x); 

(vi) x", F,(n; n+1; "τα" “τὶ ett. 


0 


(12) Prove that the equation 


5 Mae Md 
ὃχ2 k at 
possesses solutions of the type 
V=Ct",F 4 “| 
= a | 
( 4k 


m where and C are constants. 


HYPERGEOMETRIC FUNCTIONS 47 
(13) Show that the Schrédinger equation 
Vet (- t) w=0 
r 
possesses a solution of the form 


e™* Ρ, (- Ms 1: ikr—ikz 


(14) The Schrédinger equation governing the radial 
poe functions for positive energy states in a Coulomb 
eld is 


1d/(_aL\ , [8x?m zz'e*\ πίπ- 1) 
ee) * be ‘ain er Jeo 


Show that it possesses a solution 
L=r'e™ ,F,(ia+n+1; 2n+2; —2ikr) 
where k? = 8n?mW/h?, « = 4n?mzz'e?/k. 
(15) Show that the equation 
2 2 
possesses a solution 
y=xe"* ~, F,(4n+4—if; 2imx) 
and hence that a solution of equation (1.9c) is 
R = οἷ e7*F(in+4; n+1; 2imo). 
(16) Show that 


ἀῶ ν΄, ἀνα 
᾿ eet 


oo I: kK | 
= B(1, m) p44F | pital ἐκ } 
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ΟΞ 
ὦ [ x χ ἢ -ἀὐ nf Te " 
5: 


1s ** 


| αι, ony Ley 25% 
=B 1, : F 1 ἢ ΞΡ ᾿ . 
( M) p+1 41 ie a Bos nl 


, 1-x 
ai) | (1--χΧ fe ee By | 
» Be 


Ls s+ 
les. nang ths ee 
= BS, m) p41 ἔχει pa a B., 2m; } 


(17) Prove that 


(i) |". F| os Pe erated 
is «+» Pas | 


* Pas <0 Ba’ 


ΓῚ 
i? oe | ἢ} 


oo ἢ ya 
(iii) | Pf i: oe; re +O" 1 ὰ 
P en 


= T(4p) Ε ἴᾳ rong Ony $y; dail 
= “a eo Pri q i 
Lis 


2p" “cng Mas 
(18) By equating coefficients of x in the relation 
(1—x)*~?-? F(a, Bs γ; x) = 2Fi(y—% γ-- β; γ; Χ) 


Tw or i. tapi He} 
a’ ae ae 


TW) OF μι τὸ ἂρ ἐμ, utd ~~ 
p pt2*¢ Γ 
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prove Saalschutz’s theorem 
JF, in B, —n; 1 |- (γ,--αἡ,()-- β), 
, Ys 1+a+B—y— (y),(9— —a—f), 


Hence prove that 


LF, gl =(1-x)"*,F, - ἀμ 1 a 


1+a—f; 
if | x |<3—2V2. 
(19) Show that 
2F (a, B; 1+a—; 2) 
Ξε (1 -- 2) 52} (ῥα, 4+4a—B; 1+a—B; ἢ 
where ζ = —4z(1—z)~?. Hence deduce the value of 
2F, (a, B; 1+a—B; —1) from Gauss’s theorem. 


(20) Show that 
a, B, ee 1 
ὦ | ὃ, 8; 
_ TOO) ς bone &—a, σὶ ‘ 
Γ(Γ(σ Ἐ β)Γ(σ-Εγ}) *Lo+B, o+y 


where σ = δτε--α-- β--γ. 
Hence, using Dixon’s theorem, prove Watson’s theorem: 


F, * B, 93 1 | 
ΠΣ 21 +0+ B), 27; 


_ T@IG+yTG+44+3A0G—4e—-F8 +7) 

P($+3o)°G+4AFG—40+y)G—FB +7) 

and, using Watson’s theorem, deduce Whipple’s theorem 
that, ifa+f = 1, ande+6 = 2y+1, 


3 ἢ ᾧ τ Ἐΐ ‘| 
nT (O)I(e) 


~ PIT Gat {OM Gat sO GR +45 GP +40). 


CHAPTER III 


LEGENDRE FUNCTIONS 


§ 13. Legendre polynomials. If A is a fixed point with 
coordinates (a, β, y) and P is the variable point (x, y, 2), 
then if we denote the distance AP by R, we have 


R? = (x—a)’+(y—B)*? +(z—)’. 
Furthermore, we know from elementary considerations that 
1 


R 


is the gravitational potential at the point P due to a unit 
situated at the point A, and that this must be a 
particular solution of Laplace’s equation. 

In some circumstances it is desirable to expand yw in 
powers of r or r~* where r = (x?+y*+2z7)? is the distance 
of P from O, the origin of coordinates. This expansion 
can be obtained by the use of Taylor’s theorem for func- 
tions of three variables but it is much more suitable to 
introduce the angle @ between the directions OA, OP (cf. 
fig. 4) and write 

R? = r?+a*—2ar cos θ. 


The expression for yw then becomes 
abl a3 laa 
J (a? —2arp+r’) 


where p: denotes cos 0, and this can be expanded in powers 
of γα when r<a and in powers of a/r when r>a. If we 


(13.1) 


{ 


qu θπσν--------........................ ὁ “.ὕ 
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denote by P,(u) the coefficient of h" in the expansion of 
(1—2yh+h?)~* in ascending powers of ἡ, ice. if 
1 oo 
= P(w)h", 13.2 
/(1—2ph +h?) 2 ) ane 
then the potential function (13.1) can be expanded in 


P 
Τ Ε 
εη΄Ζ 
0 a A 
Fig. 4 
the forms 
4 3 4 Pu) r<a; (13.3a) 
da=0\a/ " Γ ; 
lL S&S far PL, 
- ὦν (2) (μ), γΡα. (13,30) 


It is clear from the definition (13.2) that the coefficients 
P,(u) are polynomials in μ, The first one or two can 
readily be calculated directly from the definition. By the 
binomial theorem we have 


(1—2yh+h?)-* | 
= 14+(—4)(—2ph+h?)+ (ΣΧ Ὁ. yh b+... 


= 14+ ph+4(3p? —1)h?+4(5p3 —3y)h? +... 
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so that 


P(w=1, Pio) =4, 

P(u) = 4(3p?—1), P3(u) = 4(Su°—3y). (13.4a) 
We shall show below that, in the general case P,(j) is a 
polynomial in μ of degree n; it is called the Legendre 
polynomial of order π. 

The expression for the general polynomial P,(y) can 
be derived by the method employed to obtain the simple 
expressions (13.4a). 

Expanding (1—2yh+h?)~* by the binomial theorem 
we have 

—— 2\-+ — . (3), ων 
(1—2ph+h’)*#= Σ τ (ἐμὴ h*) 
r=0 Pr 


and the coefficient of h” in this expansion is the coefficient 
of h" in the expansion 


y 3 (3), + Quh— h*s = ." ee (2uh— —h?y"~@ 
ye ὧν (μῆ-- hye 
“See eo 
since by example 9 (ii) of Chapter I 


(- .6),-.- ἴθ. 
° G-n), 
Now the coefficient of 4” in the expansion of 


= Ὁ (2h —h2y'~¢ 
—@)! 


ὍΝ * 
o'(n—20)! 
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and, by the duplication formula for the gamma function, 


πὶ -28 
(n—2g)! opener he 
. τ ΤΕ ων. 
Γ(π Ἐ32--οὐ Γ(π Ἐ1 .--ο) (ἐ -- ἀπ). (--ὄ ὅδ}, 
so that 


Ὃ oh n τ' G—4n){—4n), +) 
Ρι(μ) = ee) 4 οἷ -- n), (4 . 


a result which may be written in the form 
2 , 1 
Pn) = CHD se, (4am, —4nsd—m 5). (3.4 
Putting «= 1 in equation (13.2) and equating co- 
efficients of h” we find that 


P(1)=1 (13.5a) 


for all values of ἡ". Similarly if we put « = —1 in (13.2) 
we derive the result 


P{—1) =(—1) (13.5b) 
which is a particular case of the result 
P,(—p) = (—1)"P,(#). (13.6) 


Equation (13.4) gives P,(cos 0) as a polynomial in cos 0 
of degree n so that it should be possible to express 
P,(cos @) in terms of cosines of multiples of 0. Instead 
of attempting to do this by substituting the appropriate 
expression for cos’ @ in (13.4) we begin afresh with the 
definition (13.2). Writing (1—2 cos @h+h7) in the form 
(1 —e!®h)(1 —e™ !*h) we find that 


y P,(cos 0h" = (1—he®)~*(1—he~*)-# 


- I(r +4)I(s+4) ht tspl(r—s0) 
r=0s=0 ΓΟ ΓΑ)» Is! 
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Equating the coefficients of h” we find that 
> Γ(Σ ἘΡΓ(Σ ῬΗ -- τῇ yar-no 
P.(cos @) = ~ = $e 
sna % P@)r@yria—r)! 
Using the duplication formula, we see that 
Pg+nrgtn—r)_ 1 (2n—2r)\(2r)! 
Γ(Ω ΓΟ) 225. ri(n—r)! ᾿ 
so that 
Σ (2n—2r)\(2r)! οἰ(ῶν -- πϑ 
a, =o (rl)? {(n— I? © 
from which it follows immediately that 
P,,(cos @) 
ὃ {(2n)!}? | 1 "ἢ (n—2r)\(2r)! 
2*"(n!)* 241,5 (r!)?{(2n—r)!}? 


P,(cos θ) = 


cos(2n—2r), 


᾿ (13.7) 
an 
P2544 (COS 0) 
{ 4 (4n+2—2r)!(2r)! 
2 πε ἢ ap Gari nye on nar +e (13.8) 


From these last two equations we may derive a general 
result of some importance. We may write 


Piva) = Σ ς, οοβ(η -- 2))0, (13.9) 


where p = 4n or 4n—4 according as n is even or odd, In 
particular 


ΡΞ 7. δ,» 


Ρ (cos θ) 
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Fig. 5. Variation of P,() with p. 


| P,(cos 0) ΞΞ 1. (18.10) 


0-5 


Fig. 6. Variation of P,(cos 6) with 0. 
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so that from (13.5) 
P 
iw Ὑ ὁ, 


γιῷ 
Now, from equation (13.9) we have 
P 
| P,(cos 0)|S Σ᾽ ο;. 
r=0 
and therefore 
| P,(cos θ) 3 1. (13.10) 


The variation of P,() with » for a few values of n is 
shown in fig. 5. Since, in most physical problems, the 
Legendre polynomial involved is usually P,(cos 0) we have 
shown in fig. 6 the variation of this function with 0. 
Numerical values may be obtained from Tables of Associated 
Legendre Functions (Columbia University Press, 1945). 


ὃ 14, Recurrence relations for the Legendre polynomials. 
If we differentiate both sides of equation (13.2) with respect 
to h we have 


ἘΠΕῚ; ὟΝ = Σ nh"! P(y), Ιμ [«!, 
(1—2ph+h?)? =o | 


which may be written in the form 
(u—h) 3 h"P,(u) = (1—2uh+h?) ἊΣ nh™~*P(u). (14.1) 
n=0 n=0 


Equating coefficients of ἢ" we have 


uP, (pu) — P,-s(4) 
= (n+1)P,4 (4) —2nuP,(u)+(n—1)P,- (0), 

which reduces to 
(n+1)P,+;(4)—(2n+ 1)eP,(u)+nP,- (4) =0. (14.2) 
This relation has been proved to hold for |u| <1 but 


a 
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since the left-hand side is a polynomial in p, it must hold 
for all values of p. 

On the other hand, if we differentiate both sides of 
equation (13.2) with respect to 4 we obtain the relation 


h —— 
(1—2ph+h?)? — 2h +h = 2, h®P*(n). (14,3) 


Combining equations (14.1) and (14.3), we have 


(u—h) 2. μπΡι(μ) = Σ, nh"P,(), 
A = A= 
so that equating the coefficients of ἢ" we obtain the relation 


nH) — P,- χ(μ) = nP,(u), (14.4) 
and since each side is a polynomial in p this relation holds 
for all values of μ. 

If now we differentiate equation (14.2) with respect to 
μ we obtain the relation 
(n+1)Py41(H)—(2n+ IP, (4) 
—(2n+ DpPi(u)+nP,_(u) =0. (14.5) 


Eliminating P,() from (14.4) and (14.5) we see that 


Pi+1(H) — Ῥὰ- ((μ) = (2n+ 1)P,(u). (14.6) 

Subtracting (14.4) from (14.6) we obtain the recurrence 
relation 

Pot ι(μ) — μΡι(μ) = (n+ 1)P,(u). (14. 7) 


The differentiations with respect to A and mw under the 
summation sign is justified by the fact that the series on 
the right-hand side of equation (13.2) is uniformly con- 
vergent for all real or complex values of ἢ and μ which 
satisfy the relation | 4 |<1, Ph |<./2-1. 


᾿- 15, The formulae of Murphy and Rodrigues. From 
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the expansion (13.2) it follows immediately that 


~ d’ ἡ" ΤῊΝ 
ht — = -[1.. 28. 2) 
eo age ign tah ΕἸ) 


= 2h’ ———= (1—2uh+h 3, 
(4) en 


Substituting the value μι = 1 we see that 


. pe pn = or LOD h'(1 —h)y7 2th 
Σ ρθῶ = 2 HED yam 


= 2" P(r+3) hr τ Mi+2r+s) ἢ" 
ΓᾺ) s=o F(1+2r)s! 
where ΡΟ) denotes d"P,(u)/dy’. 


Equating coefficients of h" we see that P6°(1) = 0 if r>n, 
as is obvious from the fact that Ρ. (μ) is a polynomial of 
degree n in p, and that 


T(r+4) T(i+n+r) 
ra) Γ(Ἐ 2ν)(α -- υ) ἢ 
From the duplication formula for the gamma function 
etd) 5 Ae k, 
T()r(i+2r) or !2" (1),2” 
and from example 9 (iii) of Chapter I, 


POM) = 2 


PUSRE 2 (174 =n), 
(n—r)! 


so that 


POL) = _ yy @+DA—™), ) . 
I > τὸν (15.1) 
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Now, by Taylor’s theorem 
P= Σ SO poy, 
r=0 rt 


Substituting the expression (15.1) in this expansion we 
obtain the relation 


_  (—n)(n+1), {-πιΚ 
P= 2. avr Pst) 


which gives Murphy’s formula, 


for the Legendre polynomial P,,(). 
If now we put ἃ = | in example 5 of Chapter II we see 
that equation (15.2) is equivalent to 


2a 7 

Pu) = —— — (p*—1)", (15.3) 

2"n! du 

which is Rodrigues’ formula for the Legendre polynomial. 

Rodrigues’ formula is of great use in the evaluation 

of definite integrals involving Legendre polynomials. 
Consider, for instance, the integral 


1 
I= | S(x)P,(x)dx. (15.4) 
-1 
By Rodrigues’ formula we may write this integral as 


| fy & (x2 -17"dx, 
ie dx" 


2"η! 


and an integration by parts gives 


== Ἐπ]. f'(x) 6 τ᾿ {(x? —1)"}dx 
2"n! J_y ent 
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The square bracket vanishes at both limits so that we have 


tags [ γὼ 2 2? -1y}ax 
Sal a | 
Continuing this process we find that 
a& “anh (2—1'f%)dx. (15.5) 


2"n! 
For example if f(x) = P,,(x), m<n, f(x) =0 and so 
J=0. In other words 


| TRAP «ὦ tae αἰτῶ 
-1 
If fix) = P,(x) then 
oy) eee ap 
fe) = 5 -1y 
(2n)! 
2"n! 
Hence : 
2n)! : 
__ PGP ax = ae (1—x?)"dx 


_ 2)! TRT@+)) 
P(n!)? T(n+3) 
Making use of the duplication formula for the gamma 
function we can reduce this to the form 
ὃ 2 
P{x)}*dx = ——. 15.7 
[- (paeoydx = Σς (15.7 


A convenient way of combining the results (15.6) and 
(15.7) is to write 


2 
[ Pa(a)Pyxddx= 7-5, (158) 
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where 6,,,, is the Kronecker delta which takes the value 
0 if m # n and the value | if m = n. 
Similarly if f(x) = x”, where m is a positive integer, then 
T(m+ 1) Ὗ 
---΄ὦοὃ- “ x™"", if) 
f(x) = Γῴὼω-τη Ἑ Ὁ) 


0, if m<n, 


n2n, 


and hence, if m>vn, 


| ᾿ x"P(x)dx = eet x™-"(1 —x?)"dx. 
=i =< 


2°T(m—n+1)n! 


If m—n is an odd integer the integral on the right is zero 
while if »2—7 is an even integer it has the value 


r(4m—4n+h)0(n+1) 


1 
ΞΗΓῚ . νῷ = 
2} x" "(1 —x*)"dx rGm+in+3) 


so that, if m is an integer, 


Ι 
| x™P(x)dx 
=] 
0, if m<n, 


- __miPGm—int+d) if m—n2O is even, (15.9) 
2"(m—n)!0(4m +4n+4+ 3) 


0, if m—n>0 is odd. 


If » = n the result is 
1 1 1 
| x"P.(x)dx = — | (1—x*)" dx 
-1 2” J-1 


ΓΦ ΓΙ ἢ) 
25 T(n+3) 
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which, on account of the duplication formula, is equivalent 


to 1 P a 2"*1 (nt)? 
ΕᾺ om ΤῊΝ 


§ 16. Series of Legendre polynomials. In certain 
problems of potential theory it is desirable to be able to 
express a given function in the form of a series of Legendre 
polynomials. We can readily show that this is possible 
in the case in which the given function is a simple poly- 
nomial, For example, from the equations (13.4a) we have 

1 = Po(u), 
= P τ(μ), 
μὴ = 44+3P2(u) = $Po(u)+3P2(u), 
μ᾽ = 3u+FP3(u) = 3Py(u) +3P3(H), 
so that any cubic cou +c,u?+c,4+c3 can be written as 
the series 


(15.10) 


2c, 


3 Ρ;(μ)Ὲ (= +2) Px (7 τος) Po(p). 


It is obvious that we could proceed in this way for a 
polynomial of any given degree ἡ though if m were large 
the arithmetic involved might become cumbersome. Since 
Ρ(μ) is a polynomial of degree n in p, it emerges as a 
result of an extension of the above argument that any 
polynomial of degree » in 4 can be expressed as a series 
of the type 


me Ρι(μ)ν 


Σ ¢P(u), (—iSpsl). (16.1) 


The problem which now arises is that of expressing any 
function f(z), defined in the interval —1<yu<1, as a series 
of Legendre functions of the form 


Ss Pu). (16.2) 


r= 0 
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If it is assumed that the infinite series (16.2) converges 
uniformly in the range (—1, 1) to the sum f(u), we may 
multiply the terms of the series by P,(4) and integrate 
term by term with respect to 4 over the range (—1, 1) 
to obtain the relation 


[ ΤΡ, (μ)άμ = ᾿ C, [ Ρ(Ρ,(μ)άμ 
' - 1 os uf 


and by equation (15.8) the sum on the right-hand side is 
equal to 


= 2 2c, 
ἔ,----- ὕ, κα ™ : 
“δι 2n+1 ° 2n+1] 


which shows that the series 


oO Ι 
2, (n +P | 7}, (ν)άν (16.3) 
r= -1 

converges uniformly to the sum f(z) in the range (—1, 1). 
The series (16.3) is called the Legendre series of the 
function f(z). We shall not discuss here the conditions 
which must be satisfied by the function Χμ) if this series 
is to be uniformly convergent; for such a discussion the 
reader is referred to Chapter VII of E. W. Hobson, The 
Theory of Spherical and Ellipsoidal Harmonics (Cambridge 
University Press, 1931). 

The possibility of expanding a function in the form of a 
series of type (16.2) is a consequence of the relation (15.6). 
In the theory of special functions we frequently encounter 
sequences of functions ¢,(x), $2(x), ..., %,(x), ... which 
have the property 


[ φ,, (χ)φ,(χ)γάχ =0, (mn). (16.4) 


We then say that the functions ¢,(x), (r = 1, 2, ...), form 
an orthogonal sequence for the interval (a, 5). If, in 
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addition, the functions are such that 


[ {¢,(x)}7dx = |, (16.5 ) 


for all values of n, we say that the functions of the sequence 
are normalised, and form an orthonormal set. Given a 
set of orthogonal functions it is obviously a simple matter to 
construct a normalised sequence. For example we see from 
equation (15.8) that the sequence of functions P,(x), 
(n = 0, 1, 2, ...) is orthogonal but not normalised. By 
multiplying each function by ,/(n+4) we find that the 


functions 
| (n+4)*P,(x) (16.6) 
form a sequence of normalised orthogonal functions in 
the interval (—1, 1). 
There is another property of importance which such a 
sequence of functions may possess. If there is no integrable 
function W(x), different from zero, such that 


b 
ψ(χ)φ, (χ)άχ = 0 (16.7) 


for all values of n we say that the sequence is a complete 
orthogonal sequence. It may be shown in the case of 
the functions (16.6) ¢ by considering the Fourier coefficients 
of W(x/n)P,(x/2) in (—z, 2) that if (16.7) holds this function 
is a null-function and hence that W(x) is a null-function 
in (—1, 1). In other words it can be shown that the 
functions (16.6) form a complete sequence of normalised 
orthogonal functions. 


ἢ 17. Legendre’s differential equation. If we write 


υ = (u?—1)" 
then it is readily shown that 


(1 -- μῦ) ue +2punv = 0 
du | 


+t E. W. Hobson, op. cit., p. 40. 


§ 17 LEGENDRE FUNCTIONS 65 


and if wé differentiate this equation n+1 times using 
Leibnitz’s theorem we find that 


d"**y τυ d"v 
1—p*) —— —2u —— +n(n+1) — =0 
( μ ΤΩΣ Ya ( dae 


which when written in the form 
d? d d"v 
1—p?) — —2p — +n(n+1)}| —] =0 
{ a a } @) 
shows that d"(u?—1)"/du" is a solution of the differential 
equation 
12 
(1—p?) mae <n 2 nine iy = 6 (17.1) 
ἐμ ἀμ 
so that we conclude from Rodrigues’ formula (15.3) that 
when ἢ is an integer P,(y) is one solution of the equation 
(17.1). This equation, which we shall now consider in a 
little more detail, is called Legendre’s differential equation. 
We saw in example 1 of Chapter I how such an equation 
arises in the solution of Laplace’s equation when solutions 
of the type R(r)O(cos 0), 1.6. m = 0, are sought. 

It is obvious by inspection that the point μ = O 15 an 
ordinary point of the equation (17.1). Writing the equation 
in the form 

d*y 2u dy n(n+1) 
—1)? —= +(n—1) ——— — = (p—-Dy = 0 

(μ-- 1) 42 ( Pree al (u— Dy 
and observing that in the notation of equation (3.1) with 
a= il, 

1+(p—1) 2 
en et Pd 1) — tigi)... 
Ρῷ) τί -ἢ Hu—1)—Hu—1) 


au) = —An(n +1) A 


= —jn(nt+- {(u-1)—HU-1)? +...}, 
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we see that jt = | is a regular singular point with indicial 
equation οἷ = 0. 
Furthermore in the notation of equation (4.2) 
—2 n(n+1 
at Buy = MED, 


α(μ) = 
so that as u- 00 


α(μ).- 7 piy)~ — Med. 


showing that the point μ = οὐ is a regular singular point 
with indicial equation {g—(n+1)} (e+) =0. We thus 
see that the equation is defined by the scheme 


—j wo 1 | 
y=P{ Ὁ ntl O 4g}. (17.2) 
0 -n 0 
If, however, we put 
x = ΚΙ --οὖ 
in equation (17. : we find that it reduces to the form 


x(1—x) τὰ ots 2x) 2 hata lyy 0, (17.3) 


which is oie (8.1) Dong a=n+l, B= —n and 
y = I, so that the scheme (17.2) is equivalent to the scheme 


So ο 1 
y=P;0 n+l 0 }2-]. (17.4) 
0 -" 0 
It should be noticed that the values along the top row are 
those assumed by 4—4y, not by μ. 
We consider first the solution corresponding to the 
Singular point at infinity. We write 


yi(H) = μ' ᾿ Cyn” 


yi(4) =e — 
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which on substitution into (17.1) leads to the recurrence 
relation 


(n—v+2)\(n—v+1)c,_, = —v(2n+1—y)ec,. 
On taking co = 1 we obtain the solution 
m(n=1)_n-2., Mn—1)(n—2)(n—3) στὰς 
2.(n—1)’ 2.4. κ- ὕχπ.-3) νὴ 
.}.. (—39)G-3n) 1 
-- 41: -----:--ς-ς- 
ἱ 1.4 - κ' 


+ (dint NGI G—3n+1) 1 +4 
1.2(4—n)($—n+1) μ' 


which may be written in the form 
: 1 3 
yi) = HF (-¥ 4— jn; 4-2; 3) (17.5) 
Also, if we write for the second solution 


yw) =e" 2, ἀμ" 


we find that 
= port, @EDNt2) | -n-3 
yu) Ξ μ΄ + 2.(2n43) μ 
(n+ 1). |5Ἐ2.(π} 3Χ,π4) .-.- 5. 
2.4.(2n+3)(2n+5) 
γε: EAR (an-+4 n+1; n+; +) (17.6) 
μ μ 


provided m is any number other than a negative integer 
or half a negative integer. 

These solutions are valid for a// values of n for which 
the ,F, series have a meaning, not only for integral values 
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of π. If ἢ is an integer the series y,(u) terminates—in 
other words, y,() is a polynomial of degree m in μ. If we 
multiply this polynomial by 

(2n)! 

2"(n!)? 
we obtain the Legendre polynomial of degree n (equation 
(13.4b) above). 

On the other hand the series for y,() does not terminate 

when n>—1 so there is no point in restricting n to be 


an integer. This series solution when multiplied by a 
factor 


P(3)0 (n+ 1) 
21 (n+3) 
gives the function 
ο,(μ) 
οὐ ΝΡ -,- | οὐ ee 
= ΕΙΣ Ω oF, (4 +4, 4n+1; n+; +} (17.7) 
The function 


_ saath . | ‘ 1 
PAM) = IMT (n+ 1)}2 Κ 2; (-¥ 4—4n;4-n; )(78) 


is a solution of the Legendre equation (17.4) even when 
nis not an integer and it reduces to the Legendre polynomial 
when 7 is an integer. We shall continue to denote it by 
P,(u) but when ἢ is not an integer shall refer to it as the 
Legendre function of the first kind of degree m. The 
function Q,(4) defined by equation (17.7) will be referred 
to as the Legendre function of the second kind of degree πὶ 
even when ἢ is an integer it is not a polynomial. 

With these definitions we may write the solution of 
Legendre’s equation (17.1) as 


y = AP,(u) + Β0,(μ) | (17.9) 
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(| |>1). In some problems we know that the solution 
should be a polynomial in μ; in that case we must take the 
solution to the form y = AP,(y). 

The variation of Ο,(μ) with μ may be computed easily 
from equation (17.7) when μ᾿]. Tables calculated in this 
way are contained in the volume quoted at the end of 
8. 13. The following fig. 7, which was prepared from these 

Ξε: ᾿ sal ery 


0-4 


—_—t— 


Fig. 7. Variation of Q,() with μ. 


tables, shows the variation of O,() with μι for a few values 
of n. 

Since Legendre’s equation has a regular singular point 
at 4 = 1 we may on the basis of equation (3.8) take the 
second solution of Legendre’s equation to be proportional to 


Ρι(μ) log (u—1) ἘΣοίμ-- 1. 


The coefficients c, can now be obtained by substituting 
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this expression into the differential equation (17.1) and 
equating to zero coefficients of successive powers of 
(4-1). Instead of proceeding in this way we shall derive 
this second solution by means of a method due to F. E. 
Neumann. 


§ 18. Neumann’s formula for the Legendre functions. 
Let us now consider the integral 
{ P,(@)dé 
=-1 μα 
where | |>1, and ἡ is a positive integer. Expanding 
the denominator by the binomial theorem, we have, for 
the value of the integral, the series 


a 1 1 ᾿ ᾿ 
2, aa | ξ Ρ (ξ)άξ. 


From equation (15.9), it follows that the integrals in this 
series are zero if sSn, or if s = n+2r where r is a positive 
integer, so that the above series is equivalent to 


oo | 1 
Σ μῆτε] or@e 
which, by the same formula, is equal to 
L$ (nt 24d) (1 Ἱ 
μπ τ᾽ πο 2.2) Γ(η -Ἐσ 3) λμ2)᾽ 
Now from the duplication formula 
(n+2r)G+r) _ TGn+3+r)lGn+1+r) 
2"(2r)! r! ἢ 
so that the series reduces to 
1. 8 MGn+34+nPGnt+i+n (4) 


pet? » So C(n+3+r)r! μ 
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which is equal to 
PGn+d0Gn+1) , 5 
Ta+pe oF; (an +4, 4n+ 1 ἢ n+4; +) 
Noticing that 
Pgn+d0Gn4+ D/T(n+9 = TD (n+ D/2'T(n+-®), 


and comparing the result with equation (17.7) we see 
that this series is merely 20,(4). Hence we have shown 
that if | »|>1 


1 
Ο,(μ) = | PAS) ge, (18.1) 
«ἰὸς 


a result which is known as Neumann’s formula. Equation 
(18.1) holds not only for real values of y greater than 1, 
but for all values of 4 which are not real. For this reason 
equation (18.1) may be regarded as defining the second 
solution, Q,(u), of Legendre’s equation. 

Certain related formulae, due to MacRobert, follow 
readily from this result. If | «|>1 and m is a positive 
integer then 


1 em Ι Ὦ ἘΠῚ 
wro,un—3 | EBD ἐς -α [" PE ρα 
uy prs “ὦ f-—¢ 
and the integral on the right is equivalent to the finite sum 
m=-1 i 
ΣῈ με τ | δΡ (δας. (18.2) 
Ρ- -1 


If msn it follows from equation (15.9) that each term of 
this series vanishes so that we have 


Ι ΠΊ 
μ"Ὁ ,(μ)- 4 | gh (s) dé (18.3) 
=, pe 


provided m, n are integers and mn. 
On the other hand, ifm = n+1 the series (18.2) reduces 
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to the single term 


— _ 2(n!)? 
i Α eP(Oag = FO 


ee ore -Ο TO, 
e+ μὲς of (2n+1)! ae 


Now if m is an integer ah is a polynomial of degree 
m in μ 80 that it follows from (18.3) by finite summation 


that, if m<n, 
P,(1i)Q,(s1) = 4 | OPES ἀξ. (18.5) 
cal | 


Similarly from equation (18.4) and the definition of Ρ,, χ(μ) 
we have the formula 


ΡΥ = τ “ΚΟ ΚΟ δι 1. (80 


so that 
n"**0,() = 3, 


If we replace n in equation ἰὸν - by +1 and m by ἡ 
and subtract from equation (18.6) we obtain the relation 


Pas ι(μ)0,.(μ) — P,(0)0,, + ι(μ) = —~. (18.7) 


Other formule of a similar nature are contained in ex, 24 
below. 

We shall now make use of Neumann’s formula to derive 
the form of the second solution of Legendre’s equation in 
the neighbourhood of the points 4 = +1. From (18.1) 
we have the result that if p>1, 


Q,(1i) = 4P,(ui) log At -W,-(u), (88) 


where W,,_ ,(u) denotes the integral 
ἢ Le Mo PA "ὦ dé. 
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Now when a is an integer, P,() is a polynomial of degree 
nin μ so that {P,(u)—P,(6)}/(u—€) is a polynomial of 
degree n—1 in μι. Hence W,_,(u) is a polynomial of 
degree n—1 in μ. 
If we substitute from (18.8) into Legendre’s equation 
(17.1) we find that W,_, satisfies the differential equation 
(1—p?) Wy, — 2H, +n(n+1)W,_-, = 2P,(u). (18.9) 


Now since W,_,() is a polynomial of degree n—1 in p 
we may write (cf. § 16 above) 


W,-1(4) = ΣῈ c,P,(u). 
Using the fact that 
Ρ 
Ρι(μ) ἊΣ 2 (2η -- 4: -- 1)...-- )ν-τ αν 


where p = }(n—1) or 4n—1 according as nm is odd or 
even (which follows from equation (14.6)), and the result 
— WYP i(W) — 2uP (pw) + n(n + 1)Ρ,(μ) 
=(n—r)(n+r+1)Pi(), 
we find, from equation (18.9), that c,_,, = 0, (s = 0, ..., p) 
and that 
2n—4s—1 
Gt μόρα, ται, 
(2s+1)(n—s) 
Substituting these values in (18.8) we obtain the formula 


P,(u) | με1 ὦ eS, ἐρεῖ 3- ; 
Q,(u) = ΣΡ, (μ) tea πε Xo τε α- ὦ ' cm 


Another expression for Q,(u) may be derived from the 
fact that both P,(u) and Q,() are solutions of Legendre’s 
equation (17.1) so that 


On(4) < {(1=12)PHCw)}— Paw) — ((.-- μοι} = 0 
, ἐμ μ 
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which is equivalent to 


᾿ [(-- nO, P(t) — P,(u)04(1)}] = 0, 
showing that 
(1-1), PH) —PrWOKw)}=C, (8.1) 


where C is a constant. Now from equations (17.7) and 
(17.8), we can readily show that for large values of μ 


s ! 1 n+1 1 
PHCWOQ, τος κ(μ},(μ).- -- 2 
MD QWW~ τον - 12; QIU) — FT + 
so that as poo the left-hand side of (18.11) tends to 
—1, showing that C = —1. Writing (18.11) in the form 


a {ae it 
du (Pw) (μ"--1)(Ρ,(μ}}} 
and noting from (17.7) that 0,()-0 as u-0o we have 


| Pe dé 
0,(u) = P, = ST ---- 8.12 
ω ω ΞΡ; Ὁ 


ὃ 19. Recurrence relations for the function Ὁ (μ). 
Recurrence relations for the Legendre function of the 
second kind can be derived from Neumann’s formula (18.1) 
and the corresponding recurrence relations for the Legendre 
polynomials P,(u), From the recurrence relation (14.2) 
and Neumann’s formula we have 


(n+1)Q,41(u) +n0,_1() = (n-+4) [ AD ἀξ. 
=] ae 
Now 


+ | SPAS) ae = μο,(-- ἢ [ P(E)dé. 
Jay μ’--ξ a 
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1 
If we write the second term on the right as | P(E) Ρ, (ξ) ἀξ 
-1 
we see from (15.6) that it vanishes if nm #0. Hence we 
have 
(n+1)Q,41(u)—(2n+ 1)pO,(s)+n2,-()=9 (19.1) 


showing that the functions Q,(u) for three consecutive 
values of πὶ satisfy a relation of the same form as that for 
the functions P,(4:) (equation (14.2) above). 

From Neumann’s formula (18.1) we have 


Ἢ 
un)= -4 |) 78a 


and if we integrate by parts on the right-hand side we 


find that ; ϑ 
Ve © Le Cal! 6 Pu) ge. 
un 47. + SU} 4g |" Bag 

Hence 

1 , _. pt 
O5+s(8)— 91-00 = 3 | ἐπ ΟΞ ἐμ iO) a 
“ΔΝ. Ὁ | i - dé, 
-4 i- 


by virtue of equation (14.6). The integral on the right is 
20 (μὴ by Neumann’s formula so that, finally, 


O14 (Hu) — ΟἹ... (μὴ) = (2n +10, (0). (19.2) 


§ 20. The use of Legendre functions in potential theory. 
In potential theory we have frequently to determine 
solutions of Laplace’s equation V*y = 0 which satisfy 
certain prescribed boundary conditions. If we have a 
problem in which the natural boundaries are spheres with 
centre at the origin of coordinates it is natural to employ 
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polar coordinates r, 0, ¢. In cases in which there is sym- 
metry about the polar axis Ψ will not depend on ¢ so we 
may write y = wW(r, 8). It then follows from example 1 
of Chapter I that 

(A,r"+B,r~"~*)o, 


will be a solution of Laplace’s equation provided that », 
is a solution of Legendre’s equation (17.1). Taking v, to 
be P,(cos 0)+C,,Q,(cos 0) we see that we may write 


bir, δὴ Σ (A,r"+B,r-"~*)P.(cos 8) 


+ Σ, (C,r"+D,r~""")Q,(cos 0) (20.1) 
where the quantities A,, B,, C,, D,(n = 0, 1, 2, ...) are all 
constants. 

Now it is obvious from equation (18.8) that Q,(cos 0) 
is infinite when 0 = 0°, and we know on physical grounds 
that in the case of spherical boundaries y remains finite 
along the axis 0 = 0. Hence we must take C, = ἢ, = 0 
for all values of n and obtain the potential function 


Y= La (A,r"+B,r—"—*)P,(cos 0) (20.2) 


which is valid as long as r is neither zero nor infinite, i.e. 
if aSrsb where a and ἡ are finite and non-zero. If the 
region under discussion is the interior of a sphere, i.e. if 
Osrsa, then to avoid y becoming infinite we must take 
B,, to be zero to give 


‘. Bi A,7"P,(cos 6). (20.3) 


On the other hand, if the region being considered lies entirely 
outside this sphere, we must take 


y= >. B,r~"~* P.(cos 6). (20.4) 


§ 20 LEGENDRE FUNCTIONS 77 


Examples of the use of the solutions (20.2, 3, 4) in potential 
theory are given in Coulson’s Electricity (Oliver & Boyd, 
1948) §§ 75-78, 80; a further example is given below. 

The Legendre functions of the second kind, Q,(cos 0), 
which are absent from problems involving spherical 
boundaries, enter into the expressions for potential func- 
tions appropriate to the space between two coaxial cones. 
If0<a<0<f<z we must take a solution of the form (20.1). 
Suppose, for example, that y = 0 on 0 = a, and W = 2a,r" 
on @ = f then we must have 


A,,P,(cos «)+C,Q,(cos a) = Ὁ 
and 
A,P,(Cos f)+C,0,(cos B) = Gy, B, is D,, = 0, 


the latter results following from the fact that if « # β, 
Q,(cos «)P,(cos B)—P,(cos «)Q,(cos Bf) does not vanish. 
Solving these equations for A, and C, and inserting the 
solutions in equation (20.1) we find that in the space 
between the two cones 


ys δ “ae igi a)P,(cos 0)— P,(cos «)Q,(cos am (20.5) 
0,(cos a)P (COS B) - (COS x)Q,(Cos B) 


To illustrate the use of the solution (20.1) and of some 
of the properties of Legendre functions we shall now 
consider the problem in which an insulated conducting 
sphere of radius a is placed with its centre at the origin 
of coordinates in an electric field whose potential is known 
to be 


an=0 


Ἐ ας" Ρ (οο5 θ) (20.6) 
n=1 


and we wish to determine the force on the sphere. The 
conditions to be satisfied by the potential functions wy are 
(i) that w is a solution of Laplace’s equation; (ii) that y 
has the form (20.6) for large values of r; (iii) y = 0 on 
r= a, 
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The conditions (i) and (ii) are satisfied if we take 


v= Σ (x ot," + ΕΞ P,(cos θ) 


and (iii) is satisfied if we write B, = --αι αὐτο, We 
therefore have 


v= Σ α(ν 


qz"tt . 
- a P.(cos 0). 
The surface density of charge on the conductor is 
¢= - (%) ἜΝ σ᾿ δ᾽ (2n+1)a"~ 'a,P,(cos θ) 


or ΝΣ 4nn=1 


and since the force per unit area on the conductor is 2x0? 
the resultant force on the sphere is in the θ᾽ = Ὁ direction, 
and is of magnitude 


F= | 2x07 2na* sin 0 cos 0 dé 
: (20.7) 
= 1g? y .. 2 sgaling 
n=1m 
where /,,,, denotes the integral 
(2n+1)(2m+1) | cos @ sin θ P,,(cos 0)P,(cos 0)dé. 
0 


Changing the variable of integration to μ = cos@ and 
using the recurrence relation (14.2) we find that 


i _ (2m+ 1) αν 1) ᾿ PrP n+ ι(μ)άμ 
+n | Ρ m( LP, n= wan} Ι 
-} 


and by the orthogonality property (15.8) this reduces to 
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the form 
κα Sa (n+ 1)δ,,, nti 26m, nei (20.8) 


Substituting from (20.8) into (20.7) we find that the total 
force on the sphere is 


oo 


F= Σ (nt Vata", 


§ 21. Legendre’s associated functions. We saw in 
example 1 of Chapter I that the solution of Laplace’s 
equation in spherical polar coordinates reduces to the 
solution of the ordinary differential equation, 


dO. 
(1-- wy To --μ 6 εὐ 


which reduces to Legendre’s equation when m τ θ. This 
equation is known as Legendre’s associated equation. To 
solve this equation we may write 


© = (u?—1)7*"y. (21.2) 


Substituting this expression in the differential equation we 
find that » satisfies the equation 


ho 0 (21.1) 


(1-- “BK πὶ a ¥. = (1m) > +(n+m)(n—m-+1)y =0 
μ 


and tie this equation m times with respect to x 
by Leibnitz’s theorem we find that 


\a- ὑπ τρια snarl δ᾽ ae? (21.3) 


showing that if d"y/du™ is a solution of Legendre’s equation 
(17.1) the function ©, defined by equation (21.2), is a 


solution of Legendre’s associated equation (21.1). 
Similarly if we put © = (u?—1)*"y in equation (21.1) 
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we find that 
(- τὴ τι y —2(1-+m)ue ον > +(n—myntm+ ty = 0. (21.4) 


If now we ong equation (17.1) m times with 
respect to 4 we obtain the equation 


0-0) ἐς Ae δος ἐν ΠΑ +(n— —mynem Ὁ ΕἿΣ τεῦ 
| ἐμ ἐμ" 


showing that if (4) is a solution οἵ Legendre’s equation 
then 


(tie Se (21.5) 
dp” 
is a solution of Legendre’s associated equation (21.1). 
Taking the two solutions of Legendre’s equation to be 
P,(u) and Q,(n) it follows from (21.5) that the functions 


PO) = (u2—1)3" Pal) 
ἐμ" 


Oni) = oP SO (21.6) 


are solutions of Legendre’s associated equation. As a 
consequence of equation (21.3) we see that so also are the 
functions 


P>™(y) = μ3--) τ" | ᾿ | Ep | Peds” — (21.7) 
1 1 1 


and 
aos ᾷ 
2, "(μ) = ue | | ve | Q,(¢(d¢)". (21.8) 
It is an immediate generalisation of (6.5) that 


F(a, Bs ¥ 2) = Oats »F\(a+m, B+m; y-+m; x) 
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so that using Murphy’s form (15.2) for P,(j) we see that 
3 _ (— Malt Dn 
ee Py it Ne ieee 
an) = a! 

oF; (mn n+m+1;m+1; a 


Hence P?() as defined in (21.6) may be written in the 

form 

: I(n+m-+1) 9 re 

P ------- ——— (μ΄--Ι 
"Ὁ = πηγας τὴ ἢ 


κοι (m m—n, n+m+1; m+1; 13) (21.9) 


Other expressions for the first of the two functions (21.6) 
can be easily derived. If we make use of the result (7.4) 
we see that 


I(n+m-+1) 


OE df) Maite! Se a- 3m ι--- 1 
2°m'!T(n—m-+1) ee 


P(t) = 


x oF, (m= —n; m+]; "ἢ (21.10) 


and similarly, if we make use of relation (7.6) we may 
derive the expression 


_ Γ(πμ ἐπι 1) tm 
Pat) = m!T(n—m-+1) At) 


X oF, (n+, —h; m+1; 5 i) (21. 11) 


From Rodrigues’ formula (15.3) we derive the simple 
expression 


Paw) = = (yu? 1)" aw! ~1)". (21.12) 
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The simple differentiation 
ge ar _. =(- 1)” [(m+n+1+2r) —m~-a~1l—2r 
εἰμ" T'(n+1+2r) 


may (because of the duplication formula) be written in 
the form 


r(n+ Se = +4),($n+ Md 


at+1+2r 
μ 


= (—1)"T(m+n+1) amir: mma 


showing that, by term-by-term differentiation of the solution 
(17.7) of Legendre’s equation, we obtain the solution 


mE n++1) nant 
2"*1T(n+3) 


X oF; ΤΩΣ 4m+4n+1; π- κα; ‘) (21.13) 
μ 


Οπ(μ) = (-- 1)" ἰ(μῖ -- η)}5 


of Legendre’s associated equation. 
Solutions of the type (21.7) and (21.8) can be derived 
in a similar fashion. Using the result 


᾿ [- | (34) a ΩΝ, ae Tew) 


in equation (21.7), with Murphy’s expression (15.2) for 
P,(u), we derive the expression 


| (μ--1)",}, (-» n+1; m+1; ἘΞ i) (21.14) 
m! - 


μ [ἢ ᾿ 
| | τ | P,(EaE)", 
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and this yields the solution 
=m Ooty el τα : ἘΞ 'Ξ (21.15) 
ΤΡ ΡΤ ΕΝ ᾿ Ε , 3 
The solution provided by Rodrigues’ formula (15.3) can 
obviously be written as 


2 4,73" ,1π- πὶ 
SS ee. (21.16) 


In a similar way the solution 
Γ(ΌΓίη -- αν 1) (u?—1)7*" 
“πὲ 'Γί(η + 3) 7 ΔΙ 1 


Xo F, (sn —4m+4,4n—4m+1; n+3s4) (21.17) 


Pe™(u) = 


05") =(—1)" 


is derived from equation (21.8) and the result 


ron ffi. f Seen ap 


= (—1)"T'(n—m+]) (ἐπ- ἐπι δ), (ξη-- ἐπι +1), — nae +), 
used in equation (17.7). 

The four functions Ραμ), Ομ), Pr (4), Qn "Ὁμ) defined 
equations (21.9), (21.13), (21.15) and (21.17) respectively 
are therefore solutions of Legendre’s associated equation. 
They are known as Legendre’s associated functions. Al- 
though the expressions above have been found by assuming 
m and nto be integers it is readily shown that the solutions 
quoted are valid even when m and ἡ are not integers. 
Since Legendre’s associated equation is of the second 
degree it follows that only two of these four functions are 
linearly independent, and that the other two may be expressed 
simply in terms of them. It follows immediately from 
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equations (21.11) and (21.15) that if m, n are integers 


r(n—m+1) 


Pru) = 
"Ὁ C(n+m-+1) 


Prt). (21.18) 
Furthermore, if we apply the result (7.6) to the hyper- 


geometric series on the right-hand side of equation (21.17) 
we find that 


om m1 (3)T(n—m+1) μ’"πη-τ(μ2-- gyi 
Q, (H) = (—1) 21TH 3) 1)? 


x >F, ΤΩΝ ἐμιἘλη ἘΠ: n+3; 3) 


which, on comparison with equation (21.13), reveals the 

relation 

I'(n—m-+1) 

I(n+m-+1) 
It is now a simple matter to prove that when m and ἡ" 

are integers, and m is fixed, the polynomials ΡΠ μ) form 


an orthogonal sequence for the interval (—1, 1). Making 
use of the results (21.18), (21.12) and (21.16) we find that 


oe r(ntm+1) 1 
P™y)P™()du = OT Met) _ 
᾿ «(μ) Ρι(μ)άμ Som OF ae. 


2, "Ὁ) = Qn(H). (21.19) 


I. <— (μον St -1y"d 
dy" m dy π΄ +m μ 


and after integrating by parts n—m times the expression 
on the right reduces to 


(n+m+1) (- πὶ Ἃ ia -y 4 ἀπ᾿ 


ΔΑ ΠῚ ΠΡ 1) ets 
T(n—m+1) 2"*"'nin'! en (1 ~1)"du. 


This integral is evaluated by the method used at the end 
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οἵ § 15 and we find that 

I(n+m+1 = 

Sars 1 τ 
n+1) 


In many physical problems μ = cos θ so that —I1 Syl. 
It is then not always convenient to have a factor of the 
form (u7—1)*™. We use instead Ferrer’s function 


ἢ PM(u)PRu)dp = 25,» (21.20) 


see Τὰ 
Tou) = 1 — py ΠΕ ΔΟ, (21.21) 
ἐμ 
With this notation we may write (21.20) in the form 
᾿ Γι ἢ) 2 
ΤΑ Τρ μ)άμ- --------- - 6, «. (21.22 
ΒΡ WT) T(n—m+1) 2n4+1 ) 


The other formule are amended similarly. 


§ 22. Integral expression for the associated Legendre 
function. If we assume Cauchy’s theorem in the form ἴ 


: 5 LO. ge = fiw) 
πὶ 


where ὦ is an analytic obi of the complex variable 
¢ in a certain domain R which includes the point ¢ = mu 
and where the integral is taken along a closed contour C 
which includes £ = y and lies wholly within the domain Καὶ, 
then by differentiating both sides of the equation r times 
with respect to 4. we obtain the result 


ἀμ) ἐν" ἘΝ 

d 22.1 
dw “3 ε(ζ- (C—py*? A wide 
Substituting m+n for r and (¢*—1)" for f(0) in this equation 


+ E. G. Phillips, Functions of a Complex Variable (Oliver and Boyd, 
1940), p. 93 
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we find that, as a result of equation (21.1), 
(m+n)! (n*~1)™" [ (ζ 1} 
ἢ τος —-—__—__df, (22.2 

(a ᾿Ξ _ ni “πὶ Ιε(ζ-- ἘΜΈ ξ ( ) 

If «>0 we may take the contour C to be the circle 
[f-n | = | J@?-D. 

Integrating round this contour we obtain from equation 
(22.2) the equation 


2n 
af ia aT ΩΝ 


~ (n + = sin in my)PntH) 


from which follows immediately the Fourier expansion 
{ut+./(u?—1) cos (¢— γε 


= Ρι(μ) 2 me στ Ραμ) cos m(y—). (22.3) 


ae m)! 
Changing ἢ to —(n+-1) we obtain the expansion 
{u' + ψ(μ-- Ὁ cos ψ}-"-! 
= P(p')+2 ᾿ (--1Κ" ae PM) cos my. (22.4) 


Applying Parseval’s theorem for Fourier series to the 
series (22.3) and (22.4) we find that the series 


PAWPAW)+2 Y(t" PO" PrUdPra) cos mg 
converges to the sum 


Ὁ ἰσενία᾽ Doe Wry 5, 
2n Jun {u'+./(u'?-1) cos y}"F 
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This integral may be evaluated by means of Cauchy's 
theorem.} Its value is found to be 


P,{up! +./[(u? —1)(u’? —1)] cos 9}. 
Writing μι = cos 0, μ' = cos θ΄ and 
cos © = cos 0 cos 0’+sin @ sin 0’ cos @ 
we obtain the result 
P,(cos ©) = P,(cos 6)P,,(cos 6’) 


$2 = (- ἢ» ἐστον P"(cos 0)P™(cos0),cos(m@) (22.5) 
which is often of value in the solution of problems in wave 
mechanics. 


§ 23. Surface spherical harmonics. From the two sets 
of orthogonal functions T7(cos 0), cos (md) we can form 
a third set of functions 


Xn, mO; φ) 
 {2πΈ1}} {(n—m)!)* τ ἀΐ | 
af 2n ert ta xwondjocs me, (aren) 053. 


which is an orthogonal set of functions on the unit sphere, 
i.e. the functions of the set satisfy the normalisation relation 


sin 940" ΘΠ 
0 0 


In a similar way we can construct a set 


Y,, ἴθ, Φ) = (3) ἿΞ ata T,"(cos @) sin md, 


(n+m)! 
(msn) (23.3) 
+t For details, see E. W. Hobson, op. cit., pp. 365-71. 
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which satisfies the relations 


ὄπ 2π 
| sin 0d0 | γ᾿, μι ἴω", m4 = Onn Omm's (23.4) 
0 0 


π 2π 
| sin 94} Xn, ἴω, madd τε ῦ (23.5) 
0 0 
for all integral values of ἡ, η΄, m and m’ with msn, m'Sn’. 
Because of these orthogonality relationships we can 
establish an expansion theorem which is a straightforward 
generalisation of the Legendre series (16.3). It is readily 
shown that for a large class of functions f, the function 
ΚΘ, @) can be represented by the series 


y c,P,(cos 0) 


n= 0 
[- ε] 
+ YY χω, κίθ, $+ YanYosn(O, 4} (23.6) 
where the οὐδ Cys Χημν Yam 816 given by the ex- 
pressions 


τς ΣΙ | ag |’ ΚΘ, ¢)P,(cos 6) sin 040, (23.7) 
a ᾿ sin odo |” X,, (0, ΦΥ(θ, ddd, (23.8) 
0 0 


ἡ ᾿ sin 00 Ν Y, κ(θ, d)f(0, φγάφ. (23.9) 
0 0 


For any given function 310, @) the series (23.6) can therefore 
in principle be computed by a series of simple integrations. 

The functions X, ,, and Y, ων which are known as 
surface spherical harmonics, can be constructed easily from 
the known expressions for the associated functions T,”. 
We find, for instance, that 


ὅγε. 
X,, 10, Φ) = -- (2) sin 0 cos ¢, 
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δ". 
X>, (0, φ) = -- ῳ sin 0 cos θ cos φ, 


X,, εἴθ, ¢) = ᾿ sin? 8 cos 2¢, 


105 \" .in? @ cos @ cos 26, 
ὅπ 


2) sin 0(5 cos?  -- 1) cos ¢, 
Xx 3, (0, Φ) = sy 


X10, 6) = ἃ 


Χο, χ(θ, Φ) = - (23) sin? 0 cos 34, 


and the corresponding expressions for the Y, , are ob- 
tained by replacing cos (m@) by sin (m@) in the expressions 


for the Χ΄,, mn: 
The functions X,.,, and Y,,, have the important 
property that they are solutions of the partial differential 


equation 
to (sm) PX n+ 1X πὸ (23.10) 
sin 0 60 60 sin?0 ag? 
so that the function 
(Ar"+Br-"~")Xq, m(O, 6)+(Cr"+Dr-"~*)¥y, κ(θ, $)s 


where A, B, C and D are constants, is a solution of Laplace’s 
equation. It follows immediately from equations (23.7)- 
(23.9) taken with the expansion (23. 9 that the function 


Wr, 0, Φ) = pa Ch (J P,(cos 0) 


+ Σ y 1 (J {XamXn, κ(θ, 6) + Yam Yn, m8, $)} 


ΝΞ 1 πὶ ΞΞ 
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satisfies Laplace’s equation in the region 0<rsSa, is finite 
at r = 0, and takes the value 710, Φ) on the sphere r = a. 

For example, suppose we wish to find the solution of 
Laplace’s equation which takes on the value x”? on the 
surface of the sphere r = a. Here we have 


t(0, ¢) =a? sin* 0 cos” @ 
a* a? ( cos? @—1 


ξαὶ τρῶν τὰ 5 ) 4a? sin? θ cos 2¢ 


te 


= — — -- ΡΣ | —]} X, (0, 2 
3 3 2(cos 6)+ ( 3 2, χίθ, φ)α 


Thus the required solution is 
1 1 4 \? 9 
ψίν, θ, d)= 3 a? — 3 r’P2(cos 0) + (4) r’X>, (0, 9). 


Substituting the values of P, and X, 2) and transforming 
back to Cartesian coordinates we see that the required 
solution is 

Ww = 4(a?+2x?—y?—z?), 


§ 24. Use of associated Legendre functions in wave 
mechanics. To illustrate the use of associated Legendre 
functions in wave mechanics, we shall consider one of 
the simplest problems in that subject—that of solving 
Schrédinger’s equation 


2 

V+ a (W-V)y =0 (24.1) 
for the rotator with free axis, that is for a particle moving 
on the surface of asphere. In equation (24.1), W represents 
the total energy of the system, V the potential energy. 
In the case under consideration V is a constant, Vg say, 
and the wave function y will be a function of 0, ¢ only. 
If the radius of the sphere is denoted by a then equation 
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(24.1) is of the form 
1 ὃψ , coté dw 1 @ | 8x?m(W-V,). 
=~— + + — + —_—_.—— yy = 0. 
a* 607, a* ὃθ. a* sin 0 dd” h? v 
(24.2) 

If we consider solutions of the form 

ψ ἐς @e=in? 
then 

age FE 2 

@"+0' cot 04+ SEMA W 9. πὸ 6... 0. 
h sin* 0 

Substituting 


82*ma?(W—Vo) 
a oe 
we find that this equation reduces to Legendre’s associated 
equation (21.1) and hence has solution 


© = AP"(cos 0)+BO"(cos 0). 


However, for the same reason as in the case of potential 
theory (ὃ 20 above), we must take B=0. The solutions 
of equation (24.2) will therefore be made up of com- 
binations of solutions of the form 


Wm, nO, φ) onl Aah ς “ (cos θ), (24.4) 

where A,,,, is a constant. 
The physical conditions imposed on the wave function 
y are that it should be single-valued and continuous. 
Obviously then the “ physical” solutions will have m an 
integer, since W,, ,(9, @+2x) must equal y,, εἴθ, ¢). 
Further, in order that the series for Py'() should converge 
for the values μ = +1 it is necessary that it should have 
only a finite number of terms. This is possible only if n 
is a positive integer. If therefore the solution (24.4) is to 


p= cos 0, = n(n+1) (24.3) 
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be valid for 0 = 0 and 0 = x we must have n a positive 
integer. The physical conditions on the wave function 
are therefore not satisfied by systems with an arbitrary 
value for the energy W but only by systems for which 

W=Vo+ — τος _ n(n-+1), (24.5) 
where n is a positive integer. μῇ other words, the energy 
of such a mechanical system does not vary continuously, 
but is capable of assuming values taken from the discrete 
set (24.5). 


Examples ΠῚ 
(1) Show that, if m is odd P,(0) = 0, and that, if n is 
even, 
P.(0) = (- (— 1)" sn * (= 1)*"n! 
(jn)! 2"{(4n)!}? 
(2) Prove that 


[] μ" n+1 l+yu 
Σ, “aes ‘i Ht 
(3) If € = w+./(u*—1) show that 
ὦ (-hO-*(1—h/-* = ys h"P,(1); 
T(n+4) μη nA ae 
(ii) Ρι(μὴ) = TG) C" F\(4, —n; 4—n; ζ ἢ. 
Deduce that 


(iii) 2), —n; 3 τη; "ἘΠ τς 


(4) If n is a positive integer prove that 


; 2h" 
| Ρι(μχ! —2puh+h?)-*#dp = ——, 
-1 2n+ 
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and hence, making use of Rodrigues’ formula, deduce that 


, intl 2 
| i= py —2uh ἘΞ} td = - (n )* 
-1 


(2n ἘΠ" 
(5) Prove that 
Pix) = Σ 2n—4r—1)P, 2-102), 
r=0 
where p = 4(n—1) or 4n—1 according as n is odd or even. 


Deduce that for all x in the closed interval (—1, 1) and 
for all positive integers m, the values of the functions 


| P(x) |, n7?| Pid], n74| PC) |. 
can never exceed unity. 
(6) Prove that 


1 
[ Ρ(μ)άμ rs στ fe (H)— Pra (0}; 


and deduce, from example 1, that if m is an odd integer 


_ (==)! 
τ τ 


What is the value of the integral when » is even? 


(7) Using equation (14.2) and the results of the last 
example show that, if is even, 
(n—2)! 


i 
P(u)dp = (—1)**-*? —_—__——___., 
[ BP leiden = (—1) 2"(4n+ 1)'(4n—-1)! 
and that the integral has the value zero if n is odd, 
(8) If 
d 
t= | PaadP st) 
-1 μ 


prove that (1+ 1)u,.,+mu, = 2, Hence evaluate ν᾿. 
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(9) If m and n are positive integers, prove that 
nae 2mtn+1 f(y 4 n)1}? 
1+p)"*"P,(u)dp = ———_——-.. 
ix H) (αι m!Qn+2n+1)! 
(10) If m is even and m>-—1, prove that 


‘am _ __ 1Gm+prGmt+)) 
|, we AW aE = 27(4m+4n4+3)Gm—4tn41) 
Deduce that 


1 
ea 1 nee 
(1—kp?)~"~*P,,(w)du = —— κΚὮ -- ) 53. 
x μ΄) an( {aye ind ( ) 


(11) Show that 


+l μ--1} 
= | —— } ,F, | —n; —n, 1; — 
Ρι(μ) (Ξ Ὶ 2 (-»:: nt) 


and hence that 


(1-1"P, (i!) - = ἜΣ (Cyr. 
Deduce that 
P(cosh u)=1. 


(12) If μὴ = (u?—1)" show, by using Rolle’s theorems 
that f ‘(u) must have at least one zero between —1 and 1. 
Proceeding in this way deduce that f() has n zero, 
between —1 and 1. 

Hence show that when n is even the zeros of P,() 
occur in pairs, equal in magnitude but opposite in sign, 
and that when n is odd, x = O is a zero and the others 
occur in equal and opposite pairs. 


(13) If γίμ) is any solution of the linear differential 
equation 


a ΤΣ 2 + Blu) — dy ih a da 0 


LEGENDRE FUNCTIONS 95 


in which «, βὶ and y are continuous functions of 4 whose 
derivatives of all orders are continuous, prove that γ(μ) 
cannot have any repeated zeros except ome for values 
of μ᾿ which satisfy the equation α(μ) = 

Deduce that all the zeros of Ρ, (μ) ᾿ς - distinct. 


(14) Prove that the Legendre polynomial P,(x) has the 
smallest distance in the mean from zero of all polynomials 
of degree n with leading coefficient 2"(4),,/n! 


(15) If R denotes the operator 
d d 
aE 1- 3 ea a | 
dx \ ὼ a 
prove that 


| PX)R{S()}dx = —n(n-+1) " P(x) f(x)ax 
Jj—] -1 


provided that f(x) and f ‘(x) are finite atx = +1. 
Prove that if n=1 


᾿ . 2 
‘a log (1—x)P,(x)dx = "ant Ἢ 


(16) Prove that 
: 7 
(i) P,(x)dx _ 2,/2 . 
"ὦ saan 2n+1 
(ii) P.(x)dx E 2h" . 
" i. —2hx+h2)t  1—h? 
a7) If R? = 1—2hx+h’, prove that if | |<1,n21, 


ie (h—-x+R πε." ON 
(i) ¥ oe | aa } pyoods = aT 


2h" 
n(2n+1) 


(ii) | log (1—hx+R)P,(x)dx = — 


seston MR eet 
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(18) If 
gener 
᾿ ------- dé, (m>l), 
flx,a) = [ are (mm 
prove that 
1 qt 
| I(x, a)P,(x)dx SS ee ee 
5 (n+m)(2n+ 1) 
(19) If z is real and | z |<1 prove that 


| ee. Oe. 

io 1+zcos@ ./(i—z?) 

Putting z = Fh,/(u*—1)/(1—Ap) where fA is so small that 
| A{ut/(e?- 1) cos @} |< l 


(O<¢<nx), expanding both sides in powers of ἡ, and 
equating coefficients of h" show that 


Oy ie i 
Ps) = * | fut. /(u2—1) cos 4}"d¢. 
0 
Hence evaluate the sum 
Ἢ "ΟΡ (cos 0). 
r=0 


(20) Show by making the substitution 
z= th (μ᾽ -- Πάμ-- 1) 
in the formula in ex. 19 that 
d 
ΤῊΝ rf cee gyre 
(21) Making use of the integral expression for Ρ, (μ) 
derived in ex. 19 show that 


Ρι(μ) = Σ. "C.1— ὙΠ ΤΡ, (0). 


7 


4 
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Deduce that 


Pou) Pile) Ρχ(μ) 
Ρι(ω Pn) P;3(p) 
P(t) P3(u) P,(u) 
(22) Prove that if p> 1 


Ο,(ὴ = τεῦ 
and deduce that 
0,(u) = ἢ (μ-- ψιμ-- Ὁ cosh 6" 40, 
where « = 4 log (u+1D)/(u— 
Hence find expressions a Oo(u) and 0 ,(μ). 


(23) Determine the simple expressions for Qp(w), 
O, (4), O24) | and Ὁ 3(μ) by working out the values of the 
polynomial W,,_ ,(u), occurring in equation (18.8), for these 
values of ἡ. 


(24) Establish the following formule due to MacRobert: 
ὦ μρ, μοι = αὶ [ PHO ae, men: 


P0) 0 P,(0) 
0 P,0) 0 
PO) O P40) 


= (1—,")? 


(i=1F 4, 
= (u—t)’*! 3 


PAO! ge 4 
1 μ['ἷς ryt 


(ii) "(u2—1)04(n) = 4 | meee dé, 
-.1 


(δ pP,(y)O,(1i) = a 


n>m; 


iv) αὐ —DPA(WOAW) = 4 [ EVP ge, >m 
MG τ “ety that, if m is a positive integer and ¢ = n+ 
2 — 


ὉΠ ΓᾺ * P (cos 8) sin θ dé 
Qu(H) 4) 1-- 2ζ cos 0407 
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Deduce that, if | {|>1, 


Ο,(μ) = Σ ΚΕΝ | P,({cos 0) sin μιθ dé. 


By evaluating this integral show that 


in) = weg t F,(d, n+-1; π͵αξ; C72) 


(26) Prove that, if m is a positive integer, 


. ἢ -- mam (—1)"(2m) (1 — 7)?" 

ΠΗ Τα (p) = ————————.. 

εἶς nH) 2m \(1—2yh+h?)"** 

(27) Show that, if | « |>1,n>—1 then 
mc) — L(a-tm-+1) (u?—1)* q (1—€)"dé 
OF) T(n+1) gatl We (u—Eytmtt 
Deduce that, if | n+1|>2 then 

ΓΙ = Jat (n+m+i1) (μῆ -- 1)}" 
C= Ta) Gal 


oF; (n+2, n+m+1; 2n+2; a) 
+l 


Find a simple expression for Q"**(y). 


(28) Prove the following recurrence relations for 
Ferrer’s Associated ΠΝ Functions: 


(Ὁ THY CG) — T' () = (2n+ DA -- pw? AT); 


(ii) (n—m+ yr n+ 1H) —(2n-+ LT "(u) 
+(n+m)T,~ (4) = 0; 


(ii) ΤῊ ἢ) -- Tr) = (n—m)(1— 1? PT). 
(29) Derive the expressions for T;'(0, 6) and X (0, ¢) 
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form = 1, 2, 3, 4. Express the functions sin? @ sin’ φ. 
sin? @ cos* x sin @ cos θ cos ¢ in terms of surface spherical 
harmonics. 


(30) Find the function which satisfies Laplace’ 8 equation 
in the interior of the sphere x?+y?+z7 = a’ ; remains 
finite at the origin and takes the value «x?+fy?+yz? on 
the surface of the sphere. 


(31) TheJacobi polynomials are defined by the equation + 
F ἴα, b, xX) = 2F,(—m, a+m; δὲ x). 
If Ξε d/dx show that 
D™{x°*™— 1, — xt") =(b),,.x° | Fy (b-—a—m, b+m; b; x) 
Deduce that 


(i) F(a, by x) = aa ΩΝ aye}; 


(ii) ᾿ χ᾽ '(. --χ) f(x) F (α, b, x)dx 
0 


- !}. (-- ΟΣ: 11 -- etter ax: 


(iii) |, x? (1 —x)* °F,,(a, b, x) F,(a, δ, x)\dx =0, m+n; 
Jo 


1 
(iv) | χ᾽ 11 —x)*""{F¥,,(a, b, x)}? 
0 
7 I(b)l(a+1—b)(a+1—b),, “ m! 
T(a)(a),.(B) mn a+2m 
(32) The Tchebichef polynomials of the first and second 


+ It should be observed that this name is sometimes applied to 
the polynomial 


p® B) (x) = (a τ as 


F (α ἘΡ-1, o+1, 4-23). 
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kinds may be defined by the equations 


T(x) = 2F; (-» n; 4; =) 


U,(x) = (n+1),F, (- n+1; 3; 3) 


respectively. 
From (2) on p. 43 show that 
T,(cos 6) = cos ηθ, U,{(cos 6) = S2@+H8 
sin 0 


Prove also that 
() 142 Σ Tor" = (1-2-2 $2)-4, 
n= 1 
|t|<1, |x]<t1; 
(ii) 2"(4), T(x) = (-— Ὁ — xn ἐς = (=x) εἰ" 
(iii) [ T,,(X)T,(x)(1 — x7)" #dx = 4205 yn 
=4 

(iv) ¥ U(x)" = (1-2xt 2)"; 

n=0 ΄ 
(v) 2"**(4),4,U,(x) 


= (-- 70. -Ἐ1Χ1-- ἘΚ τ α- =x yrs. 


(vi) [ U,,(x)U,(x)(1 —x?)*#dx = ἀπὸ, ει: 
-1 


(vit) T(x) =U,(x)—xU,_ 40); 
(viii) (1—x*)U, _ s(x) = xT,(x)— T+ χοῦ. 
(33) The function C}(x) defined by 


(1—2xt+12)-" = a σὺ (v>0) 


is called the Gegenbauer polynomial of degree n and order y. 
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Prove the relations: 


(Ca) = EF (n+2s —n; v+4; δΞ =") 
wrt) — (— 16% 1) — Let 2). 
(ii) (: (1) = (— 1)"e,(— 1) nIFQn) ° 


(iii) Cx) = (-2)" 
Cet ΓΟΡῈ πὴ ati Zaye 
nil (vy) (2v+2n) dx" ᾿ 

γὴν Ι CH(x)C2(x)(1 x2)" τάχ = less 

-1 


πἰν ἘΠ) ΓΟ] δ᾿ 


CHAPTER IV 


BESSEL FUNCTIONS 


§ 25. The origin of Bessel functions. Bessel functions 
were first introduced by Bessel, in 1824, in the discussion 
of a problem in dynamical astronomy, which may be 
described as follows. If P is a planet moving in an ellipse 
whose focus S is the sun and whose centre and major axis 
are C and A’A respectively (cf. Fig. 8), then the angle 


A' 


Fig. 8 


ASP is called the true anomaly of the planet. It is found 
that, in astronomical calculations, the true anomaly is 
not a very convenient angle with which to deal. Instead we 
use the mean anomaly ζ, which is defined to be 27 times 
the ratio of the area of the elliptic sector ASP to the 
area of the ellipse. Another angle of significance is the 
eccentric anomaly, u, of the planet defined to be the angle 
ACQ where Ὁ is the point in which the ordinate through 
P meets the auxiliary circle of the ellipse. 
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It is readily shown { by a simple geometrical argument 
that, if e is the eccentricity of the ellipse, the relation 
between the mean anomaly and the eccentric anomaly is 


€=u-—esinu. (25.1) 


The problem set by Bessel was that of expressing the 
difference between the mean and eccentric anomalies, 
u—(, as a series of sines of multiples of the mean anomaly, 
ie. that of determining the coefficients cr = 1, 2, 3, ...) 
such that 


u-€= δ᾽ c¢,sin (rf). (25.2) 
r=1 
To obtain the values of the coefficients c, we multiply 


both sides of equation (25.2) by sin (sf) and integrate 
with respect to € from 0 to π. We then obtain 


i (u-O)sin(sal= ¥ ς, | " gin (rf) sin (s0)d€. 
0 : να, 0 
Now 
| sin (rf) sin (s¢)d{ = 426, , 
0 
and an integration by parts shows that 


[ (u—D sin (σζως 


= ak (€—u) cos eo + al (F - ἢ cos (5ζ)ζ. 


From (25.2), ζ-τ τ is zero when € = 0 and when ¢ = π, 
so that the square bracket vanishes; the integral can be 


written in the form 
al cos sf du 


D. E. Rutherford, Classical Mechanics (Oliver and Boyd, 
951. § 42. 
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and hence, using equation (25.1) we obtain the result 
Ὁ, ΞΞ Ξ | cos {s(u—e sin u)}du. (25.3) 
TS Jo 


The integral on the right-hand side of equation (25.3) 15 
a function of s and of the eccentricity e of the planet's 
orbit. If we write 


πῶ" | cos (x sin 6—n6)d0 (25.4) 
Tt Jo 


it follows from equations (25.3) and (25.2) that 
ee 5 ley ἫΝ (25.5) 


The function J,(x) so defined is called Bessel’s coefficient 
of order n. 

We shall now show that J,(x) is equal to the coefficient 
of rf" in the expansion of exp {4x(t—t~*)}; in other words 
we may define J,(x) by means of the expansion 


exp {as (-- x = : J A(x)". (25.6) 


= 


To prove this we need only show that the J,(x) of (25.6) 
can be expressed in the form (25.4). We first of all observe 
that 


Σ (vrs Σ 140)(- +) = ἜΣ J 


AEP ὦ 


since both expansions are equal to exp {4x(f—1/p)}. 
Equating coefficients of f” we have 


(—1)"J_,(x) = J,(x). (25.7) 


In the expansion (25.6) we may write 1 = e’® to obtain 
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the relation 


exp(ixsind)= ΟΣ J(xe™, 


Making use of the result (25.7) we see that the series on 
the right can be put into the form 


ΠΟ Σ Tan) 208 26) +2: Σ Jams (x) sin (2m+1)0 


so that by equating real and imaginary parts we obtain 
the expansions 


cos (x sin θ) =Jo(x)+2 δ᾽ Jom(x) cos (2m6), (25.8) 
m= 1 
sin (x sin )=2 δ᾽ Jomei(x)sin(2m+1)9. (259) 
m= 0 


If we now multiply (25.8) by cos nO, (25.9) by sin nd, 
integrate with respect to 0 from 0 to z, and use the formule 


ἣ cos (m6) cos (η0}) ἀθ = Ε sin (m0) sin (n0) d0 = = 5, 
0 0 | 
we obtain the formule 
J x) = af cos (x sin 0) cos (n@) ἀθ, (n even), (25.10) 
0 


J Ax) = Ξ al , sin (x sin 0) sin(n0)d@, (nodd). (25.11) 


Because of the periodic properties of the trigonometric 
functions we know that the integral on the right of equation 
(25.10) is zero if n is odd, while that on the right of equation 
(25.11) is zero if mis even. Thus for all integral values of ἡ, 
we have 


J A(x) = : ᾿ {cos (x sin 0) cos (10) +sin (x sin @) sin (n0)}d0 
0 


which Lis identical with the expression (25.4). 
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In particular 
Sfx) = 2 | cos (x sin 6) dé. (25.12) 
™ Jo 
In what follows we shall assume that the Bessel co- 


efficients of the first kind are defined by equation (25.6) 
or, which is equivalent, by equation (25.4). 


§ 26. Recurrence relations for the Bessel coefficients. 
If we differentiate the generating equation (25.6) with 
respect to x we obtain the relation 


(of =, τον 


which is equivalent to 
5S Git JG) τ ΣΟ Saye" =0. 


Equating to zero the coefficient of ¢” we obtain the relation 
25x) = ,. 18) —Ins 1(2)- (26.1) 


On the other hand, if we differentiate (25.6) with respect 
to ¢ the resulting equation is 


4x ( + 2) exp 1 ( - 7 -- 5 " nJ,(x)t"~* 


and this is equivalent to the relation 


ΟΣ ἀπ ρον — Σ nd,(x" =O. 


ἘΞ - ὦ em") 


Equating the coefficient of {ΠῚ to zero we obtain the 
recurrence relation 


= TX) = Jn— 10) + Sng 1)- (26.2) 
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Adding equations (26.1) and (26.2) we find that 

XJ (Χ) = χα, (x) —nJ, (x) (26.3) 
and subtracting equation (26.1) from (26.2) we obtain 

XJ ,(X) = nJ (x) — XT εἰ). (26.4) 
Putting ἡ = 0 in this last equation we have the important 
special case 

o(x) = —J,(x), (26.5) 


and putting m = 1 in equation (26.3) we find that 
; 1 
JX) = Jo(x)— es 1(*). (26.6) 


Differentiating both sides of (26.5) with respect to x and 
making use of the result (26.6) we have 


: i 
Jo(x) = —Jo(x)+ = J s(x) 
which, as a consequence of equation (26.6), may be written 
Je(x)+ = Jefe) +Jo(x) = 0 (26.7) 


showing that y = Jo(x) is a solution of the differential 
equation 


ἐν 1 dy 
--Ξ + -— = 0, 
dx? xdx bi pk 
We can show similarly that the Bessel function J,(x) 
satisfies the differential equation 


d*y idy Ἢ 
oie ae eae [1..ὄ τ νκκῇ i . (26. 
ἘΠ thon + ( Ὥ y (n integral). (26.9) 


For, from equation (26.4) we find, as a result of differ- 
entiating both sides with respect to x, that 


XJ (x) + I(x) = αὐ, (χ) — Ing 1%) — XS 1(X). 
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Now, from equation (26.4), 


n? 
nJ (x) Ἂ- x? J (x) — 1S y+ 1(x), 
and putting n+1 in place of m in equation (26.3) we see 


that 
χῦ, 1%) + (n+ LIne 1%) = XIX) 
so that 


Silo) +i) = Iya) 250) 


which shows that J,(x) is a solution of equation (26.9) 
provided, of course, that ἡ is an integer. 

As we pointed out in § 1, equation (26.9) is known as 
Bessel’s equation. What we have shown is that if the n 
which occurs in Bessel’s equation is an integer, one solution 
of the equation is J,(x). It is because of this fact that 
Bessel coefficients are of such importance in mathematical 
physics, for as we saw in ὃ 1, the equation (26.9) arises 
naturally in boundary value problems in mathematical 
physics. 


§ 27. Series expansion for the Bessel coefficients. We 
shall now find the power series expansion for the Bessel 
coefficient J,(x). If we write 


1\{ _ whl —— 
exp » (-- \ = exp (4x1) exp ( x) 


and make use of the power series for the exponential 
function we obtain the expansion 


πο 2.965 68 


r=0 27] s=o 2 Ἐς 
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By our definition (25.6), the Bessel coefficient J,(x) is the 
coefficient of ¢" in this expansion. If n is zero or a positive 
integer, we find that 
, oo ai 1) (x at+2s as 
J (x) = Sind) ὁ - 27.2 

) 2m ἐν (3) la 
and when n is a negative integer we can deduce the series 
for J,(x) from equation (25.7). 

Writing equation (27.2) in the form 


5 


J a -_ ἢ ἴἤἴ.. 1.2. 
a) 2"n! Σ, 5101. Ἔ η οἷ: 
we see that 
J,(x) = = oF y(n +1; —4x?). (27.3) 


The variation of the Bessel coefficients Jo(x), J,(x), 
J,(x) for 0SxS20 is shown graphically in fig. 9. These 
are the Bessel coefficients which occur most frequently in 
physical problems and their behaviour is similar to that 
of the general coefficient J,,(x). 

Simple relations for the Bessel coefficients may be 
derived easily from the series expansion (27.2). For 
example, since this equation is equivalent to 


d- oO (—1)*x?"*25 A aia 
EAE Ὁ si(n-+s)! (5) ac 


it follows, as a result of differentiating both sides of this 
equation with respect to x, and making use of the fact 
that (2n+-2s)/(n+s)! = 2/(n—1+45)!, that 


d a (=1)9 72" 1 n-1+2s 
ee NE τ ON Ne te 
dx IO =D, si(n—1+5)! (5) 


which, by comparison with (27.4), shows that 


d 
dx {x"J,(x)} = x"J,—1(X). (27.5) 
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If we write this result in the form 


+ © (F,(x)} = x" 108) 
x dx 


~0-5 | =e 
Fig. 9. Variation of Jy(x), /;(x) and J,(x) with x. 
we see that if m is a positive integer less than n, then 


( ZY XI (x) = χη... Οὐ. (27.6) 
x dx 
Similarly we can establish that 
4 {5,0 }}Ξ — xn 1) (27.7) 
dx 


or, which is the same thing, 
( 4) {x~"J,(x)} = --αῖ "Jai (X), 
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a result which may be generalised to the form 


& A =. -A-m 
(: 4) {x~ "JF u(x} = (— 1) St aX). 
x dx 
In particular we have the relation 


a) Jol), (27.8) 
dx 


which shows how the Bessel coefficients J,(x) may be 
derived from Jo(x). 

Another interesting property of the Bessel coefficients 
also follows from the power series expansion (27.3). This 
concerns their behaviour for small values of the argument x. 
Since 


σπου 


lim of (n+1, —4,?) =] 
ποτῷ 
it follows from equation (27.3) that 


Mee 1 
lim x~"J,(x) = ——. : 
xO 0) 2"n! ad ) 
In other words, for small values of x, the Bessel coefficient 
J (x) behaves like x"/2"n!. 


§ 28. Integral expressions for the Bessel coefficients. 
We have already derived one integral expression for the 
Bessel coefficient of order m (equation (25.4) above). In 


_this section we shall consider other simple integral ex- 


pressions for these coefficients. 
We shall consider the integral 


of | 
I= | (1—17)"" teat 
=1 


in which n>—4. If we develop exp (ixf) in ascending 
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powers of ixt we see that the value of this integral is 


y al (1-12) 8dr. 
ι 


s=o 5! 


If s is an odd integer then the corresponding integral 


occurring in this series is zero, and if s is an even integer, 
2r say, then the integral has the value 


[ σεν ,.. HED D 
0 | I(n+r+1) 
so that 
pe S (HDs Pent ΓΟ ἢ 
r=o (2r)! T(n+r+1) 
| = (—1)'x?" 
= T()r(n+4) 2, rit(n+r+ 12” 


since, by the duplication formula for the gamma function, 
T(4)(2r)! = 22" Γ(γ -Ἐ 3). It follows immediately from the 
series expansion (27.2) for J,(x) that 


- (4x)" ‘ __ 42\n- 4, iat, 1) 
Le)= ey | doer ted 8 


and it is easily shown that this is equivalent to the formula 


J (x) = iran) [ (1—1*)""* cos (xt)dt. (28.2) 
In particular 
2 [* cos (xt) | 
J (x) = τ]. JQ—") di. (28.3) 


The result (28.2) may be expressed in a slightly different 
form by means of a simple change of variable. It we put 
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t = cos θ we obtain the integral expression 


J,(x) = cos(xcos@)sin*"0d0, (28.4) 


x” r 
2. "THT (n+4) Jo 
while if we make the substitution 1 = sin@ we get the 
formula 
ἣν 


5}. ——_+__ 
= ποτ ΓΟ ΤΡ 
The particular forms appropriate to n = 0 are 


| cos(xsin@)cos?"0d@. (28.5) 
0 


ἐπ > [ἐπ 
J o(x) = 2 | cos (x cos θ) εθ = Ξ | cos(xsin@)d0. (28.6) 
T Jo π 


0 


§ 29. The addition formula for the Bessel coefficients. 
In certain physical problems we have to reduce a Bessel 
coefficient of type J,(x+y) to a form more amenable to 
computation. We shall now derive an addition formula 
which is of great use in these circumstances. From the 
definition (25.6) we have the expansion 


exp ΩΣ y) (:- +) = 5 JAx+ yy". 
Writing the left-hand side as a product 


οἰρί ἢ fo ἢ 
and inserting the appropriate series from (25.6) we find 


that 
Σ Sry = ΣΟ Σ IGM. 


aS “- οὐ 


Equating coefficients of 7" we obtain the addition formula 


Jty= Σ FCW). ΟῸ(9.) 
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To put this in a form which involves only Bessel co- 
efficients of positive order we write the right-hand side in 
the form 


πε n a0 

Σ ἢ ΠΟ ,..)Ὲ 2, J Ax)In—- AY) + Σ, In) 
r= = r= r=n 
and note that because of the relation (25.7) the first term 
can be written as 

=i a 
ΟΣ DI = Y (WLC rv). 

Similarly the third term is equal to 


Yo) = Σ τ νον) 
so that finally we have 


I(x+y) = oa σα ,..0) 
Ἔ Σ (— 1)"{J, (x)Jnaly ) Ἕ Jn + Ax)J (ὦ) }. (29.2) 


ΡΞ] 


§ 30. Bessel’s differential equation. We showed pre- 
viously (§ 26 above) that, if m is an integer, J,(x) is a 
solution of Bessel’s equation (26.9), We shall now examine 
the solutions of that equation when the parameter 7 is 
not necessarily an integer. To emphasise that this para- 
meter is, in general, non-integral, we shall replace it by 
the symbol ν, so that we now consider the solutions of 
the second order linear differential equation 


d*y 1 dy. νξ 

--Ξ + = 1— —)y=0. 30.1 

τῶν +( Ὥ m ae) 
Writing the equation in the form 

vty ; +x2 +(—v?+x?)y =0, (30.2) 
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we see that the point x = 0 is a regular singular point and 
that in the notation of ὃ 3, pp = 1 and 40 = —v*. The 
indicial equation (cf. (3.5) above) is therefore 


0*—y? = ἢ 


and this has roots ἡ = +yv. 
First of all we shall suppose that ν is neither zero nor 
an integer. Then the first solution is of the form 


yu > ex". (30.3) 
r=0 


Substituting this series in equation (30.2) we see that the 
coefficients c, must be such that 


oo 
Σ {vtr)\vtr—D+4n—vex7t"+ LY ox ΤΥ ἘΣ = 0. 
r=0 r=0 
Hence we must have 
c,{(v+1)?—v?} =0 
and in general 
c,{(r+v)?—v7} = —c,_., (r=2,3,...). (30.4) 


We must therefore take c, to be zero and hence, in order 
that (30.4) may be satisfied for all γ 2, we must take 


Cor44 = 0 
and 
(- 1)"co 
(2r+2v)(2r+2v—2)...(2v+2)2r(2r—2)...2’ 


an expression which may be put in the form 


= Co (- i) 
r\(v+1), 47" 


Taking cy = 1/2"v! we see that the basic solution of type 


C2, i 
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(30.3) may be taken as 


Χ τ Sant. τ ἢ 4x7)’ (30 5) 
2°T(v+1) r=0 riv+1), 


Comparing this series with the series (27.2) we see that it 
is of precisely the same form as that equation, the only 
difference being that n is replaced here by v. If we take 
the series (27.2) to define the Bessel function of the first 
kind of order mn, even when ἢ is not an integer, then we 
may write the solution (30.5) in the form 


7- J v(x). 
Similarly, if we substitute a series of type 


y= 


y= 2. as 


r=0 


to correspond to the second root of the indicial equation, 
we find that it must be of the type 


: ἀῶ = ey —4x*) | 
Y= arya) reo Fave, ta 


and with the extension of the definition (27.2) to non- 
integral values of v we may write this solution in the form 


y = J_ (>). 
Thus when ν is not an integer we may write the general 
solution of equation (30.1) in the form 


y = AJ (x) + BJ_ (ΑἹ (30.7) 
where J ,(x) is defined if the equation 
I(x) = Ξε 1 WD mimeo te 1s —15), (30.8) 


It should be observed that the results of δὲ 27, 28, with 
the exception of (27.8), are true when n is not an integer, 
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since they were derived directly from the definition (27.2), 
which is equivalent to (30.8). The transition is effected 
merely by replacing n! by Γίν- 1). 

When ν is zero or an integer we know from equation 
(25.7) that the solutions J,(x) and J_,(x) are not linearly 
independent. We must therefore use the formule (3.8) to 
calculate the second solution. 

We shall consider first the case in which ν = 0. If we let 

ο 


Ἐπ ὦ, — 


then in order to satisfy the recurrence relation (30.4) we 
must have 
waxt y (ce) (30.9) 
r=or '(o + 1), 
and putting g = Ὁ we obtain the first solution 


Wo = Jo(x). (30.10) 


Using the result 


a eae = a 1 
de (01), (01), (s=1 ρῈ5 
we see that 
dw 2 (—dx) ae | 
— =w log x—x* ——?. 
de 6 “ὅς r'(o+1), ΡῈ 05 


Putting @ = 0 and substituting the value (30.10) for wy 
we find that the second solution (dw/dg), = o is 


Yo(x) = Jo(x)logx— Σ le’ a or), (30.11) 


where 
¢(r) = 2, =, (30.12) 


The function Y,(x) so povine is called Neumann’s 
Bessel function of the second kind of zero order. Obviously 
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if we add to Yo(x) a function which is a constant multiple 
of Jo(x) the resulting function is also a solution of the 
differential equation 


-—-+y=0. (30.13) 


In particular the function 
lb 
Yo(x) = 3 {Yo(x) —(log 2—y)Jo(x)}, 


where y denotes Euler’s constant, will be a second solution 
of the equation. Substituting from equation (30.11) for 
Yo(x) in this equation we obtain the expression 


| mt 2 τυ = avr 
Yo(x) = 5 {log(x)+7}Wolx)— = Σ (Ξ ΕΣ gor) (30.14) 
π tr=1 (γἢ 
where φί(γ) is defined by equation (30.12). 

The function Y,(x) so defined is known as Weber’s 
Bessel function of the second kind of zero order. 

Thus the complete solution of the equation (30.13) is 
where A, B are arbitrary constants and J,(x), Y,(x) are 
given by equations (30.8) and (30.14) respectively. 

It can be shown by an exactly similar process that when 
v is an integer the complete solution of the equation 
(30.1) is 


y = AJ,(x)+ BY,(x) (30.16) 
where A, B are arbitrary constants, J,(x) is defined by 
equation (30.8) and Y,(x) is given by 


v-1 δ. ». 1 ik mathe 
μος ἔσειφλου- 1 $8 Ot! (2) 


-1 5 CVG)" (46+) +40}. (80.17) 


πεξο τίν)! 


§ 30 BESSEL FUNCTIONS 119 


The function Y,(x) so defined ¢ reduces to the Y,(x) of 
equation (30.14) as v0 and is known as Weber’s Bessel 
function of the second kind of order v. 

The variation of Y,(x) and Y,(x) for a range of values 
of x is shown graphically in fig. 10. 

The functions J,(x) and Y,(x) are independent solutions 
of the equation (30.1), but in certain circumstances it is 


Fig. 10. Variation of Y,(x) and Y,(x) with x. 


advantageous to define, in terms of them, two new in- 
dependent solutions. If we write 
Hx) = J(x)+iY,(x), (30.18) 
ΗΟ) = J(x)-i¥,), (30.19) 


then it is obvious that we can take the general solution 
of Bessel’s differential equation (30.1) to be 


y = A, Hx) + A, HOD), (30.20) 


where A, and A, are arbitrary constants. The functions 
H (x), Hx) defined by equations (30.18) and (30.19) 
+ This function is denoted as N(x) by Courant and Hilbert. 
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are called Hankel’s Bessel functions of the third kind of 
order v. 

The Hankel functions H“ (x) and H(x) bear the same 
relation to the Bessel functions J (x), Y,(x) as the functions 
exp(+ivx) bear to cos vx and sin vx, and are used in 
analysis for similar reasons. It should also be noted that 
the Bessel functions Y,(x), Hx), H{?(x) satisfy the 
same differential equations and recurrence relations as the 
function J ,(x). 


§ 31. Spherical Bessel functions. A problem which 
arises in mathematical physics is that of the solution of 
the wave equation in spherical polar coordinates 

2 
My, 2, 1 9/5, gay 
ol r ὃν r’*sin@ 40 \ a0 
lL δ 1284 
r* sin? @d¢? οἷ dt? 
If we take a solution of this equation of the form 


W = Yn, εἴθ, P(r)e™, (31.2) 


where Y,, ,(0, @) is the surface spherical harmonic defined 
by equation (23.3) and wW(r) is a function of r alone which 
satisfies the equation 


ay 2d¥ 5 }Ὁ 1) 
EEN i ἘΣ w / 
rm + a "2 nr or 0. (31.3) 
Now putting 
Y=r-4R (31.4) 


we see that equation (31.3) becomes 


ΔῈ 14 {wm = (n+4)’ 
er oa pate ol 
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whose general solution is readily seen to be 
R= Ad 4 ε(ων 6) Ἔ BJ _,~3(r]c). (31.5) 
Hence the function 
W=rFY,, (0, amy (Orie (31.6) 
is a solution of the equation (31.1). 

The functions Js (,+4)(x) which occur in the solution 
(31.6) are called spherical Bessel functions. We shall now 
show that they are related simply to the circular functions. 
First of all we consider the Bessel function J,(x). If we 


let ν = 4 in equation (30.6) and make use of the duplication 
formula for the gamma function we obtain the result 


+ es + o (=1)3"** 
1)= (2) ΣΦ, Grad 


which shows that 


2 
J,(x) = {| —] sin x. (31.7) 
MX, 
Again, if we put v = —4 in equation (30.6) we obtain 
the relation 
2.2 | 
J _4(x) = (2) cos x. (31.8) 
MX, 


The other functions J,,(x) where m is half an odd integer 
may be worked out in a similar fashion. It is left as an 
exercise to the reader to show that 


J m(X) = (2 2 { Sn(X) Sin χ -- ιν(χ) COS x}, 


J _{x) = (2 . 2\ (—1)""*{g,,(x) sin x +/,,(x) cos x}, 


where the functions /,,,, g,, are given in Table 1. 
I 
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These functions which arise in the way described have 
been tabulated in Tables of Spherical Bessel Functions 
2 vols. (Columbia Univ. Press, 1947) prepared by the 
Mathematical Tables Project of the National Bureau of 
Standards. 


Table 1 


§ 32. Integrals involving Bessel functions. In this 
section we shall derive the values of some integrals in- 
volving Bessel functions which arise in practical applica- 
tions. In the first instance we shall consider definite 
integrals. 

From equation (27.5) we have the relation 


Ϊ ΧΡ, - χιχ)άχ = [χυ (Χ)]. 
0 
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If n>0, x’J,(x)0 as x-+0 so that the lower limit is zero 
and we obtain the integral 


| x"J,—1(x)dx = αὐ, (α)}ὺ, (n>90), (32.1) 
0 
which, by a simple change of variable, gives the result 

| »,-κ(ξηγάν = LIfEa), (n>0). 022) 

0 
A particular case of this result which is of frequent use in 
mathematical physics is obtained by putting n = 1 in 
equation (32.2). In this way we obtain the integral 
| rJ o(Er)dr = τ κα). (32.3) 
0 

Further results may be obtained from (32.2) by familiar 
devices such as integration by parts. For example, making 
use of equation (27.5) we may write 


| " 8Jo(Er)dr = [ ? = © Grenjar, 
0 υ ¢or 


and, integrating by parts, we see that the right-hand side 
of this equation becomes 


3 a 
tse = | ry (Erde 
ξ ξ.}ο 
which reduces, by virtue of (32.2), to 
3 ,,.,,2 
© J6a)- Ξ J,(Ea). 


Now by the recurrence relation (26.2) we have the expression 


Ja) = 51 a(Ga)—JolEad 
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so that finally we have the result 
es 2a? 2%, 
| rJo(ér)dr = ἘΡΊΩΝ (308 = =) J ὯΝ (32.4) 
0 ν | ag 
Combining this result with equation (32.3) we obtain the 
integral 
a Pr: 
| ra? —r Ml Er)dr =F Jy(Ea)— ἜΣ Joba). (32.5) 
0 


The most commonly occurring infinite integrals are 
most easily evaluated by means of substituting the formula 
(27.3) in parts (ii) and (iii) of example 17 of Chapter II. 
From part (ii) of that example we see that 


a” oo 
PT W+)) [ of (ν Ἐ da*x*)x"*#e- *dx 


Γίμ-Ἐν- 1)α" 
2T(v-+ 1p τ] 


oF; (sudo $ut+dv+1;v+1; - 5) (32.6) 


If we make use of equation (30.8) on the left-hand side of 
this equation and of equation (7.4) on the right-hand side, 
we see that this result is equivalent to the formula 


= = T(u+v+1)a” 
J αχ) χε dx = —$<—<— ’ΞΕἃὯὯ'. 
; ΠΣ 2T(v-+1)(a? + pyres 
2 
αμην τ: ἦν--ἐμ; v+1; = :) (32.7) 
a“ +p 


where p>0, μὲ v>0. 
_The hypergeometric series occurring on the right-hand of 
this equation assumes a particularly simple form if either 
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a = vor # = v+1, and we obtain the formulae 


[ stoxwermde = TOAD ς 
0 


TG) | (a?+p*y*® (32.8) 


° πὰ TROD pa’ 
| J (ax)x"* *e~*dx = ra) Gia py 


Two special cases of the formula (32.8) which occur 
frequently are 


(32.9) 


on ree 1 
| JefasyeMax= Tae, G20 
‘a ΑΙ - | 
| xJ 1 (ax)e dx = (a2 +p» (32.1 1) 


Integrating both sides of equation (32.11) with respect to 
p from p to oo we find that 


oo 5 Ι Ρ 
J (axje""*dx = -- — —_——.. 32.12 
A special case of (32.9) which is often needed is 
= ~ px Ρ 
xJ o(ax)e™ *dx = —_——_.. 32.13 
| o( ) (a? + p*)* ( ) 


If we let p tend to zero on both sides of equation (32.7) 
we find that we can sum the hypergeometric series by 
Gauss’s theorem (7.2) provided that i ul<|v+1]. We 
then have the result 


ϑ τς ΖΡΓ(Σ ἘἐμΈ 3) 
| J (ax)x"dx = TG —Iu4d) (32.14) 


Similarly from part (iii) of example 17 of Chapter 1 
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we have the equation 


oo 
| oF (v+1; —4a?x?)e~ Px" t¥- Igy 
0 


- ΓΘμΈ3ν) Ὁ 
9 per ii 
which, because of (27.3), is equivalent to 


4u+4y; v+1; -- So 


2 
<) 


2 
| J (ax)e~?**x#~ dx 
0 


_ aT utd) 
ove 1 ptt ‘T(v 4+ 1) 


From parts (i) and (ii) of example 11 of Chapter II we have 
the special cases 


2 
iF; (su--4 v+1; -- } (32.15) 


"on v,~ 43,463 
Ee Shoals Pde ττοπτοὦὶ (Στ 
J, Bone rax = 2 2.16) 


and 


oe : ‘ 
νῈ3 Non pita, da a —g?/4p? 
| τ }, (αλ)ε dx = prema (?+1- oye eee 


0 
(32.17) 
of which the most frequently used are 


a 
= | 
xJ o(ax)e~?"*"dx = -" e~ 414 32.18 
[ 0 573 ( ) 


and 


x*Jo(ax)e~?*dx = _(1- a e~@/4* (32.19) 
0 2p* 4p? | 


§ 33. The modified Bessel functions. By an argument 
similar to that employed in § 1 we can readily show that 
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Laplace’s equation in cylindrical coordinates 
ὃν 1d) τὸν , oy 
do? 900 0° dg? az? 
possesses solutions of the form 
Ψ = et POA κι) 
where R(g) satisfies the ordinary differential equation 
d’R 1.4.8 2 “) | 
--- + -— —|m*+ —|R=0. (33.1 
ἀρ ρ dg ( Q? ay 
Writing x in place of mg we see that this equation is 
equivalent to the equation 
@R τὰκ ν᾽ 
— + --- -{1+-—]R=0. 33.2 
dx” χ ἂχ ( *) eed 
If we proceed in exactly the same way as in § 30 we 


can show that if ν is neither zero nor an integer the solution 
of this equation is 


R = Al ,(x)+ BI _,(x) (33.3) 
where A and B are arbitrary constants and the function 
I(x) is defined by the equation 


i= ἐδ οι 


2T(v+1) ro εἰν. 1), 
¥ 3 
------- 9 F,(v+1; 33.4 
ΓΕ) οὔιίν $x") (33.4) 
Comparing equation (33.4) with equation (30.8) we see that 
L(x) = i‘ (ix) (33.5) 


a result which might have been conjectured from the 
differential equation itself. 
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If v is an integer, n say, then I_,(x) is a multiple of 
I(x) so that the solution (33.3) in effect contains only 
one arbitrary constant. By a process similar to that 
outlined in § 30 we can show that in these circumstances 
the general solution of equation (34.2) is 


R= Al,(x)+ BK,(x) (33.6) 
where the function K,(x) is defined by the equation 
K,(x) = (—1)"""1,(x){log (4x) +7} 


<2 y (— oe r—1)! (4x) "+20 


r=0 


..1}9" +2r 
+4(-1) x aie Tn Ἐν aosiemule (33.7) 
The functions /,(x), K,(x) defined by equations (34.4) 
and (34.7) respectively are known as modified Bessel 
functions of the first and second kinds. 

The result (33.5) is very useful for deducing properties 
of the modified Bessel function J,(x) from those of the 
Bessel function J,(x). For instance, when n is an integer 
it follows from equation (25.7) that 


T_,(x) = 1,(x) (33.8) 
and from equations (26.1) to (26.5) respectively that 
21,(Χ) = Ly— (4) +I τ), (33.9) 
PE 1g(3) = Ips) Invi (33.10) 
xT (x) = xI,,_ .κ(χ)-- αἱ (χα), (33.11) 
ΧΙ (Χ) = nI,(x)+ x14 χα), (33.12) 
I(x) = 1,(x). (33.13) 
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Similarly equations (27.5) and (27.7) imply the relations 


£. {x"I,(x)} = x"I,~1(x) (33.14) 
dx 


{x~"I,(x)} = x7 "Tp a(x). (33.15) 


All of these relations can, of course, be derived directly 
from the definition (33.4) of J,(x) and it is suggested as an 
exercise to the reader to derive them in this way. 

It should also be observed that K,(x) satisfies the same 
recurrence relations as /,(x). 

The variation of J,(x), /,(x) and J,(x) with x is shown 
graphically in fig. 11 and that of K(x), K,(x) and K,(x) 
is shown in fig. 12. 


ἢ 34. The Ber and Bei functions. If we wish to find 
solutions of the form 


Ψ = R(g)e" 
of the diffusion equation 
aN Loy _ 1 oy 
Go* oo «Kat 


we have to solve the ordinary differential equation 

d*R Ἢ 1dR_ ἴω 

ἀρ edg κ 
On changing the independent variable to x = (ω κ)ξρ we 
see that this latter equation is equivalent to the equation 


dR τὰκ 
dx* x dx ae, 


Formally we may take the independent solutions of 


e~*1 (x) 


5 


-- ἢ --ὦ 


Fig. 12, Variation οἵ e*K,(x), e*K,(x) and e*K,(x) with x. 


$34 BESSEL FUNCTIONS 131 


this equation to be J,(i*x) and K,(i*x). Kelvin introduced 
two new functions ber(x) and bei(x) which are respectively 
the real and imaginary parts of J,(i*x), 1.6. 


ber (x) +i bei(x) = Ip(i*x). (34.2) 
From the definition (33.4) of J)(x) we see that 


05 am 2\2s 
ber(x)= Dy ἘΝ ΜΗ (34.3) 


Fig. 13. Variation of ber(x) and bei(x) with x. 


and that 
; ag (— 1)*(4x7)*5*! 
bei(x) = wo πτ 
1G) ἧς (2s +1)!? 
The variation of the functions ber(x) and bei(x) with x 
is shown diagrammatically in fig. 13. 
In a similar way the functions ker(x) and kei(x) are 
defined to be respectively the real and imaginary parts of 


(34.4) 
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the complex function K,(i*x), i.e. 
ker(x) +i kei(x) = K,(i*x). (34.5) 


oe the definition (33.7) of Ko(x) we can readily show 
t 


ker(x) = — {log (4x)+y} ber(x)+42 bei(x) 


oD 1 4r 
Ῥ ΕΣ τ φί(2ν), (34.6) 


Fig. 14. Variation of ker(x) and kei(x) with x. 
and that 


kei(x) = — {log (4x)+-y} bei (x)—4 ber(x) 
[:] (-- 1)'(4x)*"*? 
" » (2r+1)!? 


Fig. 14 shows the variation of the functions ker(x) and 
kei(x) over a range of values of the independent variable x. 


o(2r+1). (34.7) 
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The four functions ber, bei, ker and kei are used more 
often in electrical engineering than they are in physics or 
chemistry. For a full account of their properties and those 
of their generalisation to higher order and of their applica- 
tion to engineering problems the reader is referred to 
N. W. McLachlan’s Bessel Functions for Engineers (Oxford 
University Press, 1934). 


§ 35. Expansions in series of Bessel functions. We 
know from § 30 that 


ἢ ἧς +5 = = +(x 2_ "ἦν In(Ax)=0, (35.1) 


͵» £ +x — 2 442s} J (ux) = 0, (35.2) 
so that multiplying equation (35.1) by J,(ux)/x, (35.2) 


by J,,(Ax)/x, integrating with respect to x from 0 to a 
and subtracting we find that 


(A? — 11’) Ρ XJ (λχ), (μχ)άχ + (n? — m*) [ (Χμ, (μχ) a 
0 0 
= αἴ μι), (λα); (μα) --- }ῦ, (μαλ...(λα}}, (35.3) 
ifn>—1l,m>-—1. 
Putting m = n in this result we find that if 2 τὶ p, 
[Π xsaxsuods 
0 


5 Le (Aa) J,(ua)—AJ,(ua)J,(Aa)]. (35.4) 


ae 


The corresponding expression for A = p is obtained by 
putting μ = A+e, where e is small, using Taylor’s theorem 
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and then letting ¢ tend to zero. We find that 


| x{J,(ax)}? dx 


; 
= }q? | (ica)? (1- ΕΞ] {J,(2a)}?. (35.5) 


Suppose now that A and μὶ are positive roots of the trans- 
cendental equation 


hJ,(4a)+kAaJ{(Aa) = 0 (35.6) 
where ἢ and k are constants. It follows then that 
| χ (Χμ (μχ)άχ = ¢,0,. ἃ (35.7) 
0 


where 


cy = Way {k272a?+h?—k2n?}, (858) 


If we now suppose that we can expand an arbitrary 
function f(x) in the form 


f(x) = Σ ad (A,X) (35.9) 


where the sum is taken over the positive roots of the 
equation (35.6) then we can determine the coefficients a; 
as follows: Multiply both sides of equation (35.9) by 
xJ,(A;x) and integrate with respect to x from 0 to a; then 


|. xf(xJ,(ajxdx = Σ α; | ΧΙ (Δα), (λικ)άχ 
ΠῚ i 


0 
from which it follows that 


peels ᾿ Xl (Δ κ)άχ. (35.10) 
C; Jo 
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Because of its similarity to a Fourier series, a series of the 
type (35.9) is called a Fourier-Bessel series. 
In particular if the sum is taken over the roots of the 
equation 
Ji(ja) =0 (35.11) 


then the coefficients of the sum (35.9) are given by 
2A7 
{J {J,(4,a)}? (44a? 


Similarly if the sum is taken over the positive roots of the 
equation 


aaa |. xf (x)J,(Ajx)dx. (35.12) 


a;= 


J,(Aa) = 0 (35.13) 
we find that the coefficients a, are given by the formula 


a, = ———— |] x/(x)J,(A;x)dx. (35.14) 
ag 2 a*{ ἘΠΣ 2 a)? [ j 

In this section no attempt has been made to discuss 
the very difficult problem of the convergence of Fourier- 
Bessel series. For a very full discussion of this topic the 
reader is referred to Chapter XVIII of G. N. Watson’s 
A Treatise on the Theory of Bessel Functions, 2nd edition 
(Cambridge University Press, 1944). 


§ 36. The use of Bessel functions in potential theory. 
As an example of the use of Bessel functions in potential 
theory we shall consider the problem of determining a 
function w(o, z) for the half-space azo20, z20 satisfying 
the differential equation 


δ'ψι Lop , dy 
de” @6Q az? vA} 
and the boundary conditions: 
(i) YW = Κρ), on z = 0; 
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(ii) w0 as z> 0; 
τ ὃ 
(iii) id +x =O0oneg=a; 
6@ 
(iv) w remains finite as 00, 


We saw in § 1 that a function of the form Ww = R(@)Z(z) 
is a solution of equation (36.1) provided that 


a’*Z 


ἜΞ —42Z =0 (36.2) 
and that 
2 
= + ΠΡ +/2R =0 (36.3) 
dg’ do 


where A; is a constant of separation. To satisfy the 
boundary condition (ii) we must take solutions of equation 
(36.2) of the form 


Ζεε 
and to satisfy the condition (iv) we must take as the 
solutions of equation (36.3) functions of the form 

R = 70(λι0) 


since the second solutions Y,(A,0) would become infinite 
in the region of the axis 9 = 0. 

The differential equation (36.1) and the boundary con- 
ditions (ii) and (iv) are satisfied by any sums of the form 


ψίρ, 2) = » ae” ΛΑ] (4,0) (36.4) 


where the a; and /, are constants. But if we are to satisfy 
the boundary condition (iii) we must take the sum over 
the positive roots of the equation Ὁ 


iJ (Aya) + KJ (Aja) = 0. (36.5) 


+ For properties of the roots of this equation see example 16, p. 144. 
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The solution is determined therefore if we can find constants 
a; such that condition (i) is satisfied, i.e. such that 


ΚΟ) = Σ, αμ]ο(λι0). (36.6) 


From equations (35. = and (35.8) we see that we must take 


= sare a Ἢ o'fleWole’\de". (367) 


Hence és required solution is 


2 αὐ ie | “Με 
W(p, 2) = ἘΣ, (A? +x*){ J o(A,a)}? ice at (36.8) 


where the sum is taken over the positive roots of the 
equation (36.5). 

If, instead of the boundary condition (iii), we had the 
condition wy = 0 on g = a then it is easily seen that the 
solution would have been 


e”**#Jo(4i0) [“ 
,. ον (λιογάρ' (36.9) 
#2 Taal? {J(Aa)}? a 
where the sum is taken over We positive roots A; of the 
transcendental equation 
J (Aa) = 0. (36.10) 
For example suppose that w satisfies the conditions 
yw =0 on g@=4, ¥0 as z>0, Ψ -- Ψψρ(α"--ο7) on 
z = 0, 0SeSa, then the solution to the Bay is given 


by equation (36.9) with f(o’) = ψο(αῦ -- ο΄). By equation 
(32.5) we have 
| o’fle’o(e’)de’ = “FJ (a) ame Tol) 
0 
- "9 "Δα. 


138 SPECIAL FUNCTIONS OF PHYSICS AND CHEMISTRY § 37 


since A; is a root of equation (36.10). Thus the required 
solution is 


, - ὅϑψο e~** J (1,0) 
MY BT (aa) 


Tables of the first forty zeros €; of the function J,(é) 
with the corresponding values of J,(é;) are available + so 
that it is convenient to express results of the kind (36.11) 
in terms of them. It is readily seen that in this case 


2) = 8a? e~ J 4(a€;) 
Me Re Ke a. ae eS 


(36.11) 


(36.12) 
where { = z/a and ἃ = o/a. 


§ 37. Asymptotic expansions of Bessel functions. In 
certain physical problems it is desirable to know the value 
of a Bessel function for large values of its argument. In 
this section we shall derive the asymptotic expansion of 
the Bessel function of the first kind J,(x) and merely 
indicate the results for the other Bessel occurring in mathe- 
matical physics. 

We take equation (28.1) as our definition of the function 
J,(x). Applying the theory of functions of a complex 
variable it is readily shown that this definition is equivalent to 


εὐ. _ ttt | eer 
"0 το τῶ ie Tn Ore dats we eat} 


(37.1) 


where L, is the straight line A(t) = —1 in the upper 
half of complex f-plane and L, is the corresponding part 
of the straight line A(t) = +1. By changing the variable 

+t A. Gray, G. B. Mathews and T. M. MacRobert, A Treatise 


on Bessel Functions and Their Applications to Physics, 2nd edition 
(Macmillan, 1931). 
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from ¢ to u = ix(l—¢) in the first integral and to 
u = --ἰχί! --- ἢ in the second we see that 


Je) = | = τ πο (37.2) 


: ; 1 ix—(4n+4)zi |. —u,n-4 ( ἫΝ 
(x) = ------- e e~“u 1 --ἡ du 
iW) = τα τῇ , 2χ 


and j.(x) denotes its complex conjugate. Expanding 
(1+iu/2x)"~* by the binomial theorem and integrating 
term by term we find that 


lx) = fe*~ G0" Y™, Fy ὩΣ Σ--η; x) (37.3) 
ix 
If we adopt Hankel’s convention of writing 
-- r(4—"),G+n), 
(n, r) =(—1) a ies 
in equation (37.3) and substitute the result in equation 


(37.2) we find that for large values of x the asymptotic 
expansion of the Bessel function J,,(x) is 


J(x)~ [{ {eos (x—4nn—4n) Ἢ ως, —— 


—sin(x—4inx—4n) : a = eth (37.4) 


The corresponding expansion for the Bessel function of 
the second kind is found to be 


Y, (x) ~ ‘2 sn (x—4nn—4n) Σ (- 1), 27) a 


+008(x—4nn— ἐπ) Σ᾽ (He seal (37.5) 
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Substituting these asymptotic expressions in equations 
(30.18) and (30.19) we find that as x00, 


+ + 
ΓΝ (2) je), H2)~ (2) i) 616) 
TX tx 
where j,(x) is given by equation (37.3). 
In certain problems only a very crude approximation 


to the behaviour of the Bessel function is desired. In these 
circumstances the following formule are usually sufficient: 


J,(x)~ [= cos(x- 4nn—4n), 
¥,(x)~ [3 sin(x—}nn—4n); (37.7) 
TX 


HOXx)~ 2 gt ὐσπι- ἐπὶ 
\) mx 
HONx)~_ [2 e~*+imi-481, (378) 
TX 


Similar formule exist for the modified Bessel functions. 
Proceeding in the same way as in the establishment of 
equation (28.1) we can show that 


cet td dec 
I(x) = πατῇ (χ}" | τ 
which becomes 


= Leyte: Se fae Ὁ τὰς » εἰ ἫΝ 
I(x) πα Γ(ι 4) ἰ Ξ é€ 1+ ox du 


verona) 


by a simple change of variable. By a method similar to 
that employed above to obtain the asymptotic expansion 


et (1 — 22)" 41 
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of J,(x) we can then show that if —4nx<argx<4z, 
ἔξ (=r), eH ὁ (7) 

ΤῸ mx) r=0 (2x) (xx) r’=0 (2x) 


and that if —3n<arg x<4n the factor exp {—x+(+-4)zi} is 
replaced by exp { — x—(n+-4)zi}. Thecorresponding formula 
for the modified Bessel of the second kind is 


_(x (n, τ) 
bia é ye "eae (2x)' 


I,(x)~ 


as X¥- 0. 


Examples IV 


(1) Making use of example 2 of Chapter II and of 
the expansion (25.8), expand cos (x sin θ) as a power series 
in sin@ in two ways. Hence by equating powers of 
sin?* @ show that if s is a eT yaar 


Derive the corresponding result for ane and show that 
the two results may be combined into the single formula 


aS CRNA (r+ 2m 22) 
(ν τα 1,2,3,...... 

(2) If uw and ζ denote the eccentric and mean anomalies 
of a planet show that 


cos(nuj=n δ᾽ ~ m-n(me) cos (mf), 


sin (nu) =n Σ 1 J m=n( ime) sin (mf). 
m 


n= -c¢ 
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(3) Making use of the expression for P, (cos 0) given 
in example 19 of Chapter ILI show that 


- K [ cos ἢ i 
2, = (cos 0) = 6" 8° J,(r sin 8). 


(4) Show that 

(i) ϑυ (2) = ,. 3(2)— 3... 12) +3) n41(2)— Ins 32); 
(ii) 4J9(z)+3J9(z)+J3(z) = 0. 

(5) Prove that 


a gas _ (-1)* [*cos(N+4u ὧδ 
gs α PILOTS π | cosdu  ,/(x?—u?) 


and deduce that 
: + Σ, (—1)'Jo(rx) = 0. 
(6) Show that the curve with freedom equations 


x =f-—sint, y = cos ¢/(1—cos ἢ) 


may be represented in the interval 0<‘<z by the Fourier 
series 


yas ΕΣ J,(n) cos (nx). 
(7) Prove that 
elt con (=) a (2n+ le" J, . ,(kr)P,(cos 8). 
(8) Prove that 


n=? oO 


lim P, (cos *) εξ J,(x). 
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(9) Show that 


(i) Σ = Ja) = Jo{./(a?—2ax)}: 


ii) ¥ = απο, (2a) = χα τα) τ, f2/(x-+a)} 
a=0 ἢ: 
(10) Prove that if m is a positive integer 
(2/3) = (- γ᾽ x" Z Jo). 
(11) If x>a show that 


ἢ e* 5 cos (x sin θ)εθ = Jo{,/(x?—a’)}. 


0 
(12) Prove that, if -—l<x<l, 
1 oe 
————— =4n+n J (mm) cos (mmx). 
Tana 7 HE, ola) cos (max) 
Deduce that 


Jats) = B21 428° ἢ eee 


m=1 x?—m?n* | 
(13) Prove that 
F Π] Α 
ΘΝ - (x) ΓΕ ἡ αν (x Ξ a 

where A is a constant, and, by considering the series for 
J (x) and J,(x) when x is small show that A = (2/7) sin (yz). 
(14) Show that the complete solution of Bessel’s 

equation may be written in the form 


AJ,(x) + BJ,(x) Γ 


dé 
E{J,(@)}? 


where A and B are arbitrary constants. 
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(15) Show that the complete solution of the differential 
equation 


is 
y = /x{AJ,(€)+ BJ_4(®}, 

where €* = 4x°/27 and A and B are arbitrary constants. 

(16) If a and 6 are real constants show that the roots 
of the equation 

axJ}(x)+ bJ,(x) = 0 

are simple roots except possibly the root x = 0. 

Show also that the equations J,(x) = 0, J,(x) = 0 have 
no roots in common except possibly x = 0. 

(17) If x>1 and m+n+1>0, prove that 


ἷ πὰ ΚΑ ΟἸ ν = iE (x? —1)~2"O7(x), 
Ἵ] 


where Οὐχ) denotes the associated Legendre function of 
the second kind. 
(18) Prove that 


L(x+y)= Σ 1,Ἱ],- Ὁ) 


0 
+ Σ {1 ,ς, ) ἘΠῚ (Χ)],)}. 


(19) Show that 


1 : 
| e*P (μ) ἀμ = [Ξ 10). 
«a4 x 


(20) Prove that 
Σ σὴν = εξ Σ (x) 
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and deduce that 
@ Ire) Σ Jpan(eerew sin 
m= ἢ ἢ]. 


τ" oo ἫΝ 

(il) I(x) — 2, ml Intnl). 

(21) Using the expansion of the last question prove 
that J,(ae'+ be'*) is the coefficient of ΖΓ in the expansion of 


exp - esinartsn δὶ , Σ εἰνατίκβ} (a) (byt 


t 


n=-oOn= - 
By putting R= acosa+bcosf, 0=asina+bsin f, 
B = «+0—nx prove Neumann’s addition theorem 


— howe \an oo 
JR) = tas ΟΣ “μενα, (6) "", 


where R? = a* +b? —2ab cos 0. 
(22) Prove that 
Jo(ax\lo(bx) -- |. Jo(Rxd0, 
' 0 
where R? = a*+b?—2ab cos 0, and deduce that 
3 J o(ax)J o(bx)e~ “dx = pee K(k) 
0 τι (ab) 

where 


κ(ὸ = Si ee 


r Delos, wigtsg act, 
0 V(—k? sin? $) (a+b)? +c? 
Prove, in a similar way, that 
2 


| J \(ax)J ,(bx)e" “dx = ak (ab) {(1—4k*)K(k) — E(k)} 
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where k and K(k) are as defined above and 


E(k) = | Ahn 8? ein? Ode 
0) 
(23) Show that 


χα : eXn—mai |” 7 (Χ sin Bye *™ do 
ῦ 


and hence that 
{I,(x)}* = an J»,(2x sin 0)dé. 


Deduce that 


τ Fates! 


θα οι δι ἡ πεν 
> slants)nts)! ° ΣΧ . 


(24) Prove that 
2 
ἢ Jo(ax)J o(bx)e~?*’xdx = » exp ( ~ a*—b >) I (32) 


4p? 2p 


(25) We define the Bessel-integral function of order n 
by the equation 


Jilin) « -| 2) ay n't, 2, ἃ. 
one 
Prove that if m is even 
jitoa* | cos (nd)ci(x sin 0)d0, 
o 


where ci(x) denotes the cosine integral, and derive the 
corresponding expression when ἢ is 
Show that: 


(i) Life = SY Pix); 


ΝΞ “οὐ 


(ii) Ji,(x) = J,(x)/x; 
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(iit) (a — 1)Ji,— (x) —(n + IJ i, 4 4) = 2nJ I(x); 
(iv) ci(x) = Jig(x)—2Ji,(x) + 2Ji4(x)—...; 
(v) si(x) = 2Ji,(x)—2Ji,(x)+2Ji,(x)—... 
(vi) Jig(x) = y +log (4x) — = Fill, 1; 2, 2,2; —x?/4). 
(26) If 
= x Ὁ. 2m+1 | 
“Ta ee 


where A is a zero of Jp(x), show that 


4m? 4*m*(m—1)* 
ΠΣ re + are te 
—1)?...22.42 
Show that 


i—x™= ¥ aJ,(4,x), 0sxsi, 
i=1 


where 4,, A2, ... are the positive zeros of J,(x) and 


| 1-—J] 
“as saat " 


Putting m = 1, 2, 3, ..., and rearranging the series show 
that 


JoAxX) 1, _ γ2. 
Bi) 8 Ὁ 


oe Δ 1 aya any, 
2, Bi)” 126° C-*) 


TolArx) - ὃ “-«ὩΝ10.. νυ. x4 
2 77,0) - 7 x*\(19 — 8x* + x"). 
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Deduce that 


eS 4 oe | ee ae 19 
Σ 


Β ei DEEL ER σὺ, Sie? ἔνι. eee 
2%, A? J ,(A;) g” hi APJ) 128’ τοι AIA) 4608 


(27 If a, a, ... denote the positive zeros of J,(x) 
show that 


. a ἘΠῸΝ [. xf(x)J (ax)dx = 4 | x” *f(x)dx. 
i= 1 Od y41(%) Jo 0 


Taking f(x) = x’, x"** show that 
S102 4(v+1) ist αὐ 16(v+1)*(v+2) 
(28) Putting ν = 0, f(x) = x” in (27), show that 
= I, 1 
2, A, 4m+1) 

where /,, Az, ... are the positive zeros of Jo(x) and J,, is 
defined in (26). 

Substituting the value for J,, derived in (26) and using 
the fact (obtained by putting ν = 0 in (27)) that ¥ 1 = 4, 
show that the sum 


= 1 
Ἐπ Σ, _ 


| satisfies the recurrence formula 


ς Ἢ y συ Se 
ee ee 4 [m—r+)!7 


“ir 
4 m\(m +1)! 


Shoat set ὁ Sa τὰ 
Σ᾿ ae 535 Py te ey me, λδ 12288 
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(29) Show that 
a ee 
2 Bra) 218 


Multiplying both sides of this equation by x and integrating 
term by term, show that 


> Fix) = 4x(1—2 log x), 


11 ARITA) 
and hence that 
ὦ ἦ 
2 Bray? 


CHAPTER V 


THE FUNCTIONS OF HERMITE AND LAGUERRE 


ἢ 38. The Hermite polynomials. The Hermite poly- 
nomial H,(x) is defined for integral values of ἡ" and all 
real values of x by the identity 


oo ςς x a ”. (38.1) 


If we write 
f(x, t) = e2tx-? me ee 


then it follows from Taylor’s theorem that 


ΟὟ ee 
A(x) = | — =| ee ee 
“ (Z) ἜΣ: |. 


Now it is obvious from the form of the function 
exp {—(x—1)”} that 


ὃ' —(x—t)? Ee on » oe ~x? 
KE ΜΒ a eos a 


and so we have the form 
Beg eer | | 
Hf = | --Ἠ1 ‘a —— a sa * 
Xx) Ξε (-- ἡ a (e*) (38.2) 


for the calculation of the polynomial H,(x). 

It follows from this formula that the first eight Hermite 
polynomials are: 

H)(x) = 1, 

H,(x) = 2x, 
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H,(x) = 4x? —2, 

H,(x) = 8x*—12x, 

H,(x) = 16x*—48x? + 12, 

H,(x) = 32x*—160x* + 120x, 

H (x) = 64x° —480x* + 720x? — 120, 

H (x) = 128x’—1344x* + 3360x* — 1680x. 
In general we have 


H,{x) = (2x)"— nD ax? 


n eee (2x)"- τ... 
or, in the notation of Section 12, 
H,{x) = (2x)" 2 Fo (a 4—4n; -- 5) (38.3) 


Recurrence formule for the Hermite polynomials follow 
directly from the defining relation (38.1). If we differ- 
entiate both sides of that equation with respect to x we 
obtain the relation 


γιο2 χιπτι — ¥ A, (x) {" 
n=o n! 
from which it follows directly that 
2nH,_ x = Hi(x). (38.4) 


On the other hand if we differentiate both sides of the 
identity (38.1) with respect to t we obtain the relation 


-- 2ικχ--ἰ3 τὸν < H,{x) -1 
a8 Ὡς π᾿ 


which can be written in the form 


Ὁ H,{x) n μι Η, (ΧῚ πὉ'ὁὌο. < H,{x) - 
- ὌΝ ΝΠ ors ΝᾺ a Ae ὧν (n—1)! τς 


152 SPECIAL FUNCTIONS OF PHYSICS AND CHEMISTRY § 39 
to yield the identity 

2xH,(x) = 2nH, ~ 1(%) + Hy + 1(%) (38.5) 
by the identification of coefficients of £". 


Eliminating 27H, _ ,(x) from equations (38.4) and (38.5) 
we obtain the relation 


H,(x) = 2xH,(x)— Η, + 1(x). (38.6) 
Differentiating both sides of this identity we find that 
Hy (x) = 2xH,(x) + 2H, (x) — Hy +(x) 
and, by equation (38.4), H,+,(x) = 2(n+1)H,(x) so that 

Hi"(x)—2xH(x)+2nH,(x) = 0. (38.7) 


In other words y = H,(x) is a solution of the linear differ- 
ential equation 
y"~2xy'+2ny = 0. (38.8) 


§ 39. Hermite’s differential equation. We saw in the 
last section that H,(x) is a solution of the differential 
equation (38.8). Replacing the integer nm in that equation 
by the parameter vy we obtain Hermite’s differential equation 


d"y dy : 
—— —2x — +2vy = 0. 39.1 
dx? dx ii iss 
If we assume a solution of this equation in the form 
y= Σ α,χ 
γι ῷ 


and substitute in the equation (39.1) we obtain the recur- 
rence relation 
2(r+e-—v) 
a,4.2 = ———————_ 4, (39.2 
2 (r+0+2r+o+1) iio 


on equating to zero the coefficient of x’**. Equating to 
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zero the coefficient of x*~* we obtain the indicial equation 

o(e@—1) = 0. (39.3) 

Corresponding to the root 9g = 0 we have the recurrence 

relation 

wf I 

(r+1)(r+2) * 
which gives the solution 


(39.4) 


Gyp42 


oe a oe =v 
yi(x) mig; (1- ἣν CO) OW 564...) 


2! 4! 6! 


(39.5) 


where a, is a constant. 

Similarly, corresponding to the root g = 1 of the indicial 
equation we have the recurrence relation 
- 2Arti-v) a, (39.6) 

(r+3)(r+2) 
from which is derived the solution 
y2(x) = Asx (: = ἘΠῚ x? + ven) x*+ 7,92 
where a, is a constant. The general solution of Hermite’s 
differential equation is therefore 
y = yy(x) + yo). (39.8) 


For general values of the parameter v the two series 
for y,(x) and y,(x) are infinite. From equations (39.4) 
and (39.6) it follows that for both series 


ap43 


G,42™ 5 d,, a5 roo. (39.9) 
Γ 


If we write 
exp (x?) = bg +b2x7 +... +b,x" +b, 42x"77? +... 
L 
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then 
ἜΚΕΣ = AS αὐτὰ (39.10) 


Suppose that ay/by is equal to a constant y, which may be 
small or large, then it follows from equations (39.9) and 
(39.10) that, for large enough values of N, dys. m/bnsom™): 
In other words the higher terms of the series for y,(x), 
y2(x) differ from those of exp (x) only by multiplicative 
constants 7, 72» So that for large values of | x |, 


yiX)~ne™,  yal(x)~72€* 
since for such values the lower terms are unimportant. 

We shall see later (§ 41 below), that in quantum 
mechanics we require solutions of Hermite’s differential 
equation which do not become infinite more rapidly than 
exp (4x*) as |x|-0o. It follows from the above con- 
siderations that such solutions are possible only if either 
y,(x) or y2(x) reduce to simply polynomials and it is 
obvious from equations (39.5) and (39.7) that this occurs 
only if ν is a positive integer. For example, if ν is an 
even integer n we get the solution 


y(x)= cH,(x), (39.11) 
where ¢ is a constant, by taking a, = 0, 


7 n! 
a, =(—1)"——-e. 
(fn)! 
Similarly, if v is an odd integer n we get the solution 
(39.11) by taking a, = 0 and 


a, = (—1)#-# 55’ 


(4n—})! 
Hermite’s differential equation therefore possesses solutions 
which do not become infinite more rapidly than exp (4x?) 
as | x | oo if and only if ν is a positive integer n. When 
this is so the required solution of Hermite’s equation is 
given by equation (39.11). | 


Cc. 
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ἢ 40. Hermite functions. A differential equation closely 
related to Hermite’s equation is 


2 
ἐν (-- χὴν τῷ. (40.1) 
dx? 
If we transform the dependent variable from wy to y where 


=e *'y (40.2) 


and put A = 1+2y then it is readily shown that y satisfies 
Hermite’s equation (39.1). The general solution of equation 
(40.1) is therefore given by equations (40.2) and (39.8) 
with y,(x), y2(x) given by equations (39.5) and (39.7) 
respectively. 

The argument at the end of the last section shows that 
the equation (40.1) possesses solutions which tend to zero 
as |x [- τοῦ if and only if the parameter A is of the form 
1+2n where n is a positive integer. When A is of this 
form the required solution of (40.1) is a constant multiple 
of the function ‘Y,(x) defined by the equation 

P(x) = e #°H,(x) (40.3) 
where H,(x) is the Hermite polynomial of degree n. The 
function ‘Y,(x) is called a Hermite function of order n. 

The recurrence relations for ‘Y,(x) follow immediately 
from those for H,(x). For instance equation (38.4) 15 
equivalent to the relation 


2n'¥,-1(x) = ΧΟ +P.) (40.4) 
and equation (38.5) is unaltered in form so that 
2x'P (x) = 2n'P,,— (x) Ἐν]. τί). (40.5) 


Eliminating 2n'Y,,_ ,(x) from equations (40.4) and (40.5) we 
have the relation 


Ψ,(χ) = ΧΨ, ()--, 10). (40.6) 
From the point of view of mathematical physics the 
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most important properties of Hermite functions concern 
integrals involving products of two of them. In establish- 
ing most of these properties the starting point is the 
observation that the function ‘Y,(x) satisfies the relation 
YP", +(2n+1—x*)¥, = 0, (40.7) 
as is obvious merely by substituting 2n+1 for 4 in equation 
(40.1). 
Writing down the corresponding relation for Ψ',,, 
yr +(2m+1—x7)¥,, = 0, (40.8) 


multiplying it by Ψ,, and subtracting it from equation 
(40.7) multiplied by ‘Y ,, we obtain, as a result of integrating 
over (— οὐ, 00), the relation 


χα) [ WW dx = [ (P,P — BW"), 


Now an integration by parts shows that the right-hand 
side of this equation has the value 
oa oo oo : 
| eas -| ΟΡ, PE ax 
“οὖ =o 


and, if we remember that, for all positive integers n, 
Ψ,()-0 as | x [- τοῦ, we see that this has the value zero. 
Hence if we let 


i ὦ [ ΨΟΟΨ,Ολάν 
we see that 
In, « = 0, if m #7. (40.9) 


In particular 
I,- i,n+1 = 0 


so that from equation (40.5) we have 


ἧ 2x, (0), -(x)dx = 2nl,1,n-1- (40.10) 
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Now if in equation (40.3) we substitute for H,(x) from 
equation (38.2) we have 


W(x) = (-- 53» ᾿ (e-*) (40.11) 
x 
so that the left-hand side of equation (40.10) is equal to 
Ἢ τ πὸ ὦ 
- ixe® Ξ- te" *)d. 
€ £ eax 


and an integration by parts shows that this is equal to 
Fa, nt Tne, καὶ 
i.e. to J, ,. Hence from equation (40.10) we have 
᾿ τ 
Repeating this operation n times and noting that 


lo,o = ᾿ 6. ἀχ = Jn 


we find that + 

Ty, np = 2'n!,/n. (40.12) 
Combining equations (40.9) and (40.12) we have finally 

In, » = 2"! [5 pan: (40.13) 


The evaluation of more complicated integrals can be 
effected by combining this result with the recurrence 
formule we have already established for the Hermite 
functions, For instance, it follows from equation (40,5) 
that 


| XP (x)P(x)dx τα nlp. 514m. καὶ 
showing that 


| ” B_(x)¥,(x)dx =Oifm¢ént1 (0.14 
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and that 
| XY,(x) Yaa (x)dx = 2"(n+1)!/x. (40.15) 
-ω 

Similarly, making use of equation (40,4) and equations 

(40,.13-15) we can show that 
be OifmAnt+1 
| (xVM x)dx = 25 nlx if m=n—1 
“=? —2"(n+1)!./x if m=n+1., 


§ 41. The occurrence of Hermite functions in wave 
mechanics. The Hermite functions which we have dis- 
cussed in the last section occur in the wave mechanical 
treatment of the harmonic oscillator.t Although this is 
a very simple mechanical system the analysis of its properties 
is of great importance because of its application to the 
quantum theory of radiation. 

The Schrédinger equation corresponding to a harmonic 
oscillator of point mass m with vibrational frequency ν is 

d*y ἐπ tn 
dx re ΟΝ 
where W is the wi energy of the oscillator and A is 


Planck’s constant. The problem is to determine the wave 
functions wy which have the property that 


(i) y0 as | x |>00; 
(ii) [Ψ [dx = 1. 


If we let 


(W—2n?mv?x2)~ =0, (41.1) 


1 [ἢ 


x=— | — 


2n 


ἘΝ, F. Mott and I. N. Sneddon, Wave Mechanics and Its 
Applications (Oxford, 1948), p. 50. 
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then the equation (41.1) becomes 
a"y 2 
0 : 
": ε( -2 )v- (41.2) 
and the conditions (i) and (ii) become 
(1 0 as | € |>00; 


(ii’) [ [Ψ [βαξ = 2π Ἐ 


The argument given at the beginning οἵ ὃ 40 shows 
that equation (41.2) possesses solutions w which satisfy 
the condition (ii’) if and only if the parameter (2W/hy) 
which occurs in the equation takes one of the values 
1+2n where ἢ is a positive integer. In other words 
solutions of this type, which are known by the probability 
interpretation of the wave function y to correspond to 
stationary states of the system can exist if and only if 


W= hv(n+}) (41.3) 


where ἢ is a positive integer. When this is the case the form 
of the wave function w is known from ὃ 40 to be 


wy = CY,(¢) (41.4) 


where C is a constant. Applying condition (11) and 
equation (40.12) we see that 


πὸ ‘am 1 
NR) apt 


Thus the wave function corresponding to an admissible 
energy (n+4)hy is 


4nmy . Y(6) an’ my 
Vn ( i nie C= 2n [ax (41.6) 


In quantum theory the matrix elements (n | x | p) 
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defined by the equation 
Leliam | Sita bias 


are of considerable importance in the case of the harmonic 
oscillator. In terms of the variable € we have 


(n|x|p)= ΩΣ, Ἢ SWS) (ode, 


so that substituting from equation (41.6) we have 


(n|x|p) = Ἢ — ἢ ἐψ (ξ)Ψ (ξ) ἀξ + {2 * 27(n!)*(p!) 3}. 


It follows then from equations (40.14) and (40.15) that 


acs (n|x|p)=0 ifp#ntl (41.7) 
and that 
ih) + 
(n|x|n+1)= ἘΞ 3 (41.8) 
nh 1} | 
(n | x | "- 1) ΞΞ ἘΠ . (41.9) 


ὃ 42. The Laguerre polynomials. The Laguerre poly- 
nomials L,(x) are defined for m a positive integer and x 
a positive real number by the equation 


exp (- ΞῚ = a=) Σ, on ". (42.1) 


Expanding the exponential function we see that 


eee (ie i ee et 
oe exp( £0 ri(l—ty** 


1-—t 
= ᾿ς Ὺ (-θὙ 1} pts 
r=0s=0 ris! 
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The coefficient of ¢" in this expansion ts 
(<1 + Dane 
r=o ri(n—r)! 
Using the relations 
! (-1% (-n»)r 
r+1),_,= = ——— = ---- 
Sil r! (n—r)! n! 
we see that this sum can be written in the form 
(—"), 
Σ ae (r!)? 
Identifying the coefficient of κ᾿ with L,(x)/n! and adopting 
the notation of § 11 we see that 
L,(x) = n! ,F,(—n; 1; x). (42.2) 
It should be observed that L,(x) is a polynomial of degree 
nin x, and that the coefficient of x” is (— 1)”. 

A useful formula for the polynomial L,(x) can be 
obtained by finding a new representation for the confluent 
hypergeometric function on the right-hand side of equation 
(42.2). By Leibnitz’s theorem for the nth derivative of a 
product of two functions we have 


e*D"(x"e"*)=e* Ὁ (1) <0 (Dere(D'e™) 
r=0 Γ᾽ 
where D denotes the operator d/dx. Using the relations 


e*D'(e~*)=(—1), δ χ" =n!x'/r! 
we see that 


e*D"(x"e"*) = n! Fs cp? my (xy. (42.3) 
It follows immediately from equation (42.2) that 
L(x) = & -- i! | (42.4) 
x 
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The first five Laguerre polynomials can be calculated 
easily from this equation; we find that 

Lo(x) =I, 

L,(x) bes l—<x, 

L(x) = 2-—4x+x?, 

L3(x) = 6—18x+9x? —x’, 

L4(x) = 24—96x+72x? — 16x? + x*. 

Equations (42.4) can be used to show that the functions 


φ,(Χ) = — =e L(x) (42.5) 


form an orthonormal system. From (42.4) we have as a 
result of mm integrations by parts 


|e  € Libis « i x” (xe “τὰν 
Jo 


= (—1)"m! [ΖΞ qm 


and this is zero if n>m. Since L,,(x) is a polynomial of 
degree m in x it follows that 


ἱ: e7*L,(x)L,(x)dx = 0 if m #n. (42.6) 
0 


Since the term of degree nm in L,(x) is (—1)"x" it follows 
that, when m = n, 


Ν" εὐ) βάχ = (—1)" ᾿ e-*x"L,(x)dx 
0 0 


= | nix"e”-*dx 
0 


= (n!)? 
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Combining this result with equation (42.6) we find that 


| ~ Oeil A (42.7) 
αἱ ἢ 


showing that the ᾧ᾽5 form an orthonormal set. 

Recurrence formule for the Laguerre polynomials may 
be derived directly from the definition (42.1). Differ- 
entiating both sides of this equation with respect to f we 
obtain the identity 


~G-) ee yee 


n=o (n—1)! azo an! 
which may be written in the form 
ae: 20 Lot" 
‘ Sara mee Σ, (n— 1)! 
ὟΝ L(x)t" _ 
-(1- Σ ee 


n=0 


Equating to zero the coefficient of ¢” in the expansion on 
the left we obtain the recurrence relation 


Ly (xX) +(x—2n—1)L,(x)+n7L,_ (x) =0. (42.8) 
Similarly if we differentiate both sides of (42.1) with respect 
to x we obtain the identity 


ΕΣ πρὸ pr 4-(1 t) τ: =a pp 


which yields the recurrence relation 


L(x) — nL, ~ (x) + nL, - (x) = 9. (42.9) 
Differentiating equation (42.8) twice with respect to x and 
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replacing ἡ by n+ 1 we find that 
Di +2(x) +(x —2n —3)L,, + (x) 
+(n+1)*Li(x)+2L,4,(x) =0. (42.10) 

Now from (42.9) 

Πα) = οι 1) 1.0) -L,09}, 
and hence 

D+ α(Χ) = (n+ 1{ L(x) — L009}. 


A similar expression for I,,.(x) in terms of L,(x) and its 

derivatives can be readily obtained. Substituting these 

values of L,.,(x) and L(x) in equation (42.10) we find 
that 

xL(x)+(1—x)L,(x)+nL,(x) = 0. (42.11) 

ἢ 43. Laguerre’s differential equation. Equation (42.8) 


shows that y = AL,(x) is a solution of Laguerre’s differ- 
ential equation 


x2 .ᾳ- x2 = x+y =0 (43.1) 
in the case in which ν is a Sih integer n. If we put 
y = 1, « = —v in equation (11.2) we see that it takes 


the form (43.1) so that it follows from equation (11,4) 
that one solution of equation (43.1) is 

y;(x) = ,F,(—v; 1; x), (43.2) 
Similarly we see from equation (11.8) that the second 
solution is 


y2(x) = γι). . log x-+ ᾿ ox! (43.3) 


where the coefficients c, are defined by equations of the 
type (11.9). The general solution of Laguerre’s differential 
equation may therefore be written in the form 


y = Ay,(x)+ By2(x) 
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where A and B are constants and y,(x), v(x) are defined 
by equations (43.2) and (43.3) respectively. 

If we are interested only in solutions which remain 
finite at x = 0 it is obvious from equation (43.3) that we 
must take the constant B to be zero. Further if a, is the 
coefficient of x” in the series expansion for y,(x) it is easily 
shown that a,,,/a,~r~*. As the result of a discussion 
similar to that advanced in ὃ 39, it follows that if ν is 
not an integer 

γιμι(χ)» εἰ as x00. 
If, therefore, we are looking for solutions which increase 
less rapidly than this we must take v to a positive integer, 
in which case y,(x) reduces to a polynomial. 

The equation (43.1) possesses a solution which increases 
Jess rapidly than οὐ as x-0o if and only if the parameter v 
occurring in it is a positive integer, n say. Ifit is also required 
that the solution shall remain finite at x = 0, the solution 
is of the form 

y = AL,{x) (43.4) 
where A is a constant and L,(x) is the Laguerre polynomial 
of degree ἡ. 


§ 44. The associated Laguerre polynomials and func- 
tions. If we differentiate Laguerre’s differential equation 
m times ni respect to x we ἫΝ that it becomes 


ἐξ. nt  +(n—m) ΤΣ ᾿ 


Ἁ dx m+2 dxmti 
which shows that L7(x) _ by 
U(x) = — =i, (x), (n2m) (44.1) 

is a saad af the il Sats equation 


xo ὭΡΗ, Y +(m+1- 2 +(n—m)y = 0. (44,2) 
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The polynomial L(x) defined by equation (44.1) is called 
the associated Laguerre polynomial. It follows from 
equation (42.2) that it may be presented as a series by 
the equation 


C(x) = re ,Fy\(—n+m; m+1;x) (n2m). (44.3) 
Similarly equation (42.4) leads to the formula 

L(x) = ὡς ἷν τ᾿ wren} (44.4) 

The simplest associated Laguerre polynomials are: 

Li(x) = -1, 
L\(x) = —4+2x, L3(x) Ξ 2, 
Li(x) = —18+18x—3x*, L2(x) = 18—6x, L3}(x)= —6, 
Li(x) = —96 + 144x—48x?7+4x3, L3(x) = 144—96x+12x? 


Li(x) = —964+24x, L3(x) = 24. 


The definition (44.1) for the associated Laguerre poly- 
nomial is the one usually taken in applied mathematics. 
In pure mathematics the function 


L(x) = MEM! Pn; mti;x) (44.5) 
πη] 
which is a solution of the differential equation 
d*y | dx . 
x —> +(m+1-—x) — +ny =0 
ΖΞ ( ) i y 
is often taken as the definition of the associated Laguerre 
polynomial + so that care must be taken in reading the 


+ See, for instance, E. T. Copson, Functions of a Complex Variable 
(Oxford University Press 1935), p. 269. 
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literature to ensure that the particular convention being 
followed is understood. 

It is readily shown from equation (42.1) which defines 
the generating function for Laguerre polynomials that the 
associated Laguerre polynomials may be defined by the 
equation 

(—1)"t™ exp (- =) =(i-1"*' } =e) (". (44.6) 
1- n=m ἢ! 

This identity can then be used to derive recurrence relations 

for the associated Laguerre polynomials similar to those of 

equations (42.8) and (42.9) (cf. examples 6 (ii), (iii) below). 

The Laguerre functions R,,,(x) are defined by the equation 


R(x) = ex Last (x) (n21+1). (44.7) 


If in equation (44.2) we replace m by 2/+1, m by n+1 
and y by e**x~'R we see that the function R,,(x) is a 
solution of the ordinary linear differential equation 

aR .2dR i n , Wl+1) 


——_ SEE ...........ΨὍ0. τ 


dx* x dx 


* ee x? 


In most physical problems in which this equation arises 
it is known that / is an integer. By reasoning similar to 
that outlined in §43 it can readily be shown that the 
equation (44.8) possesses a solution which is finite at 
x = 0 and tends to zero as x--0o if, and only if, the para- 
meter 1 which occurs in it assumes integral values. When 
this does occur the solution is AR,,x) where the function 
Κι, (ΑἹ is given by equation (44.7) and A is an arbitrary 
constant. 

We shall now evaluate the integral 


I, = ἢ x7{Rui(x)}72dx 


0 
which arises in wave mechanics. From equation (44.6) we 


Σ =0. (44.8) 
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have the identity 


τ. Lee OLY) nat +1 
Otis Pt (n+1)\(n'+D! ie 


Y—t L—T) 945 2141 
δὼ Ὁ: 174204 —¢)7!*2 t τ 


from which it follows that 


oo oo 


πο ἜΠ 
εἶν π' ἫΝ (n+I)\(n'+D! = 


χ t2 21} (x) P!*1x)dx 


εἰ (it ars δ 2126 τα τ χί XT | 
(—p? 212)? |, οὐ. πῶς ree er dx 


This last integration is elementary and gives for the expres- 
sion on the right 


(21: 2)}{π|2 5 (.--ἡᾷ -- τὴ 
(1 —{t 21+3 


and by means of the binomial theorem we may expand 
this function in the form 


(i-t—t-t) ¥ (21+ Fi (τ) Ὅτε, 
r=0 rt 


Now the coefficient of (¢t)"*' in this expansion is 
2n{(n+1)!} 
(n—I—1)! 

and this is equal to J,)/{(n+1)!}*. Hence 


es 2, 2n{(n+1)!}* ἘΞ 
᾿ x*{R(x)}*dx = f-lepr ; (44.9) 


me -᾿ --- ---- Ὁ --- 
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ἢ 45. The wave functions for the hydrogen atom. We 
shall conclude this chapter by discussing the motion of a 
single electron of mass m and charge —e in the Coulomb 
field of force with potential 

Ze’ 


Vir ΠΕΣ «ταν ὙΠῸ 
r 
due to a nucleus of charge Ze. In a first approximation we 
may treat the mass of the nucleus as infinite and in this 


case the wave function w of the system is governed by the 
Schrédinger equation 


vay 4 am (ws 5: ψ-ὸ (45.1) 


and the conditions: 
(i) ψίγ, θ, φ- 2π) aa Wr, 0, ?) for all r, 8, φ; 

(ii) y& must be bounded in the range 0<@<7 for all 
r, 93 

(iii) ~-0 as r> 00; 

(iv) w remains finite as r0; 

(v) yw is normalised to unity, i.e. | Ψ ?dz = 1 where 
the integral is taken throughout the whole of space. 


W is the total energy of the system. 
The equation (45.1) may be solved by setting 


Ψ = Κι) Θ(θ)Φ(φ), (45.2) 
where, by the method of separation of variables we have 
ve +u’® = 0 (45.3) 
dg? ; 
1 ἃ d© 
— — θ — K(i+1 
τ 9 ἢ (sn iO) + εἰ (I+ aa ©=0, (45.4) 
M 
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1d(.dR\ . (8n2m = M+) 
r? “(1 Ἢ εἾ h? ( Ῥ + (45.5) 


To satisfy condition (1) we must pers as a solution of 
(45.3) a function Φίφ) such that ®(¢+2z) = Φί(φ) for 
all φ. Thus u occurring in equation (45.3) must be an 
integer and a convenient solution will be 

© = Ae (45.6) 
where A is an arbitrary constant. Equation (45.5) is the 
well-known equation of which the solution is the associated 
Legendre polynomial Pi(cos@). If / is integral and 
I=|u|, then Pi(cos@) is the only solution which is 
bounded in the range 0<@Szx and is therefore the only 
one leading to a wave function ~ which satisfies condition 
(ii) above; if / is not an integer no bounded solution exists. 

To solve equation (45.5) we write 


2_ _ 8x°mW 4n*mZe? 
Cc=- ---- -ς-τ v= —— 
h? ho 
and change the independent variable from r to x where 
x = 2ar. We then find that 


2 
Ry 2 ἢ. MADR 0 (458) 
ax? x dx 4 x x 


where, in order that conditions (iii) and (iv) should be 
satisfied, R must be such that R-0, as x-oo and as 
x0. From the arguments of ὃ 44 it follows that this 
is possible only if / is a positive integer and only if v is 
an integer, m say, which is greater than /+1. When this 
is so we may write the solution of equation (45.8) in the 
form (44.7) so that the solution of (45.5) is proportional 
to R,,(2ar). Now by the second of the equations (45.7) 
we have 


(45.7) 


—-2n?Ze4 
W= — —— (45.9) 
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for the possible values of the total energy W. The wave 
functions corresponding to the value of energy given by 
the integer ἢ are 

Walt, 9, Φ) = CyyRyy(2ar)Pi(cos O)e™* (45.10) 
where C,), is a constant determined by the condition (v). 
In polar coordinates dr = r* sin @ drd@ dd so that this 
condition gives 


1 = C3, | ee sin 0{P%(cos 6)}2d0 | r?R2(Qur)dr 
ῃ 0 


0 
so that it follows from equations (21.20) and (44.9) that 


= on (—D-wiin-I- 14 NF 
Cu (2x) Sxu(I-+u)!{(n+1!}° ᾿ (45.11) 


Now from equations (45.7) and (45.9) we find that 
α = 4n?Ze*m|(h7n) so that if we introduce the Bohr radius 
a by means of the equation 
h? 
a= ------ 45.12 

4n* me? ( ) 
we find that α = Z/an. Introducing this result into (45.11) 
and substituting the value obtained for C,,,, into equation 
(45.10) we find that 


Writs 9, φ) 


(a) (+a) {n+} 8a Ryi(2ar)P; (cos θ)ε 
(45.13) 


with uS/Sn—1 are the wave-functions corresponding to 
the energy (45.9) of the hydrogen atom. 

If we write Wy = —2n?Ze*m/h* then corresponding 
to the energy W, we have the wave function 


ΕΖ. τῇ 
Wrool!, 0, φ) ean YE (2) é μ (o = Zr/a), 
J/n\a 
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and to the energy level W,/4 we have the three wave 
functions 


Vro0l?s 9, ¢) = — (2) (2—g)e-#, 


4. [Ὡπ)λα 
' 2 
Worol?, 0, Ψ) = <a () oe” *¢ cos 6, 
Wari(t, 0, 6) = ses) oe * sin Oe*'?, 
8./(2)\ a) ᾿ 


From the last of these functions we can construct two 
functions 
Ἢ 
- (2) ge~# sing ¢. 
4./2n \.a cos 
Similarly corresponding to the energy W,/9 we have the 


six wave functions 


δ δὸς τῆς (2) (27-180 +29%)e¥#, 


r, 0, ) = 2 (2) τ 
Wsi0(r, θ, d) δι π - (6—@)oe ** cos 0, 


+ 
ψοικίν, 8, φ) = ae (2) (6—g)oe™* sin be*", 


4 
Ws20(r, 9, φΦ) = ao (?) oe *%(3 cos? θ-- 1), 


a —1_(2¥ 020-80 53 +i 
Waa1(7, 9, 9) ns) o-e” ** sin 8 cos θε "ἢ, 


: 4+ 
Wa22(r, 8, 6) = aa?) οε΄ ὃ sin? Oe*?", 


162,/n\a 
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In wave mechanics | ¥(r, 0, ¢) [5 dt represents the 
probability that the electron whose wave-function is 
Wr, 8, Φ) is to be found in a small volume dr centred at 
the point whose polar coordinates are (r, 0, φ). To make 
the total probability unity we must have 


Ι Wr, θ, φ)} de -- 1, 


where the integral is taken throughout the whole space. 
Since, in polar coordinates, dr = r* sin 6 drd θ ἐφ it follows 
that the probability that the electron is at a distance 
i r—tdr and r+4dr from the nucleus is ¢(r)dr 
where 


o(r) =r? i. do Ι΄ sin 0 dO .| (r, 0, φ) |?. 
0 0 


Hence for an electron in the state defined by quantum 
numbers n, /, u (equation (45.13) above) we have 


= (22)°_@=I-1)! 
br) = (2) oP (RelQa)) 


The mean value of any function f which depends on r 
alone is given by 


a [λόφον 
0 


that is, by the formula 


. (22 (α--1--}} [ὁ β 
᾿ = aed | 0 POY Rar) yar. 


For examples of the use of this formula see example 17, 
p. 178. 
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Examples V 
(1) Prove that if m<n 
my! 
5 (AAD) = Maal 
(2) If 
K(x, = _Y Val WOW 
where the w,(x) are the orthonormal set of Hermite functions 
defined by the relation 
W(x) + 27 (nl) ἐπ  Ὑ,0), 
prove that 


Ode? +2ax=2? —4y? + 2Aty— Ate? 
e * K(x, ys t)dx gi . = 
-ο 


Assuming that K(x, y, ἢ is of the form 
A exp (Bx? + Cxy+ Dy’) 
where A, B, C and D are functions of ¢ prove that 
1 [Ξ γι--ἀξ- ἐγ Ὁ Ἢ 
fre hen Or rae 
J{x(1—-1)*)} 2(1--ἢ 
(3) Show that 


᾿ phew Δί ως Οὐδ 
ὦ ΡΣ τ" ja—4% ἘΠῚ 


(4) The Schrédinger equation for the three-dimensional 
oscillator in cylindrical coordinates (ρ, @, z) is 


K(x, y, ὃ = 


2 
V+ ἘΞ {W—2n?m(v707 +w7z”)}y = 0. 


Show that the solutions yw of this equation such that 
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0 as ὁ and zoo and p is finite at the origin are of 
the form 


C εχρί iu —4079?—48727} f,), (Q)H,(z), 


where C is a constant, /, , |u| are positive integers, 


c= 2z,/(mv/h), β πες 2z,/ (mo/h), and I ἢ (ἔ is ἃ Ιν- 
nomial of degree 2] in ᾧ which satisfies the odd site 


d? 1\d u? 
th (he ool ας Bina 


Show also that the corresponding values of W are given 
by the equation 


W = (2/+|u|+lhv+(n+dho. 
(5) Prove that 


Laie F ἐπ 


m= 0 m\(n—m)! L,,~»(X). 


(6) Prove that 
. ἃ mi — [πῈ} . 
(i) ay 5)" En); 


(ii) Ls s(x) +(x —2n— 1) 24x) + mE") + πα, () = 0: 


(iii) Μα(χ) -- np _ γ(Χ) ἘΠῚ {(x) = 0. 
(7) If mSn prove that 
| © Info“ Hd = (— 15 =I! gant 1g-opmcgy 
Jo n! οὐ 
(8) Prove that if m<n, a>0, 
e~ MIM dx Εἰς (-1).ὰὰ !)?(a—1)"-™ 
0 (n—m)!a"*! 
and deduce that if n=/+1 


[ x't1R (α)άχ = {(n+1)!}? (—1)rtigzt2, 
0 


(Ἢ -- 1 --ὄ ]}ὺ} 


i] 
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(9) Show that 
(n!)?p!_ κα ptt, 
x1 —x)?L"(ax)dx = (—1)?** ο΄ τ L454 1(@)- 
I a {(n+p)}}? "" 
(10) Prove that the function 
x*e7 * Li 4 p(X) 
is a solution of the differential equation 
d*y a—y+idy 1 2p+a+1 ao 
OY, SOE a + fe + ἘΠ΄ - SS =. 
dx? τ χ tx 4 2x 4x? : 


(11) The potential energy for the nuclear motion of a 
diatomic molecule is closely approximated by the Morse 
function | 

V(r) = De~**"—2De™™. 
Show that the spherically symmetrical solutions of the 
Schrédinger equation with this potential are 


Wir) = 3: fexp(— be *)}(2br)k-"ALH=H; (be), 
r 
(O<nsk-3) 


where ὁ = 2n(2mD)/(ah), C,, is ἃ normalisation constant, 
and the corresponding values of the energy are 


oe eh a. eed OF 
(12) Prove that the normalisation constant C, in the 
previous example has the value 


2ba 
Ν N, 


N, = {T(2b—n)} : P(2b—2n+r—1) 


where 


r! 
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(13) Show that in parabolic coordinates €, ἡ, ¢ defined 
by the equations 
x = (ξη)} cos φ, y = (n)*sin φ, z = 4(E—n) 
Schrédinger’s equation for a hydrogen-like atom of nuclear 
charge Ze takes the form 


a/v), δ, av), 1/1, 1δῳ 
δὲ (: 4 τ On ( 4 + " * ἢ rye 
2π2 


+ ἜΣ {W(E+n)+2Ze"}p -- 0. 


Show that this equation has solutions of the form 
W(é, ἢ, ΦῚ Ξε (#4) εἰ" Gt 2aays Ἢ (£) ἽΒ (1) 
n*a na 


na 


where a is the Bohr radius, n = k+/+a+1 and W= 
—2n*me*z?|(n*h’). 

Determine the constant C so that w is normalised to 
unity. 

(14) In the theory of the rotation and vibration spectrum 
of a diatomic molecule there arises the problem of solving 
the Schrédinger equation with potential energy 

B Zé 


μι) Ξ -- - ----- 
©) γ 3 


Show that the energy levels are given by 
W = —2n*mZ*e*|(h’o") 


where o = n+4}+.,/{6+(/+})*} with n, / integers and 
b = 8n?mB/h*, and that the corresponding wave functions 


Cus (27 ΟΝ (Ξ) P¥(cos O)ei™, 
od od 


where a is the Bohr radius, « = 2,/{b+(1+4)7} and C,, is a 
normalisation factor. 


PS 
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(15) Show that the constant C,, occurring in the last 
example has the value 


2 Vf (21+1)n\(1—u)! \" 
loal(n+a+1)} l4n(Qn+«+1)(1+u)! 
(16) If 
Id; n3s)= [ χ ΞΕ. (Χ)}3άχ, 
} 0 


prove that ΚΙ; n; s) + {(n+J)!}* is the coefficient of (tr)"*! 
in the expansion of 


(1 - ἡ} 11 —tyt! > (2l+s+r+2)! ΔῈ 
Ρ r! 


Deduce that 
hn: _ f+)" 
 G=l=f 
oe _ {(n+D!}° 
ΚΙ; πὶ 1) τ τ ici —2K(1+ 1}. 


(17) Show that in a hydrogen-like atom of nuclear 
charge Ze the average distance of the electron from the 
nucleus, in the state described by quantum numbers /, 


n, is ν᾽ 
πἶα 1f,_ Kl+1) 
: κεῖ ten ᾿ 


Find the average value οὗ 1/r and show that the total 
energy of a hydrogen atom is just one-half of the average 
potential energy. 


APPENDIX 
THE DIRAC DELTA FUNCTION 


ἢ 46. The Dirac delta function. In mathematical 
physics we often encounter functions which have non-zero 
values in very short intervals. For example, an impulsive 
force is envisaged as acting for only a very short interval 
of time. The Dirac delta function, which is used extensively 
in quantum mechanics and classical applied mathematics, 
may be thought of as a generalisation of this concept. 

If we consider the function 


1 
5,(x) = ie te (46.1) 
0, | x |>a, 
then it is readily shown that 
oo 
| 6,(x)dx = 1. (46.2) 


Also, if f(x) is any function which is integrable in the 
interval (—a, a) then, by using the mean value theorem 
of the integral calculus, we see that 


ἰ I(x)6,(x)dx = 1 [(χ)χ = f(@a), | 0 |S1. (46.3) 


We now define 
d(x) = lim δι(χ), (46.4) 


Letting ἃ tend to zero in equations (46.1) and (46.2) we 
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see that d(x) satisfies the relations 
d(x) = 0, if x £0, (46.5) 


᾿ dew 1: (46.6) 


The “function” 6(x), defined by equations (46.5) and 
(46.6), is known as the Dirac delta function. It is unlike 
the functions we normally encounter in mathematics; 
the latter are defined to have a definite value (or values) 
at each point of a certain domain. For this reason Dirac 
has called the delta function an “ improper function ” and 
has emphasised that it may be used in mathematical 
analysis only when no inconsistency can possibly arise 
from its use. The delta function could be dispensed with 
entirely by using a limiting procedure involving ordinary 
functions of the kind 5,(x), but the “ function ” d(x) and 
its “* derivatives ”’ play such a useful role in the formulation 
and solution of boundary value problems in classical 
mathematical physics as well as in quantum mechanics 
that it is important to derive the formal properties of the 
Dirac delta function. It should be emphasised, however, 
that these properties are purely formal. 

First of all it should be observed that the precise variation 
of d(x) in the neighbourhood of the origin is not important 
provided that its oscillations, if it has any, are not too 
violent. For instance, the function 


x)= tin sin (27nx) 
n+ co TX 
satisfies the equations (46.5) and (46.6) and has the same 
an properties as the function defined by equation 
If we let a tend to zero in equation (46.3) we obtain 
the relation 


᾿ f(x)6(x)dx = (0), (46.7) 
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which a simple change of variable transforms to 
oD 
f(x)d(x—a)dx = f(a). (46.8) 


In other words the operation of multiplying f(x) by 6(x—a) 
and integrating over all x is merely equivalent to sub- 
stituting @ for x in the original function. Symbolically 


we may write 
J (x)0(x — a) = f(a)o(x—a) (46.9) 


if we remember that this equation has a meaning only in 
the sense that its two sides give equivalent results when 
used as factors in an integrand. As a special case we have 


xd(x) = 0. (46.10) 
In a similar way we can prove the relations 
6(—x) = d(x), (46.11) 
ni SH, a>; (46.12) 
a 


§(a?—x?) = {6(x—a)+6(x+a)}, a>0. (46.13) 
Let us now consider the interpretation we must put 
upon the “ derivatives” of d(x). If we assume that 6’(x) 
exists and that both it and d(x) can be regarded as ordinary 
functions in the rule for integrating by parts we see that 
a oo 
[7 sere eoas =| fo9800 
- |" repaeodr = 1°. 
Repeating this process we find that 


4 fO)5%x)dx =(—1)7%0). (4614) 
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The statement is often made that the Dirac delta 
function is the derivative of the Heaviside unit function 


H(x) defined by the equations 
1, [[χ»0; 
Η —_ 3 ἢ 
©) ἢ if x<0: 


and it is easy to see on geometrical grounds that there are 
reasons for conjecturing such a relationship. To make it 
precise we note that if, in the definition 7 of the Stieltjes 
integral 


᾿ flce)d F(x), 


we take F(x) to be the Heaviside function H(x), we find 
immediately that, for any integrable function f(x), 


| ” fxdH(x) =f). (46.15) 


Comparing equation (46.15) with equation (46.7) we see 
the relation between H(x) and d(x). It may be seen from 
these equations that d(x) is not a function but a Stieltjes 
measure, and that the use of the Dirac delta function could 
be avoided entirely by a systematic use of Stieltjes in- 
tegration. 


+ The simplest definition of a Stieltjes integral (“ f(x)dF(x) is as 


the limit of approximative sums 2 /(¢,)[F(x,) γ- πὰ ο — the 
x,, the points of subdivision of (a, δ) and €,, lies in the interval (x,., x,). 
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